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On the Smarandache additive sequence 


Fanbei Li 


School of Mathematics and Statistics, Inner Mongolia Finance and Economics College 


Hohhot 010051, Inner Mongolia, P.R.China 


Abstract In this paper, we use the elementary method to study the properties of the Smaran- 
dache additive sequence, and give several interesting identities for some infinite series involving 


the Smarandache additive sequence. 


Keywords Smarandache additive sequence, infinite series, convergent properties. 


81. Introduction and results 


For any positive integer m > 2, let 1 < di < dg < --- < dm are m positive integers, 
the Smarandache multiplicative sequence A,,, is defined as: If dj, dz, ---, dm are the first 
m terms of the sequence A,,, then dy, > dg_1, for k > m+ 1, is the smallest number equal 
to dj) - d5?---d°», where a; > 1 for alli = 1, 2, ---, m. For example, the Smarandache 


multiplicative sequence A, ( generated by digits 2, 3, 5, 7 ) is: 
2, 3, 5, 7, 210, 420, 630, 840, 1050, 1260, 1470, 1680, 1890, 2100, ---. 


In the book “Sequences of Numbers Involved Unsolved Problems” , Professor F.Smarandache 
introduced many sequences, functions and unsolved problems, one of them is the Smarandache 
multiplicative sequence, and he also asked us to study the properties of this sequence. About 
this problem, Ling Li [3] had studied it, and proved the following conclusion: 

For any positive integer m > 2, let 1 < d, < dz <-:- < dm are m positive integers, and 
Am denotes the Smarandache multiplicative sequence generated by dj, dz, ---, dm. Then for 


any real number s > 0, the infinite series 


1 
ae 


ne€Am 
is convergent, and 
1 7 1 “1 
> A-lgathe 
Ss S 
n€Am ue i=1 dj — 1 i=1 di 


In this paper, we define another sequence called the Smarandache additive sequence as 
follows: Let 1 < dy < dg < --- < dy» are m positive integers, we define the Smarandache 
additive sequence D,, as: If di, dz, --- , dm are the first m terms of the sequence D,,, then 
dy > dx_1, for k > m+1, is the smallest number equal to a, - dj + ag: dg+-----: + Am -dm, 
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where a; > 1 for alli =1, 2, --- , m. For example, the Smarandache additive sequence Dp ( 
generated by digits 3, 5 ) is: 


5.5, & 11, 13) 14, 16,17; 18, 19; 20, 21, 22, 23, 24, 25, 96, 27, 28, +, 


About the properties of this sequence, it seems that none had studied it yet, at least we 
have not seen any related papers before. The problem is interesting, because it has the close 
relationship with the positive integer solutions of the indefinite equation 


dy -%,+dg-%9+d3-43+---+dm_—1°Lm_-1 + dm + tym =N. 


The main purpose of this paper is using the elementary method to study the convergent proper- 
ties of some infinite series involving the Smarandache additive sequence, and get some interesting 
identities. For convenience, we use the symbol a,,(n) denotes the n-th term of the Smarandache 
additive sequence D,,. In this paper, we shall prove the following: 

Theorem. For any positive integer m > 2, let 1 < dj < dg <-:: < dm are m positive 
integers, and D,,, denotes the Smarandache additive sequence generated by d;, do, ---, dm. 
Then for any real number s < 1, the infinite series 





(1) 


is divergent; For any real number s > 1, the series (1) is convergent, and 
= 4 


De =9+ YY ag-bLe 


neEDm n= am(n)<M  ™ n<M 








fan 
=) 

3% 
= 


where ¢(s) is the Riemann zeta-function, and M is a computable positive integer. 
From this Theorem we may immediately deduce the following three corollaries: 
Corollary 1. Let Dg be the Smarandache additive sequence generated by 2 and 38, then 


we have the identity 
2 


3 1 =3 1 157 
ba er ee az(n) 6 144° 





Corollary 2. Let D2 be the Smarandache additive sequence generated by 4 and 6, then 
we have the identity 





Se 1 x? 157 
n2 4+ aq2(n) 24 576 


Corollary 3. Let D3 be the Smarandache additive sequence generated by 3, 4 and 5, then 
we have the identity 





ee eee ei ee ee 
n2- £~a2(n) 6 4 36 49 64 81 100 121 169 196° 
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§2. Proof of the theorem 


To prove our Theorem, we need a simple Lemma which stated as follows: 
Lemma. Let a and b are two positive integers with (a, b) = 1, then for any positive 
integer n > ab —a — 6, the indefinite equation 


ax + by =n 


must has nonnegative integer solution. 
Proof. See Theorem 19 of Chapter one in reference [3]. 
Now we use this Lemma to prove our Theorem. First note that for any positive integer 
n >m, we have 
AQm(n) = ay + dy + Ag+ dg+-+-+Qm + dm, 
where a; > 1 for alli=1, 2, 3, ---, m. 


From Lemma and its generalization we know that there exist a positive integer M such 


that for any integer n > M, the indefinite equation 
dy + % + dg-%q+dg-%3 +-+-+dm—1+Lm—-1+dm + Lm = 7 


must have positive integer solution. So from the definition of the sequence D,, we know that 
for any integer n > M, there must exist one term ay,(k) in sequence D,, such that n = a,,(k). 
Therefore, 








1 oil 1 1 
ee TAS ~ Dy a(n) = as (n) 


s 
n€Dm nim Am (n)>M oo ™ Am (n)<M™ 


= 1 1 1 
= ee ee ea > as, (n) 








n=1 n<M am(n)<M ™ 
1 1 
= s)- 0 “a Ds as,(n)’ 
n<M Am (n)<M 


co 
where ¢(s) is the Riemann zeta-function, and the infinite series S- —, is convergent if s > 1; 
n 


n=1 


It is divergent if s < 1. This completes the proof of Theorem. 

Now we prove our Corollary 1. Taking a = 2, b = 3, note that a-b—a—b = 2-3-—2-3=1, 
so from Lemma we know that for all integer n > 2, the indefinite equation 2x + 3y = n has 
nonnegative integer solution. Therefore, for all integer n > 7, the indefinite equation 2x+3y = n 
has positive integer solution. So M = 7, from our Theorem we may immediately deduce that 


“1 1 1 
Lam = + 2 aa de 


n=1 Am(n)<7  ™ 








2 1 1 1 1 
= 2454-45-57 
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2 
where we have used the identity ¢(2) = 7. This proves Corollary 1. 
To prove Corollary 2, note that the greatest common divisor of 4 and 6 is (4, 6) = 2, so 


from Corollary 1 we have the identity 


a | 1 1.2. a 1 
Bec = [@eteted oa} 


n=1 








mw? 157 
24 576 
Now we prove Corollary 3. Note that for all n > 15, the indefinite equation 


30+ 4y+5z=n 


has positive integer solution. So we can take M = 15. Therefore, 


Lao = + Yaw Law 


a3(n) 











n=1 a3(n)<15 n<15 
_. 1 1 1 1 1 1 1 1 1 1 
6 4 36 49 64 81 100 121 169 += 196 
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Contributions to classical differential 
geometry of the curves in FE? 


Melih Turgut and Suha Yilmaz 


Dokuz Eylul University, Buca Educational Faculty 


Department of Mathematics, 35160, Buca-Izmir, TURKEY 
E-mail: melih.turgut@gmail.com 


Abstract In this paper, a system of differential equation whose solution gives the components 
on the Frenet axis in E® is established by means of Frenet equations. In view of some 
special solutions of mentioned system, position vector of rectifying curves, osculating curves 
with constant first curvature, normal curves and special cases are presented. Moreover, a 


characterization for inclined spherical curves is expressed. 


Keywords Euclidean space, Frenet frame. 


81. Introduction 


It is safe to report that the many important results in the theory of the curves in E? were 
initiated by G. Monge; and G. Darboux pionnered the moving frame idea. Thereafter, F. Frenet 
defined his moving frame and his special equations which play important role in mechanics and 
kinematics as well as in differential geometry (for more details see [2]). 


In the Euclidean space E®, it is well-known that to each unit speed curve with at least four 
continuous derivatives, one can associate three mutually orthogonal unit vector fields T, N and 
B are respectively, the tangent, the principal normal and the binormal vector fields[3]. At each 
of the curve, the planes spanned by {T, N}, {T, B} and {N, B} are known respectively as the 
osculating plane, the rectifying plane and the normal plane[5]. The curves y: I C R — E® for 
which the position vector y always lie in their rectifying plane, are simplicity called rectifying 
curves [1]. Similarly, the curves for which the position vector y always lie in their osculating 
plane, are for simplicity called osculating curves and the curves for which the position vector 
always lie in their normal plane, are for simplicity called normal curves [4]. 


In this paper, position vectors of some special curves in EL? are investigated. First, a system 
of differential equation is established. In view of some special solutions, characterizations of 


some special curves are presented. 
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§2. Preliminaries 


To meet the requirements in the next sections, here, the basic elements of the theory of 
curves in the space EF? are briefly presented. (A more complete elementary treatment can be 
found in [5].) 

The Euclidean 3-space E° provided with the standard flat metric given by 


(,) = dx? + dx3 + dz3, 


where (21, 22,73) is arectangular coordinate system of E*. Recall that, the norm of an arbitrary 
vector a € E? is given by ||al| = \/(a,a). ¢ is called an unit speed curve if velocity vector v of 
y satisfies ||v|| = 1. For vectors v,w € E® it is said to be orthogonal if and only if (v, w) = 0. 
Let 0 = 0(s) be a curve in E®. If the tangent vector field of this curve forms a constant angle 
with a constant vector field U, then this curve is called an inclined curve. The sphere of radius 
r > 0 and with center in the origin in the space E® is defined by 


Ss = {p = (p1, P2, P3) € E : (p, p) = rt. 


Denote by {T,N,B} the moving Frenet frame along the curve y in the space E3. For an 
arbitrary curve y with first and second curvature, « and 7 in the space FE, the following Frenet 


formulae are given in [3], [5] 


fb 0 K 0 T 
N'’ |=] -K O T N |, (1) 
B’ 0 —r 0 B 


where 


Here, first and second curvature are defined by « = K(s) = ||Z’(s)|| and 7(s) = — (N, B’). 

It is well-known that for a unit speed curve with non vanishing curvatures following propo- 
sition holds [3]. 

Proposition 2.1. Let y = y(s) be a regular curve with curvatures « and 7, then; 
i. It lies on a sphere ($7) if and only if 





ii. It is an inclined curve if and only if 


* = constant. (3) 
: 
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§3. Position vectors of some special curves in E* 
Let y = y(s) be a unit speed curve in E3. We can write this curve respect to Frenet frame 
{T, N, B} as 
p= y(s) => mT + m2N + m3B. (4) 


Differentiating (4) and considering Frenet equations, we have a system of differential equation 
as follow: 


m),— max —-1=0 
m,+mik—m3T=0 ?- (5) 
ms + mT =0 
This system’s general solution have not been found. Due to this, we give some special values 


to the components. 
Case 3.1. m, = c;=constant. From (5); we have second component 


m2 =— (6) 


And (5)2 gives us third component as 


gee” (*) ta%. (7) 


aT K 


mg and mz satisfie (5)3, if we rewrite it, we have following differential equation: 


(2) ] += GY-E-0 " 


Considering (8) and Proposition 2.1., suffice it to say that following theorem holds: 





Theorem 3.1. Let y = y(s) be a unit speed curve and first component of position vector 
on Frenet axis be constant in E°: 
i. If first component is zero, then y is a spherical curve (lies on $7). 
ii. If y is an inclined curve, then y lies on 9%. Therefore, y is an inclined spherical curve. 
iii. All normal curves in the space E® are spherical curves. 
iv. Position vector all of normal curves in E%, can be written as follow: 


p=-iw-1(2\ (9) 


K T K 


Here, it is not difficult to see that statement ii is a characterization for inclined spherical 
curves in FE. 
Case 3.2. m2 = cg = constant # 0. From equation (5), we obtain 


/ 
mi, —1 





From (5)2, we write 
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Now, let us suppose y is an inclined curve. In this case, If (10) and (11) are used in (5)3, we 


have a differential equation as follow: 
,/K T 
Se es hy 12 
m(=+=)-7=9 (12) 


From solution of (12), we have the components 


TS 
2 
—K 
er re ay) 
and : 
Ke TS 
oat 1 
™3 aa +1), ( 5) 
where 1 € R. 


Corollary 3.2. Let y = y(s) be a unit speed curve and second component of position 
vector on Frenet axis be constant in E%. Then position vector of y can be formed by the 
equations (13), (14) and (15). 

Case 3.2.1. m2 = 0. From (5), we have, respectively, 


m=stl (16) 
and 
kK 
m3 = —(s+1) = constant, (17) 
= 


where /€ R. 
Considering (16) and (17), we give following theorem. 
Theorem 3.3. Let y = (s) be a unit speed rectifying curve in E%. Then; 


i. There is a relation among curvatures as 
K 
—(s +1) = constant. 
= 
ii. Position vector of y can be written as follow: 
K 
p= 9(s)= (st OT + — (s+) B, (18) 


where / € R. 

Remark 3.4. The case m3 = constant ¥ 0 is similar to case 3.2. 

Case 3.3. m3 = 0. Thus 7 = 0. And here, let us suppose k = constant. Using these in 
the system of differential equation, we have second order homogenous differential equation with 
constant coefficient as follow: 


my + monk? +K=0. (19) 
Solution of (19) is elementary. Then, we write the components 


my, = Ao sinks — A; cos Ks (20) 
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and 
. 1 
mg = A,cosKs+ Agsinks + —, (21) 
K 


where A; and A» are real numbers. 
Corollary 3.5. Let » = y(s) be a unit speed osculating curve with constant first curvature 
in E3. The position vector of y can be composed with the equations (20) and (21). 
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Almost everywhere stability of the 
third-ordered nonlinear system 


Kaiming Wang 
School of Science, Chang’an University, Xi’an, Shaanxi, P.R.China 


E-mail: kaimingwang@yahoo.com.cn 


Abstract In order to study the stability of the third-ordered nonlinear system, density function 
is constructed, and an inequality condition is obtained. When the considered system satisfies the 
inequality for almost all initial value, its equilibrium is asymptotically stable almost everywhere. 


Keywords Amost everywhere stability, density function, nonlinear system. 


81. Introduction 


Stability is a very important factor in the research of high-ordered nonlinear systems. The 


techniques currently available in the literature are mostly based on Lyapunov’s direct method 


[1], contraction mapping method [2], nonlinear measure method [3,4], nonlinear spectral radius 


method [5] and so on. All these methods make the research of stable property everywhere. And 


if the equilibrium of the considered system is asymptotically stable, there exists a neighborhood 


of the equilibrium which is called attraction region. In fact, many real nonlinear systems do not 


need its equilibrium acting as an attracting actor attracting any points of its neighborhood, and 


it only needs attract “enough many” points in this neighborhood. That is, the equilibrium is 


asymptotically stable almost everywhere (the equilibrium attract all points of this neighborhood 


except a set with zero Lebesgue measure). In 2001, Anders Ranzter [6] proposed the concept 


of a “dual” to Lyapunov’s stability theorem, and he called the stability derived by this way as 


stability almost everywhere. In control theory, most control systems are in high order, such as 


robotic multi-fingered grasping system [1], neural networks and so on, and up to their research 


background, it’s unnecessary or impossible for the equilibrium to attract every points in its 


neighborhood even it’s very small. So it’s valuable to study the almost everywhere stability of 


high-ordered systems. 


In this paper, the almost everywhere stability of the third-ordered system is studied. Den- 


sity functions are constructed as a tool to analyze the stability of the given system. 


For convenience, the following notations were used in this paper: 
(6) fe) fe) 

(1) Vp = gradp = (# Pe, es £-), ps R= R: 

(2) Ve F=divF = oe | oe Alea pea, F:R° = R’,; 


(3) “e” represents inner product. 
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§2. Main result and its proof 


Consider a third-ordered nonlinear system 


, 


est w(a,x ja + f(a,2 ) = 0, f(0,0) =0, (1) 


or its equivalent system 


/ 


at y 
] = z = F(a, y, 2). (2) 
z —Y(z,y)z _ f(a, y) 


It is assumed that solutions of (1) exist and are unique. 

Theoretically, this is a very useful nonlinear system since (1) is a rather general third- 
ordered nonlinear differential equation. For example, many third-ordered differential equations 
which have been discussed in [8] are special cases of (1). However, most available research all 
focused on the study of everywhere stability and it is not necessary or impossible for many real 
nonlinear equations. 

In this section, the almost everywhere stability of third-ordered system is studied. Density 
functions are constructed as a tool to analyze the stability of the given system, and the condition 
was obtained that the given nonlinear system is asymptotically stable almost everywhere. That 
is, we will prove the following main result: 

Theorem. If there exists a positive real number a, such that the following inequalities 

(1) af(z,y)>0, yf(x,y) >0; 

(2) U(a,y) > 0; 

(3) fo (f(@,0) + yea, v)u + yfe(a, 0) + 20(a, y)) dut 

+o-hy(o,y) (Je fu, O)du+ [2 (ow(a,v) + f(e,v))dv)) 
< W(a,9)2(1 — 507?) 
hold for almost all x, then the equilibrium of the nonlinear system (2) is asymptotically stable 
almost everywhere. 

First, we need a basic Lemma that is necessary in proof of Theorem. 

Lemma. Given the equation «(t) = f(a(t)), where f € C1(R”, R”) and f(0) = 0, suppose 
there exists a non-negative scalar function p € C'(R” \ {0}, R) such that p(x) f(ax)/|a| is inte- 
grable on {a € R” : |z| > 1} and [Ve (fp)|(x) > 0 for almost all x. Then, for almost all initial 
states (0) the trajectory x(t) exists for t € (0,00) and tends to zero as t > co. Moreover, 
if the equilibrium x = 0 is stable, then the conclusion remains valid even if p takes negative 
values. 

Proof. ( See reference [6] ). 

By this Lemma, now we shall complete the proof of Theorem. 

The proof of Theorem. 

Let £ = (x,y,z)' and 


xL y y 
va@= | Flu,0)du+ | vb(a,o)do+ f (an, v)dv + 52%, 
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then V(#) > 0 and when & # 0 we have V(%) 4 0. That is, V(Z) is positive definite. Let 
p=V~°, then p € C!(R3\{0}, R) and p(Z)F (2) /|Z| is integrable on {% € R? : |z| > 1}. 
For divF = —w(a,y) and 


[V e (Fp)|(Z) =V~°""(4)(VV e F—aVV e F)(z), 


thus the condition [V e (F'p)|(%) > 0 in Lemma will be satisfied if and only if (VVeF-—aVVe 
F)(%) > 0. From 


VV eF 
V,a' + Vyy! + Vz2' 


y y 
= Fe Oy ey | hae, vod + y | le aaa abcaae Laleaye, 


and 


a lVVeF , 
= —a!w(z, ,0)du— a w(x, ,v)d 
a we» | f(u,0)du— a ven [ vo(e,edv 


y 
—aty(ey) f° f(a,v)av— 50-1 v(e. 9), 


then by Lemma, we know that if the inequality 
y 
i (F(x,0) + ye (x, v)u + yf (a, v) + 20(a,y)) dv 
£ y 
or titeg) ( [sw oaus ov(e,n) + f(2,0)}ae)) 
0 0 


< W(2,y)220 - 507%) 


holds, the equilibrium of the considered system is asymptotically stable almost everywhere. 

This completes the proof of Theorem. 

Simplifying the condition of the inequality in Theorem, we may immediately derive the 
following two corollaries. 

Corollary 1. If there exists a positive real number a > - such that two functions f(z, y) 
and 7(x, y) satisfy the inequalities 

(1) af(e,y)>0, yf(e,y) > 0; 

(2) o(a,y) > 0; 

(3) fo (f(@,0) + yva(x, v)v + yfa(x,v) + U(x, y)z)dv < 0; 

(4) i f(u, 0)du + So (ev (2, v) + f(x,v))du < 27(a- 5) 
for almost all x, then the equilibrium of the nonlinear system (2) is asymptotically stable almost 
everywhere. 

Corollary 2. If there exists a positive real number a > ; such that two functions f(z, y) 
and w(a,y) satisfy the inequalities 

(1) af(z,y)>0, yf(x,y) >0; 

(2) (x,y) > 0; 
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(3) yf(z,0) + yev(a,y) + y? (ala, OE + fo(z,m)) < 0, VE. € [0,y]; 
(4) xf(s, 0) Ax y(hy(a, h) Te f(h, 0)) < 2(a _ 5)Vs € (0, x], he (0, y] 
for almost all x, then the equilibrium of the nonlinear system (2) is asymptotically stable almost 





everywhere. 

Since the system (1) and system (2) are equivalent systems, and the equilibrium «* = 
0 of system (1) answers to the equilibrium <* = 0 of system (2), their stability properties 
are also equivalent. Thus the equilibrium of the third-ordered nonlinear system (1) is almost 
asymptotically stable when the inequality condition of Theorem is satisfied. 
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Abstract For any fixed positive integer m, if there exists a positive integer n such that 
n =m-d(n), then m is called an index of beauty, where d(n) is the Dirichlet divisor function. 
In this paper, we shall prove that there exist infinite positive integers such that each of them 


is not an index of beauty. 
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§1. Introduction and result 


For any positive integer n, the famous Dirichlet divisor function d(n) is defined as the 
number of all distinct divisors of n. If n = p{" - p5?---p?* be the prime power factorization of 
n, then from the definition and properties of d(n) we can easily get 


d(n) = (a1 +1) - (ag +1)--- (ax +1). (1) 


About the deeply properties of d(n), many people had studied it, and obtained a series results, 
see references [2], [3] and [4]. In reference [5], Murthy introduced an index of beauty involving 
function d(n) as follows: For any positive integer m, if there exists a positive integer n such 


that 
n 


a(n)’ 


then m is called an index of beauty. At the same time, he also proposed the following conjecture: 


m= 


Conjecture. Every positive integer is an index of beauty. 

Maohua Le [6] gave a counter-example, and proved that 64 is not an index of beauty. We 
think that the conclusion in [6] can be generalization. This paper as a note of [6], we shall use 
the elementary method to prove the following general conclusion: 

Theorem. There exists a set A including infinite positive integers such that each of n € A 
is not an index of beauty. 


§2. Proof of the theorem 


In this section, we shall prove our Theorem directly. For any prime p > 5, if we taking 
m = p?-', then we can prove that m is not an index of beauty. In fact, if m is an index of 
beauty, then there exists a positive integer n such that n = p?~!-d(n). From this identity we 
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can deduce that p?~! | n. Let n = p*-b and (p, 6) = 1, then from n = p?~!- d(n) and (1) we 
have the identity 


pv?! b= (a +1) -d(d). (2) 


From (2) we may immediately deduce that a > p—1. If a= p—1, then (2) become b = p-d(b), 
this contradiction with (b, p) = 1. If a = p, then (2) become p- b = (p+ 1) - d(b), or 


b= (1 a *) - d(b). (3) 


It is clear that (3) is not possible if prime p > 5. In fact ifb = 1, 2, 3, 4, 5, 6, 7, then it is 
easily to check that (3) does not hold. If b > 8, then from (1) and the properties of d(b) we can 


1 

deduce that b > ; -d(b) > (: + —]-d(b). So (3) does not also hold. So if a = p, then (2) is 
Pp 

not possible. If a >p+1, let a—p=k>1, the (2) become 


pet! .b = (p+k+1)-d(b). (4) 


Note that b > d(b) for all positive integer b, so formula (4) is not possible, since p*t! > p+k+1 
for all prime p > 5 and k > 1. Since there are infinite prime p > 5, and each p?~! is not an 
index of beauty, so there exists a set A including infinite positive integers such that each of 
n € A is not an index of beauty. This completes the proof of Theorem. 
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Abstract For any positive integer n, let ¢(n) denotes the Euler function. That is, @(n) 
denotes the number of all positive integers not exceeding n which are relatively prime to 
n. g(n) denotes the smallest positive integer m such that n | ¢(n)-m. The main purpose 


of this paper is using the elementary method to prove that for any positive integer group 


k k 
(m1, M2, «+: , Me), we have the inequality g (11 ms) < Il g(m,), and also for any positive 
i=1 i=l 
integer k, there exist infinite group positive integers (mi, m2, --- , mx) satisfying the equation 
k k 
a (Sem) =¥- atmo. 
i=1 i=l 


Keywords Euler function ¢(n), three prime theorem, elementary method. 


§1. Introduction and results 


For any positive integer n, the famous Euler function ¢(n) is defined as the number of 
all positive integers not exceeding n which are relatively prime to n. Now we define a new 
arithmetical function g(n) as follows: g(n) denotes the smallest positive integer m such that 
n | @(n)-m. That is, g(n) = min{m: n| ¢(n)-m}. For example, the fist few values of g(n) 
are: g(1) = 1, 9(2) = 2, 9(3) = 3, g(4) = 2, 9(5) = 5, 9(6) = 3, 9(7) = 7, 9(8) = 2, 9(9) =8, 
g(10) = 5, g(11) = 11, g(12) = 3, 9(13) = 13, 9(14) =7, --- 

About the elementary properties of g(n), it seems that none had studied it before. From the 
definition of g(n) we can easily deduce that for any prime p, g(p%) = pand P(n) < g(n) < [[-. 


pin 
where P(n) denotes the largest prime divisor of n. If n > 1 and n 4 2°, then g(n) must be an 


odd number. g(n) is also a square-free number. 

The main purpose of this paper is to study the elementary properties of g(n), and prove 
the following two conclusions: 

Theorem 1. For any positive integer k and any positive integer group (m1, m2, --: , Mz), 
we have the inequality 

k 
g (11 ms) < [[ 9(mi). 
i=1 i=1 


Theorem 2. For any positive integer k > 3, there exist infinite group positive integers 
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(m1, M2, +++, Mz) such that the equation 


k k 
g (>: ms) = de g(mi). 


Whether there exist infinite group positive integers (m1, m2) satisfying the equation 


g(my + m2) = g(m1) + g(m2) (1) 


is an unsolved problem. If m, > 1 and mz > 1 satisfying the equation (1), then one of m; and 
mz must be a power of 2. 
It is an interesting problem to study the mean value properties of g(n). For any real 


number x > 1, we conjecture that 


S g(n) ~C-2?, 


n<u 


where C is a computable constant. 


§2. Proof of the theorems 


In this section, we shall use the elementary method to prove our Theorems. First we prove 
Theorem 1. For any positive integer m, since g(m) = g(m). So Theorem 1 holds if & = 1. 
Now we assume k = 2. For any positive integers m and n, if (m, n) = 1, then we have 
o(mn) = o(m)¢(n). This time let g(mn) = t, g(m) = u, g(n) = v. From the definition of g(n) 
we have mn | d(mn)-t, m | d(m)-u, n | d(n)-v, so we have 


mn | o(m)- d(n)- uv = d(mn) - uv. 


By the definition of t we have t < wv. Therefore, g(mn) < g(m)- g(n). 
If (m, n) > 1, without loss of generality we can assume that m = p® and n = p®, where 
p be a prime, a and @ are two positive integers, then from the definition of g(n) we have 


g(p*p") = p, g(p™) = p, g(p") = p. So g(p*-p*) < g(p*) - g(p"). From the above we know that 
for any positive integers m and n, we have the inequality 


g(mn) < g(m) - g(n). (2) 


Now for any positive integer k > 2 and any positive integer group (m1, m2, --:, Mk), 
from (2) we have 


k k-1 k-1 k 
g (11: =g (m II m) <g (I ms) -g(me) Se < [[ acm). 


This proves theorem 1. 
To prove Theorem 2, we need the famous Vinogradov’s Three Prime Theorem, which was 
stated as follows: 
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Lemma 1. There exists a sufficiently large constant AK > 0 such that each odd integer 
n > K can be written as a sum of three odd primes. That is, n = p, + po + p3, where p; 
(i =1, 2, 3) are three odd primes. 

Proof. (See reference [1]). 

Lemma 2. Let k > 3 be an odd integer, then any sufficiently large odd integer n can be 


written as a sum of k odd primes 
N= pit pat-+:+ De. 


Proof. (See reference [2]). 
Now we use these two Lemmas to prove our Theorem 2. If k = 3, then from Lemma 1 we 


know that for any prime p large enough, it can be written as a sum of three primes: 


P=P1+ p2+ Pps. 


Now for any positive integer k > 4 any prime p large enough, if k be an odd number, then 
by the definition of g(n) we know that g(p) = p. From Lemma 2 we may immediately get 


p= g(p) = g(p1 + po +++: + pe) = pit pet -::+pe=g(pi) +--+ +9(pr). 


If k > 4 be an even number, then k—1 > 3 be an odd number. So if prime p large enough, 
then p — 2 can be written as a sum of & — 1 primes: 


p-2=pitpet:::+ pri. 


This implies 





p=2+pi t+ pote+++ pri 
or 
p=g(p) = g9(2+pit pet: ++ pe_-1) = 2+pitpet-:-+pe_1 = 9(2)+9(p1) +--+ 9(pe-1). 


Since there are infinite primes p, so there exist infinite group positive integers (m1, m2--- , Mr) 


satisyfying the equation 


k k 
g (>: m) ~ de g(mi). 


This completes the proof of Theorem 2. 
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Abstract It is proved that product of interval number and triangular fuzzy number is a trape- 
zoidal fuzzy number. This paper transformed a multiple attribute decision making problem 
with attribute weight of interval number and attribute value of triangular fuzzy number into 
trapezoidal fuzzy number multiple attribute decision making problem, gave TOPSIS to trape- 
zoidal number multiple attribute decision making problem, and analyzed a practical problem 
by this method. 


Keywords Interval number, triangular fuzzy number, trapezoidal fuzzy number, multiple 
attribute decision making, TOPSIS. 


§1. Introduction 


The multiple decision making (MADM) problems mainly solve to rank alternatives with 
multiple attributes among finite alternatives. It is an extremely active topic in discipline re- 
search of modern decision science, systems engineering, management science and operations 
research, etc. Its theory and method are applied to many domains: engineering design, econ- 
omy, management and military affairs, etc. Reference [{1],[2],[3] studied MADM problems that 
attribute weight information is interval number, attribute value is crisp number from various 
angles. Reference [4]-|7] studied MADM problems that attribute weight information is interval 
number, attribute value is also interval number from various angles. At present research results 
about methods for MADM problems with attribute weight of interval number and attribute 
value of triangular fuzzy number are few. This paper focuses on this kind of MADM problems 
mainly. This paper will prove that product of interval number and the triangular fuzzy number 
is a trapezoidal fuzzy number and transform a multiple attribute decision making problem with 
attribute weight of interval number and attribute value of and triangular fuzzy number into 


trapezoidal fuzzy number multiple attribute decision making problem, then gives TOPSIS to 





1 This work is supported by scientific research plan foundation of Xi’an University of Technology(210406) 
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trapezoidal number multiple attribute decision making problem. Finally, an example is used to 
illustrate the proposed method. 


§2. Presentation of the problem 


Definition 1. Let a be interval number, where a = [a~, at] = {ala~ <a <at,a <at}. 
If 0 <a <at, then a is called a positive interval number. 

Definition 2. Fuzzy number p is called a triangular fuzzy number if its membership 
function pip(x) : R — [0,1] is described as 








x l 
m1” m= x € [I,m] 
p(x) = ; x €|m,ul 
m-uU m—U 
0, otherwise 


where x € R, 1 < m < u,l,u stand for the lower and upper bounds of the support of p 
respectively, and m for the most possible value of p. The triangular fuzzy number p can be 
denoted by (l,m, u). When | = m = u, p is a real number. 

Definition 3. Fuzzy number f¢ is called a trapezoidal fuzzy number if its membership 


function juz(z) : R — [0,1] is described as 








0, r<a 
r—a 
; a<xa<b 
b-a 
1G(x) = 1, b<a<ec 
d— 
a c<a<d 
0, d<a2 


where a < b < c < d, a and d stand for the lower and upper bounds of the support of t 
respectively, and closed interval [b,c] for the most possible value of t. The trapezoidal fuzzy 
number t can be denoted by (a,b,c,d). When a = b,c = d, t becomes a ordinary interval 
number. When 6 = c¢, ¢ turns into triangular fuzzy number. Interval number and triangular 
fuzzy number are special cases of trapezoidal fuzzy number. When a = b = c = d, it isa 
ordinary real number. 

The \-cut of a trapezoidal fuzzy number t = (a,b,c, d) is 


ty = [a+ A(b—a),d + A(c—d)], A € (0,1). 


Definition 4./8] Let vector T= (t1, t2,-++ ,tm) and V= (01, 02,°++ ,Um), where trape- 
zoidal fuzzy number t; = (tri, tai, tai, tai), 0; = (Wri, Dai, Vai, V1i), 1 = 1,2,--- ,m. The distance 
between JT and V is defined as 








D(T,V) = ||T-V|| = [(taz — 1a)? + (toe — V2e)? + (tas — U3)? + (bas — vae)?]- 


Ble 


i=1 
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We will use product of interval number and triangular fuzzy number in our discussion, for 
this reason we prove the theorem as follows. 

Theorem. Let a = [a1,Q2] be arbitrary positive interval number, p = (1,m,n) be arbi- 
trary triangular fuzzy number, then the product of a and p is 


a: p= [a1, a9] - (I,m, u) = (ail, arm, agm, agu). 


Namely the product of positive interval number and triangular fuzzy number is trapezoidal 
fuzzy number. 
Proof. For VA € (0, 1], the A-cut of p is 


pa = [I+ Mm—1),u+ Am — w)] 


According to extension principle we have 


ap = LU Mlar,aa][l+(m—1),u+ A(m— u)]}} 
0<A<1 
= U Alayl + A(arm — al), au + A(agm — agu)] 
0<A<1 


= (ayl,aym, agm, au). 


§3. Principle and step 


This paper considers the fuzzy MADM problems that the alternatives are X;, X2, ---, 
X,, and attributes are G,,G2,---,Gm. The fuzzy MADM problems with incomplete weight 
information, w; denotes the weight of attribute G;, w; € [a;,G;], where, 0 < a; < 8; <1, 
j=1,2,---,m, wy +w2+--:+w, = 1, is known to us. Triangular fuzzy number (1;;,mji;, ui;) 
is the value for alternative X; with respect to attribute G;, i = 1,2,---,n, 7 = 1,2,---,m. 
Principle and step for the fuzzy MADM problems that weight of attribute is interval number 
and value of attribute is triangular fuzzy number are as follows: 

Step 1. Construct triangular fuzzy number decision-making matrix A = (Lij)nxm- Where 
triangular fuzzy number 2%; = (lij,mji;,Uij) is the value for alternative X; with respect to 
attribute G;,7=1,2,---,n, 7 =1,2,---,m. 

Step 2. Normalize decision-making matrix A. 

Generally speaking, the attribute values have different physical dimension from each other, 
the decision matrix needs to be normalized so as to transform the various attribute values into 
comparable values. Therefore, we can obtain the normalized fuzzy decision matrix denoted by 
A = (jj)nxm according to formulas in reference [9]. 

Let 5 = (1 


igo Mag, Ug), then 


min ymin ymin 
lig Mig Wag yee ie 








— porary a. ey . 
Lip max? max? man el Jy, fig = ( ’ )’ 1g € Ja, 
U, Us U; Uig MiG lij 
where uj* = max ui, if 7 € Ji, ion =minl,;, if 7 € Jo, where J; and Jz are the set of benefit 
a t 


criteria and cost criteria, respectively. 
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Step 3. Construct the weighted normalized fuzzy decision matrix. 


W111 W2%1Q 0 ++ Wm®im Yu Y1i2 *** Yim 
A = W1%21  W2%22 +++ WmLam Yor Y22 *** Yam 
W1Ln1 W2En2 ve Wm@lnm Yn Yn2 wre YUnm 
where 
Yig = Wy Vay 
= [a;,8;) - lig, Miz, Vig) 


’ 
= (aajlig, GMs, BV Miz, By Mig) 
= (ai7, bij, Giz, diz), 
i= 1,2,-+: 2,7 =1,2,-°: »™Mm. 


Step 4. Determine the fuzzy positive-ideal solution X* and fuzzy negative-ideal solution 
X~, respectively, and 








+ (at at + = 
x = (U7 19a 17? + Undo = (9, .Yas"** Ym) 

where 
+ + opt ot aty 3 
” — (max a;;, max b;;, max c;;, max d;;) = (a!,b!,c’,d! =1,2,---,m 
#3 (eine Mi<i<n “1<i<n “? 1<i<n ii) = (0587565545 ),5 =1,2,--- ym, 
~ =( min a;;, min b;;, min c;;, min d;;) = (a, ,b,,c,,d,),7 =1,2,--- ,m. 
uy (aca ‘Di<icn “ici<n 9? 1<i<n wi) = NO5 PG AG oy d= Leys 





Step 5. Calculate the distance of each alternative from fuzzy positive-ideal solution X* 
and fuzzy negative-ideal solution X~, respectively. 


The distance of each alternative from X+ and X~ can be calculated as 





1 m 
Dy} = D(X, X*) = =| S [(aig — af)? + (bij — OF)? + (cag — Cf)? + (dig — aT?) , 


~a\% 


Dy = D(Xi,X~) = 5 | Dolley — a5)? + (bij — BF)? + (Cus — oF? + (dey — 45) 


j=l 














i= 1,2,--- ,n. 

Step 6. Calculate the relative closeness coefficient of each alternative to the ideal solution. 

D; 

— eee a 
D; + Dz 


4 


= 1,0)444 ni 


Step 7. Determine the ranking order of all alternatives according to the closeness coefficient 
in the descending order and select the best one. 
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84. Numerical example 


This example is taken from [10], one of the Middle East country intends to purchase a cer- 
tain number jet fighter from US, The American Pentagon’s officials have provided information 
relative to 4 kinds of jet fighter types which are permitted to sell. The Middle East country 
sends out the expert group to carry on the detailed inspection to these 4 kinds of jet fighters, 
Inspection results are shown in Table 1. Which kind of jet fighter should they purchase to make 
total utility value of decision making biggest? 


Table 1: Expert group inspection result to 4 kinds of jet fighters 








Type Maximum Cruising Peak load Price Reliability | Maintainability 
speed (Ma) | radius (mile) (pond) (/$ 10°) 
Xy 2.0 1500 20000 5.5 High Very high 
X2 25 2700 18000 6.5 Low General 
X3 1.8 2000 21000 4.5 High High 
X4 2.2 1800 20000 5:0 General General 























weight w, of attribute G; cannot determine completely (7 = 1,2,--- ,6), but according to the 
expert analysis we know wy, € [0.15,0.20], we € [0.1,0.15], w3 € [0.1,0.15], wa € [0.1,0.15], 
ws € [0.15, 0.20], we € [0.2, 0.25]. 
Step 1. Correspond to normalized fuzzy decision-making matrix as follows: 
(0.8,0.8,0.8)  (0.56,0.56,0.56) (0.95, 0.95, 0.95) (0.82, 0.82,0.82) (0.6,0.8,1) (0.7, 0.9, 1) 
(1,1, 1) (1,1, 1) (0.86, 0.86, 0.86) (0.69, 0.69,0.69)  (0,0.2,0.4) (0.3, 0.5, 0.7) 


(0.72, 0.72, 0.72) (0.74, 0.74, 0.74) (1,1,1) (1, 1,1) (0.6,0.8,1) (0.6, 0.8, 1) 
(0.88, 0.88,0.88) (0.67, 0.67,0.67) (0.95, 0.95,0.95) (0.9, 0.9,0.9) (0.3, 0.5,0.7) (0.3, 0.5, 0.7) 


Step 2. Construct the weighted normalized fuzzy decision matrix. 


(0.12, 0.12, 0.16, 0.16) (0.056, 0.056, 0.084, 0.084) (0.095, 0.095, 0.1425, 0.1425) —< 
(0.15, 0.15, 0.2, 0.2) (0.1, 0.1, 0.15, 0.15) (0.086, 0.086, 0.129, 0.129) => 
(0.108, 0.108, 0.144, 0.144) (0.074, 0.074, 0.111, 0.111) (0.1, 0.1, 0.15, 0.15) => 
(0.132, 0.132,0.176,0.176) (0.067, 0.067, 0.1005,0.1005) (0.095, 0.095, 0.1425, 0.1425) — 


(0.082, 0.082, 0.123, 0.123) (0.09, 0.12,0.16,0.20) (0.14, 0.18, 0.225, 0.25) ] 
(0.069, 0.069, 0.1035, 0.1035) (0, 0.03, 0.04, 0.08) (0.06, 0.1, 0.125, 0.175) 
(0.1, 0.1, 0.15, 0.15) (0.09, 0.12, 0.16, 0.20) (0.12, 0.16, 0.2, 0.25) 
(0.09, 0.09, 0.135, 0.135) (0.045, 0.075,0.1,0.14) (0.06, 0.1, 0.125, 0.175) 


Step 3. Determine the fuzzy positive-ideal solution X* and fuzzy negative-ideal solution 
X~, respectively, 
Xt = [(0.15, 0.15, 0.20, 0.20), (0.1, 0.1, 0.15, 0.15), (0.1, 0.1, 0.15, 0.15), 
(0.1, 0.1, 0.15, 0.15), (0.09, 0.12, 0.16, 0.20), (0.14, 0.18, 0.225, 0.25)] 


X7~ = [(0.108, 0.108, 0.144, 0.144), (0.056, 0.056, 0.084, 0.084), (0.086, 0.086, 0.129, 0.129), 
(0.069, 0.069, 0.1035, 0.1035), (0, 0.03, 0.04, 0.08), (0.06, 0.1, 0.125, 0.175)] 


Step 4. Calculate the distance of each alternative from fuzzy positive-ideal solution and 
fuzzy negative-ideal solution, respectively, as Table 2. 
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Table 2: The distance measurement 





X1 XQ X3 X4 
Dt 0.03522 0.06953 0.03124 0.05555 
D- 0.06885 0.03740 0.06758 0.03412 








Step 5. Calculate the relative closeness coefficient of each alternative to the ideal solution 


C; = 0.66157, C2 = 0.34976, C3 = 0.68387, Cy, = 0.38051. 


Step 6. According to the closeness coefficient in the descending order, the ranking order of 


four alternatives can be determined as follows: 
X3> X,> X4> Xo. 


Therefore, the best selection is X3. 

The example result by method of this paper proposal is consistent with result by method 
of fuzzy compromise decision method which is first synthesis, then weighted decision making 
pattern introduced in reference [10]. Generally speaking, except some decision-making situa- 
tions with quite obvious differences, the result which the different method obtains isn’t always 
consistent. Therefore in the actual decision-making, we should use several different methods to 
solve, synthesize and compare result from each method in order to obtain a relatively reliable 


conclusion. 


§5. Conclusion 


This paper has proven the product of interval number and the triangular fuzzy number 
is a trapezoidal fuzzy number. We have proposed new method for multiple attribute decision 
making problem with interval number attribute weight and triangular fuzzy number attribute 
value. A numerical example is used to illustrate the method. Through the example, we can 
find this method is reliable, the thinking is clear and is one practical decision making method. 
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Abstract For any positive integer n, we define a new Smarandache function G(n) as the 


smallest positive integer m such that II ¢(k) is divisible by n, where ¢(n) is the Euler 
k=1 
function. The main purpose of this paper is using the elementary methods to study the 


elementary properties of G(n), and give three interesting formulas for it. 
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81. Introduction and results 


For any positive integer n, we define a new Smarandache function G(n) as the smallest 


positive integer m such that II ¢(k) is divisible by n. That is, 
k=1 


G(n) = min{m: n | II o(k), me N}, 
k=1 

where ¢(n) is the Euler function, N denotes the set of all positive integers. For example, the 
first few values of G(n) are: G(3) = 7, G(4) = 4, G(5) = 11, G(6) = 7, G(7) = 29, G(8) = 5, 
G(9) = 9, G(10) = 11, G(11) = 23, G(12) = 7, G(13) = 43, G(14) = 29, G(15) = 11, G(16) =5 

--. About the properties of this function, it seems that none had studied it yet, at last we 
have not seen any related papers before. Recently, Professor Zhang Wenpeng asked us to study 
the arithmetical properties of G(n). The main purpose of this paper is using the elementary 
methods to study this problem, and prove the following three conclusions: 

Theorem 1. For any prime p, we have the calculating formulae 


G(p) = min{p”, q(p, 1)}; 
G(p?) = a(p, 2), if q(p, 2) < p?; G(p?) = pv’, if a(p, 1) <p? < q(p, 2); 
G(p?) = a(p, 1), if p? < q(p, 1) < 2p?; and G(p?) = 2p”, if q(p, 1) > 2p?, 


where q(p, 7) is the i-th prime in the arithmetical series {np + 1}. 
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Theorem 2. G(n) is a Smarandache multiplicative function, and moreover, the Dirichlet 


series oe el 


n 





is divergent. 
n=1 
Theorem 3. Let k > 2 be a fixed positive integer, then for any positive integer group 
(m1, M2, +++, Mp), we have the inequality 


G(mim2+++mr) < G(m1)G(m2) --- G(mz). 


Note. For any positive integer n, we found that n = 1, 4, 9 are three positive integer 
solutions of the equation G(n) = n. Whether there exist infinite positive integers n such that 
the equation G(n) = n is an interesting problem. We conjecture that the equation G(n) = n 
has only three positive integer solutions n = 1, 4, 9. This is an unsolved problem. 


§2. Proof of the theorem 


In this section, we shall complete the proof of our theorems directely . First we prove Theo- 
m s 


rem 1. Let G(p) = m, namely p | II o(k), p t II o(k),0<s<m. Ifm= qi'q5?--- qe" be the 
factorization of m into prime ogee then de a ge" (qn — 198?” (qa —1) ++ 1g — 1), 
from the definition of G(n) we can deduce that p divides one of (q;‘), 1 <i <r. So p divides 
q; ‘~} or gq; — 1. We discuss it in the following two cases: 

(i). If p | q;—-1, then we must have m = q; = lp+1 be the smallest prime in the arithmetical 
series {kp +1}. 

(ii). If p | g@*—', then we have have a; = 2, m = q? = p*. 

Combining (i) and (ii) we may immediately deduce that G(p) = min{p?, lp + 1}, where 
lp +1 be the smallest prime in the arithmetical series {kp + 1}. Similarly, we can also deduce 
the calculating formulae for G(p”). This proves Theorem 1. 

Now we prove Theorem 2. For any positive integer n > 1, let n = p{'ps?---pr be the 
factorization of n into prime powers, G(p;’) = mi, 1 <i <r, m= max{m1, m2, --:, mr}, 
then from the definition of G(n) we have pj divides @(1)¢(2)---¢(m;) for alll <i <r. So p}* 
divides $(1)(2)---d(m). Since (p;, pj) = 1,147, so we must have n = pf''p5? --- pe” divides 
(1) 0(2)---o(m). Therefore, G(n) = m = max{G(pt"), G(p$?), +++, G(per)}. So Gin) isa 
Smarandache multiplicative function. From Theorem 1 we may immediately get 

ar > DSP > Lo = 


7 oo oe p ? 


G(n) 


n 








CO 
So the Dirichlet series se is divergent. 


n=1 


Finally, we prove Theorem 3. First we prove that the inequality G(mimz2) < G(m,)G(mg2) 
holds for any positive integer m, and mz. 
Let G(m 1) = u, G(mz2) = v, from the definition of G(n) we can easily get 


my | [[ oli), ma | TI (0. 
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Without loss of generality we can assume wu < v, then 


Uv UU 


Tl 9) =[Lo- T] 9%) = [] 4) - ¢(u+1)--- (x) --- d(2v)--- o(ue) 
k=1 


k=1 k=u+1 k=1 


Notice that 


(2) | P(2v), O(3) | Pv), +++, O(u) | P(ur), P(U+1) | O(U+1), +--+, Oe) | Cr), 


this means that 


[[¢@! I] ¢@. 


i=1 k=utl 
ot UU UV UU 
[[e@][¢@ 1 [Lem 
i=1 i=1 k=1 
Hence 


mym2 | |] ¢(4). 
k=1 
From the definition of G(n) we know that G(m mz) < uv, or 
G(m mz) < G(m1)G(mz2). 
If k > 3, then applying the above conclusion we can deduce that 


G(mimg--- mz) = G(mi(mg-+-mMx)) < G(m1)G(m2-- + mx) 
< G(m1)G(m2)G(m3 +++ mz) 


This completes the proof of Theorem 3. 
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Abstract In this paper, we use the elementary methods to study some divisibility of the 


Smarandache numerical carpet, and obtain two interesting congruent theorems. 
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§1. Introduction and results 


In his book “Only problems, not solutions”, F.Smarandache introduced many interesting 
sequences and arithmetical functions, at the same time, he also proposed 105 unsolved problems 
and conjectures. One of the problems is that he asked us to study the properties of the numerical 


carpet. It has the general form 
1 


lal 
labal 
labcbal 
labcdcbal 
labcdedcbal 
labcdcbal 
labcbal 
labal 
lal 
1 


On the border of level 0, the elements are equal to“1”, they form rhombus. 

Next, on the border of level 1, the elements are equal to “a”, where “a” is the sum of all 
elements of the previous border, the “a”s form a rhombus too inside the previous one. 

Next again, on the border of level 2, the elements are equal to “b”, where “b” is the sum 
of all elements of the previous border, the “b”s form a rhombus too inside the previous one. 


And so on -:: 
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It can be described as 


1 
1 ay ay 1 
1 ay ag... a2 ay 1 
1 a, a2 a3... Gg fin, ay 1 
1 a, ag... a@Q ay 1 
1 ay ay 1 
1 


Let n, k be integers with n > 0 and k > 0, the n — th square can be presented as a set of 
sub-squares. In [2], Mladen and T. Krassimir gave and proved the general formula: 


k-2 
ax = A(k +1) J] 4n+1 - 44), 
i=0 
where integer k > 2. 
Let C(n, k) denotes a, in n-th square, then 


C(n, k) = 1, for k = 0; 
C(n, k) =4(n+1), for k=1; 
k-2 
and C(n, k) = 4(n+1) J] (4n4+1- 41), for2<k<n+l. 
i=0 
The main purpose of this paper is using the elementary methods to study the divisibility 
properties of the sequence {C(n, k)}, and give two interesting congruent theorems. That is, 
we shall prove the following conclusions: 
Theorem 1. Let n be a positive integer. Then for any prime p with 4n—4k+9 < p < 4n+1 
and p= 1 mod 4, we have the congruence 


C(n, k) =0 mod p. 





ut) < 


Theorem 2. Let n be a positive integer, p be any prime with n—k+34 | 3 


1 
pnt ES and p = 3 mod 4, then we have the congruence 


C(n, k) = 0 mod p. 
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§2. Proof of the theorems 


In this section, we shall use the elementary methods to complete the proof of the theorems. 
First we prove Theorem 1. For any positive integer n and prime p= 1 mod 4, if 4n-4k+9 < 
p<4n-+1, then it is clear that 4n —-4k +9 <4n+1 or k > 2, from this we have 


k—2 
C(n, k) =4(n +1) [[ 4n- 4441) 
1=0 
k-2 
=4(n+1) [[(4(n-4) +1). 
1=0 


It is obviously that C(n, 2) = 4(n+1)(4n+1), so if & = 2, then from 4n —4k +9 <p <4n+1 
we know that p= 4n +1, and C(n, 2) =0 mod p. This means that Theorem 1 is true. 
Now we assume that 2 < k <n, so that 








k-2 k-2 
C(n, k) = 4(n +1) [] (4(n i)+1)= (n4 1)-4*. [J (n-i4+4) mod p 
i=0 i=0 
k—-2 
= (n+1)-4*. TI (»-1+ 4) mod p. (1) 
i=0 


; ' — 3pt+i1 
where & denotes the solution of the congruence equation x -% = 1 mod p, and 4 = 7 ; 


4 
Since 4n — 4k +9 < p < 4n+1, so there must exist an integer 0 < i < k — 2 such that 








1 
p=n—it amr . Therefore, from (1) we may immediately deduce that 
k—-2 
C(n, k) =4(n +1) [][ 4(n—2) +1) $0 mod p. 
i=0 


This proves Theorem 1. 
Now we prove Theorem 2. Let n be a positive integer, p be a prime with 


Rae — 


Jsosn+ |S and p= 3 mod 4. 





3 


a 1 
Note that 4 = a then from the definition of C(n, k) we have 








k-2 k—2 
C(n, k) = 4(n+1) [[¢@ )+)= (ntit” II (n—i+4) mod p 
i=0 i=0 
7 | yaa pt 
= n+1)-( i ) I] (»-é+ 24+) moa p (2) 


1 1 
From congruence (2) we know that if n—k+24P7— <p< np Poe, then C(n, k) = 0 mod p. 
p+ 





1 
The inequality n —k+2+ pa <p<nt implies 


ak 1 
n—k+34+"— <psn4 7 
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Since prime p must be a positive integer, so we have 
—k+2 1 
m— kt 34 pa <p<nt ES 


This completes the proof of Theorem 2. 
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Abstract For any positive integer n > 2, the Smarandache Inferior Prime Part {pp(n)} 
is defined as the largest prime less than or equal to n; The Smarandache Superior Prime 
Part {P,(n)} is defined as the smallest prime greater than or equal to n. In this paper, we 
defined some determinants involving the Smarandache prime part sequences, and introduced 


two conjectures proposed by professor Zhang Wenpeng. 
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81. Introduction 


For any positive integer n > 2, the Smarandache Inferior prime part {p,(n)} is defined as 
the largest prime less than or equal to n. The Smarandache Superior prime part {P,(n)} as the 
smallest prime greater than or equal to n. For example, the first few values of these sequences 
are pp(2) = 2, pp(3) = 3, pp(4) = 3, Pp(5) = 5, pp(6) = 5, pp(7) = 7, Dp(8) = 7, pp(9) = 7, 
pp(10) = 7, pp(11) = 11, pp(12) = 11, pp(13) = 18, pp(14) = 18, p,»(15) = 18, p,(16) = 18, 
Pp(17) = 17, pp(18) = 17, pp(19) = 19, ------ : 

P,(1) = 2, P(2) = 2, 3) =3, 2) =5, 2,6) =5, £(6) = 7, P(7) = 7, B(8) = 11, 
P,(9) = 11, P,(10) = 11, P,(11) = 11, P,(12) = 18, P,(13) = 13, P,(14) = 17, P,(15) = 17, 
P,(16) = 17, P,(17) = 17, P,(18) = 19, P,(19) = 19, ------ ‘ 

From the definitions of these sequences we know that for any prime q, we have p,(q) = 
P,(q) = 4. 

In his books “Only problem, Not solutions” and “Sequences of Numbers Involved in 
Unsolved Problem” ?!, Professor F.Smarandache asked us to study the properties of these se- 
quences. About these two problems, some authors had studied them, and obtained some 
interesting results, see references [3], [4], [5] and [6]. For example, Xiaoxia Yan [6] studied the 


the asymptotic properties of rz and proved that for any positive integer n > 1, we have the 
asymptotic formula . 


DH 


*=1+0(n-*), 


| 


where In = {Pp(2) + Pp(3) +--+ + Pp(n)}/n and Sn = {Pp(2) + Pp(3) + +++ + Pp(n)}/n. 


wa 
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The sequences {p,(n)} and {P,(n)} are very interesting and important, because there 
are some close relationships between the Smarandache prime part sequences and the prime 
distribution problem. Now we introduced two determinants formed by the Smarandache prime 


part sequences. They are defined as follows: For any positive integer n, c(n), C(n) aren x n 


determinants, namely 


and 


For example, by definitions and calculating, we can find some values of these determinants 





pp (2) pp (3) pp(n-+ 1) 
ony =| Pett) pp(n +3) pp(2n +1) 
pp(n(n—1) +2) pp(n(n—1) +3) =» pp(n? +1) 
P,(1) P,(2) P,(n) 
omy = | Pent D P,(n +3) P,(2n) 
Pp(n(m—1) +1) Pp(n(m—1) +2) + Pon?) 





as the following table: 













































































n | c(n) C(n) n | c(n) C(n) 

2/1 4 3 | 14 4 

4 |0 188 5 | -96 -1424 

6 | 0 0 7 | 1.0214 x 105 37536 

8 | 0 0 9 |0 0 

11 | 7.8299 x 107 1.1478 x 108 13 | 9.8338 x 108 —2.7958 x 10° 
17 | 8.2462 x 10!4 | 1.3164x10'% =| 19 | —1.9608 x 10! | 1.629 x 1045 
23 | 24545 107 | 3.156 x 10"* 29 | 8.9308 x 10?” —6.5008 x 1077 
31 | 1.1095 x 1029 | —2.2835 x 1078 | 37 | 9.9497 x 1037 —7.4692 x 1038 
Al | —1.1502 x 104? | 3.9244 x 104° | 43 | 1.5152 x 104 —2.1671 x 10“ 
47 | 1.1606 x 10°! | 4.06 x 105° 53 | 2.9359 x 1059 5.8735 x 10°° 
59 | 3.8402 x 10®7 | —3.1043 x 10 | 61 | —3.614 x 107 | 6.9858 x 107? 
67 | 3.2341 x 108° | —9.5374 x 1078 | 71 | —2.219x 10° | 4.7688 x 108° 
73 | —3.6656 x 108 | —3.6985 x 108° | 79 | —4.2038 x 10° | —3.8762 x 1097 
83 | 1.6966 x 1014 | 2.1389 x 10!°4 | 89 | —1.0695 x 1018 | 5.7824 x 10!° 
97 | 2.8.98 x 10!24 | 1.9968 x 10!%4 





About the elementary properties of these determinants, it seems that none had studied 
them, at least we haven’t seen any related papers before. Recently, Professor Zhang Wenpeng 
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asked us to study the properties of these determinants, at the same time, he also proposed 
following two conjectures: 

Conjecture 1. For any composite number n > 6, we have the identities c(n) = 0 and 
C(n) = 0; 

Conjecture 2. For any prime q, we have c(q) 4 0 and C(q) 40. 

From the above table, we believe that these two conjectures are correct. About Conjecture 
1, we have solved it completely, which will be published in Pure and Applied Mathematics. 
But for Conjecture 2, it’s still an open problem. We think that Conjecture 2 is very interesting 
and important, since if it is true, then we can get a new discriminant method to distinguish 
prime from integer numbers through calculating these kinds of determinants. So we introduce 


the conjectures in this place, and hope more scholars who are interested in to study it with us. 
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81. Introduction 


Let x9, %1,%2,--- be integers with x; > 1, 7 > 1. We can express the original fraction Liss 
In 
in the form 2p + a n > 0. Such an expression is called a continued fraction, and 
: be a 
we denote it by (%9,21,°+: ,2n)- 


About the characteristic properties of continued fractions, it can be found in the references 
[1] and [2]. In this paper, we shall study the relationship between the Dirichlet L-functions 
and the continued fractions, and obtain some interesting identities. That is, we shall prove the 
following main conclusion: 

Theorem. Let p bea prime, a be a positive integer with (a,p) = 1 and 1<a<p. Then 
we have the identity 


S- x(a)|L(1, x)? 

x mod p,x(—1)=-1 

n?(p—1) 7 k Qn-1 . a. 3 

= 1) ng (1 1+(-1)")|, 
Dp eae ee) 








where y be the Dirichlet character mod p, L(1, x) be the Dirichlet L-function corresponding 
character x and a/p = pp /Qn = (%0,%1,°++ Zn). 
—1 
For a = 2, 3,4 and - from the above Theorem and the definition of continued fractions 


we may immediately obtain the following four Corollaries. 
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Corollary 1. Let k > 1 be an odd prime, then we have the identity 


2 $7(k) 1 3 
LILA, vl? = % | : | | 
x mod x 12 k 2 k 





Corollary 2. Let p be an odd prime and y be the Dirichlet character mod p. Then 


2 2 
- -1 
: _ "te 0) ifp=1 mod 3; 








_) 36 p? 
S> xB)ILA, P= 38 
— 1)(n — 2)(p — 
x mod p,x(—1)=—1 5 . (P-)@ 5 ip ay ifp=2 mods3. 
Pp 


Corollary 3. Let p be an odd prime and y be the Dirichlet character mod p. Then 


2 —1)?2(p-1 
7 @ (p oy ifp=1 mod4; 


48 
S 4)|L(1,x)/? = Pp 
; x(4)|Z(, x)| an? (p—1)(p?-—6p4+17) . 
x mod p,x(—1)=-1 BR P , ifp=3 mod4. 








Corollary 4. Let p be an odd prime and y be the Dirichlet character mod p. Then 


SE x(4*)waye-- eee 





x mod p,x(—1)=—-1 


It is clear that these four Corollaries are an extension of Walum [3]. 


§2. Some Lemmas 


To complete the proof of the Theorem, we need the following three Lemmas. 


Lemma 1. Let 2,21,--- ,%n be integers and p;, gq; (0 < i < n) be natural numbers, 
continued fraction (a9, %1,°+* ,%i) = Pi Then 
Gi 


PiGi-1 — Pi-19s = (-1)"**",  -1 <i<n. 








and p_2 = 0,p-1 = 1,q-2 = 1,¢-1 = 0, pj = Tipi-1 + Pi-2, G = TiGi-1 + Gi-2,0 Si <n. 
Proof. (See reference[1] ). 
Lemma 2. Let & and h be integers with k > 3 and (h,k) = 1. Then we have 





S(h,k) = =a s x(A)IL(1, x)", 
dk 


x mod d,x(—1)=-1 


. h 
where S(h,q) = S- ((<)) ((“)) is the Dedekind sum and ¢(k) is the Euler’s function. 


qd 
a=1 
Proof. (See reference [4]). 
Lemma 3. Let a and q be two positive integers with (a, q) = 1. Then we have the identity 





Lh a oe ae —1+(-1)" 
S(a, q) = 12 perms +( 1) : g , Te q aE ‘ ) ’ 
k=0 ike 
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where a/q = Pn/dn = (®0,©1,°** ;Ln)- ? 
Proof. Let continued fraction (%,21,--+,2n) = Pn and suppose Pn _ °' For the 
qd q qd 
é Pn Pn-1 : . 
successive terms — and ——., from Lemma 1 we know that 
n n—1 
PnQn-1 — InPn-1 = (-1)"*" nn = 1. (1) 
Using the properties of Dedekind sums and (1) we get 
S(Dn, In) = (Dalai tas In) (2) 


= $(G-1((-1)"*" + pn—19n); In) 
= S((-1)"Giri, an) 
= (-1)""'S (G=1, dn) 
= (-1)"*S(dn1,%); 


where G,—1 denotes the solution x of the congruence equation xgqn—1 = 1( mod q,). From the 





reciprocity formula and Lemma 1 we obtain 


2 2 

Gn + In—1 +1 1 
Ss n—1;4n S no 4n— = 3 
(Qn—1;9n) + S(Qns In-1) Tagen ri (3) 


= 1 (== dn ) +4 1 1 
12 dn dn-1 12¢ndn—1 4 
1 n— n— In n— 1 
12 dn Qn-1 dn Qn-1 4 
Hence, by expression (3) and Lemma 1 we obtain 


Sdn ai, da) (4) 
i Qn-1 | In-2 +1 [Pn — Pn-1 1 
= no 4n— n —1)” 
Sq v) - 12 ( Qn ™ dn-1 “ ) ( ) dn Qn-1 4 


+ 
= ess n n n— 1 
= —S(dn—2;dn—1) + 2 (2 sf, SO on) 1 (-1)r+2 (2 Pp :)) 


dn dn-1 dn dn-1 
1 dn—2 dn-3 Pn-1 Pn-2 1 
= n—3> 4n— + Ln 1)” + 
(9 ve 2) 12 (t= Qn—2 et a ) dn-1 dn—2 A 
sail (2 _ In-2 | tn) ,(-1yttt (2 aol 1 
12 dn Qn-1 . dn Qn-1 4 


: n= eS n n n— 
= S(Gn—3, n—2) 4 1D (2 A es + Ln tra) +(-1) +1 (2-2 2) |. 


dn dn—2 an dn—2 












































If n is an even number, then from (4) we have 











n—-1 
1 dn-1 11) n+l1 (= e) a 
S(dn—1:4n) = S$(q-1;40) + pe eal 1)'tn_4(5 
(n—1; Gn) (4-159) + 75 ( a (=1) oe DA ) ) 





i=0 dn dn qo 


7 3 Sonia. t Gn—1 } ( vad 


dn dn 
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If m is an odd number, then from (4) we have 


S(am-1d) = Slaosm) +35 (m= 2) 4 ay (2B)  F Ds (6) 


Qn 71 dn 71 


= §(1,21)4 * (a= ~) eee (& a **) Ki De 


dn Ly dn LY 




















= (a1 — 1)(a1 — 2) de 1 = 4 pres 2 xo 4 a al 


1221 12] dn In U1 i=0 





1 . i Qn-1 Pn 1 
Ss -1)'tn_; yyrti . 
1 [>t Jatnne+ 2 + (1) | ; 


i=0 my 


Now combining (2), (5) and (6) we obtain 











1 = ~19n-1 Pn as (a1) 
S(Pns Qn) = —1)**ta,+(-1)"" ~ 
(Pas Gn) = = pa yen, + (ayn tet + Be ‘ 
Note that py/dn = a/q we may immediately get 
1 L —~14n-1 a —1+ (aL) 
S(a,q 1g + | : 
(0,0) = = pa Hoy + (ayes 4 8 ‘ 


This completes the proof of Lemma 3. 


§3. Proof of the theorem 


In this section, we shall complete the proof of the theorem. Using Lemma 2 we have 


> x(a)|L(1, x)? = ED) Senay (7) 
x mod p,x(—1)=—1 Pp 


Then from Lemma 3 and expression (7) we can easily get 


» x(a)|L(1, x)? 


x mod p,x(—1)=-1 
_wrp-l1fi/f< a eee n-1dr-1 , @\ | ~1+(-1)" 
Dp a (sx pte ie :) : 8 | 


“ n— 3 
= ae Lp eee a a oy) | 5 
= Gn pp 2 














This completes the proof of Theorem. 
The proof of the Corollaries. 
For any integer k > 1 and (h,k) = 1, using the Mobius inversion formula and Lemma 2 we 


have 





SMEG? = 2? HO sn, (8) 


x mod k,x(—1)=-1 d|k 
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Taking h = 2 in (8), note that $(2,1) = 0 and the continued fraction 


k/d q2 : 2 





,2) = (0,21, £2). 


If k be a odd prime, then from (8) and Theorem we have 





xa oP = 724 > A geo, &) 























x mod k,x(—1)=—-1 d|k 

2 0(k) Md) a. 
=n S(2, 5) 

k Po d d 

1 6(k) u(d) | k41 qa , 2 
= 1 Sf \2 ae 

oP ca ents? 
_ 1°o(k) u(d) (k/d-1 k/d-1 | 2 
~ 12k 2, d ( P ? k/d Ta) 
Pll 3 
“2k 5-al- 


This proves Corollary 1. 
Using Theorem and the definition of continued fraction we can also deduce that 


; —1)*(p-1 
’ @ te ifp=1 mod3; 














_) 36 p 
—1)(p—2)(p—5 
x mod p,x(—1)=-1 le ‘ (p ad 3 \(p ) if p =2 mod3. 
36 p 
ua ae ee, ifp=1 mod 4; 
o.: 48 pe ’ ae io , 
>> x(4)|Z(1, x)| _ 1? (p — 1)(p? — 6p + 17) 


ifp=3 mod 4. 





x mod p,x(—1)=—-1 28 Po) , 


Sx (44 )eayr-- eee 


and 





x mod p,x(—1)=-1 
This completes the proof of Corollary 2, 3 and 4. 
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Abstract <A connectivity judgment method for public transportation network based on set 
idea is studied in this article. Set is successfully applied to solve the problem of connectivity 
judgment on multi-vertex and multi-edge situation. Search level control strategy is adapted 
to control the search level of route, and three buffer tables are designed to improve the search 


speed and reduce the complexity of the algorithm. 


Keywords Traffic management, connectivity, transportation network, search level control. 


81. Introduction 


Public transport is the main travel way for urban civilians and it is also a way advocated 
by many governments and organizations under the environmental protection concept. On one 
hand, along with the rapid construction of urban traffic, public travel becomes more and more 
convenience; on the other hand, problems come forth with multi-route selection. The civilians 
had to balance among fee, time cost, distance and some other areas, especially those in large 
cities or metropolis. And all this must be built on the basis of the pathway connectivity 
judgment, which are the major problem, as well as the common pathway connectivity judgment 
and optimization strategy under the circumstances of huge number of bus stations and bus lines, 


discussed in this article. 


§2. Some known results 


2.1. Problem description 

The pathway connectivity discussed in this article means the connectivity between one 
station and another, directly or through a third or more other stations, then find out the routes 
and the shortest one if they are interlinked. 

Basically, there are three types of public transport network distribution which include, all 


stations on the up-run line and down-run line are the same (also called Bidirectional Lines), 
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part of the stations are the same (also called Unilines) and Loop Lines. Meanwhile there are 
two types of fee system known as Flat Fare and Sectional Fare. Following are three lines’ 
information of a city’s public transport system. 

1135, Bidirectional Line; Sectional Fare. 





$0417 — $0272 — $3425 — $3476 — $2974 — $0234 — $0521 — 53806 — 51682 — 53559 — $0928 — $3564 





80079 — 83175 — $2866 — $3172 — $3269 — $3603 — $1633 — $2578 — $2579 — $2576 — $2577 
ligi, Uniline; Flat Fare. 
Up-Run Line: 


$0238 — $0540 — $0542 — $3024 — $0494 — $0973 — 53077 — $2082 — $2213 — $2210 — $3332 — $3351 








$1414 — $0228 — 50916 — $2811 — $0019 — $0966 
Down-Run Line: 





$0966 — $0019 — $2811 — $0916 — $0228 — $1414 — $3351 — $3332 — $2210 — $2213 — $2082 — $3077 





80973 — $0494 — $3024 — $0542 — $0541 — $0152 — $1765 — $3506 — $0238 
l590, Loop Line; Flat Fare. 





§3295 — $3265 — $1334 — $1333 — $2285 — $1341 $0327 — $0317 — $0325 — $0532 — $1734 — $0314 





$1749 — $1831 — $1337 — $0705 — $0704 — $3295 
Where I, is the n-th line’s No. of the city and s,, is the m-th station of the city. 


2.2. Specialties of the problem 

Algorithm for connectivity of graph (include Breadth First Search and Depth First Search), 
adjacency matrix and genetic algorithm are all popular methods adopted to judge pathway con- 
nectivity [1]. Outwardly, public transport network can be convert to a digraph whose vertexes 
are stations and edges are lines between two stations. However there are actually more than one 
line that go through each station in the transport network, namely, there are unfixed number 
of edges between two adjoining vertexes. Therefore, the commonly used router algorithms for 
connectivity of graph are inapplicable since the time complexity and space complexity increase 
rapidly, for instance, there are more than 850 bus lines and more than 6400 bus stations in city 
Beijing [2], since the transport network is a n-th Iterated Line Digraphs [3]. 

Moreover, on one hand, because of the characteristics of intersecting and overlapping, there 
are always more than one line between two adjoining stations. On the other hand, even if there 
are route between two stations, it may need ongoing transferring, however, it is contrary to the 
public’s travel habits which are the less transferring, the better. According to the experience 
of livelihood, the public trend to reduce their transferring frequency when other factors’ (such 
as cost, time, etc.) difference is not very significant. That is, if there are more than one route 
between two stations, direct line has priority over transferring and transferring once has priority 
over transferring more times and so on, and other factors will be taken into considering only 
when transferring frequencies are equal. Therefore, the following discussing was based on the 
public’s travel habits. 


§3. Set method for connectivity 


3.1. Set description 
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Provide elements 1; , s; and sets L;, S;, Aj, Byj,;), which are 

1;: Element that indicate line No. 7; 

s;: Element that indicate station No. J; 

L,;: Set of all lines that pass through station s;; 

S;: Set of stations that line 1; pass through; 

A; : Set of stations that can be reached start from s,; 

Boy 3 Set of all stations that can be reached with line |; start from station s,; 

The following expressions can be drown out: 

1). Direct Line: the following expression is obtained if there have direct line run from start 
station Sstar¢ to destination station Sges¢: 


Sstart — Sdest <> {Sstart, Sdest } Cc Astart a) Adest- (1) 


Where Sgtart — Sdest indicate there are direct line run from Sgtar¢ tO Sdest- 

Therefore, when {Sstart, Sdest } GC Astart M Adest is establish, namely, there are direct line 
run from Start tO Sdest, the collection C' of all direct lines are the intersection of set Lstart 
whose elements are lines that run through s.4¢,¢ and set Dges: whose elements are lines that 
run through S¢est: 


C= Lstart a Lidest- (2) 


2). Transfer Once: The circumstance of transfer once from S start tO Sdest can be converted 
to the iterate of two direct lines, that is, the direct line from the start station $.¢,¢ to the transfer 
station Stransfer and the direct line from the transfer station Stransfer to the destination station 
Sdest- Therefore the set D of transfer stations is: 


D= Astart a) Adest- (3) 


3). Transfer multi-times: Similarly, The circumstance of transfer n-times from $start to 
Sdest Can be converted to the iterate of (n+ 1) direct lines. 


Direct line and transfer-needed circumstances are schematically shown in Fig.1. 





4} Direct Line b} Transfer Gace epigusneneaed 


Figure 1: The schematic view of (a) direct line, (b) transfer once and (c) transfer twice of 
two stations 


3.2. Optimization of algorithm 
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As has been analyzed that the public always trend to choose routes with less transferring 
frequencies, therefore, direct lines should be searched in top priority; deepen search level when 
direct lines are not available, and so on and so forth, until searching accomplished (succeed 
or failed) [4]. In addition, three memory tables known as line-station hash table, line-station 
searching linked-list and station-line searching linked-list are designed to improve the searching 
speed and reduce the searching complexity. 


3.2.1. Route searching strategy 

Search-level control strategy was introduce to integrate the Breadth-First Search method 
(BFS) and Depth-First Search method (DFS). The concrete steps are list in the following: 

1). Construct set Lstart,; Laest; Astart ANd Adest for Sstart aNd Sdest- 

2). Figure out the intersection of Letarte and Laest, Dstart Laest- Then find out all identical 
elements in both set Detar and set Lges: via BFS, thus the result C is the set of all direct lines. 

a). If C is not empty, record all elements in C' which are all direct lines. 

b). If C is empty, deepen the search levels and continue with step 3); 

3). Figure out the intersection of Agstarz and Adest, Astart™MAdest- Then find out all identical 
elements in both set Agzarz and set Agest via DFS, thus the result set D is the set of all transfer 
stations. 

a). If D is not empty, record all elements in D which are all transfer stations. Then search 
direct lines with sstar_ and D, Saes¢ and D, respectively. 


b). If D is empty, keep on deepening the search levels. 


3.2.2. Hash table structure of line-station 
Hash table structure of line-station stores the relationship of lines and stations. When the 
hash value, which are corresponded a line and a station, is “1” , the corresponding line runs 
through the corresponding station, and “0” , doesn’t. The hash table structure of line-station 
is schematically shown in Fig. 2. 


81 82 S83 +++ Sn-1 Sn 
hy 0 1 O ::: 1 0 
Io 1 0 O ::: 1 0 
Is 0 0 1 ::: 0 1 

Im-1 O 1 O ::: il 1 
bees 0 O 1 ::: 0 1 


Figure 2: Schematic view of the hash table structure of line-station 
As can be clearly seen from Fig. 2 that row direction of the hash table represents the lines 


“1” in row direction marks all 


and column direction, the stations. Moreover, the hash value 
stations that the corresponding line runs through, and in column direction marks all lines that 
run through the corresponding station. As is known to all that the time complexity of search 


algorithm using hash table structure is the constant of O(1) [5]. 
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Two important logic relationships are implicated in the above hast table structure, known 


as: 

1). There are direct line between every two stations whose hash value is “1” in the same 
row. 

2). Transfer is available between every two lines whose hash value is “1” in the same 
column. 


Actually, the purpose of constructing the hash table of line-station is to perform high speed 
search utilizing these two relationships. 


3.2.3. Line-station searching linked-list 
Line-station searching linked-list stores the stations that every line runs through succes- 
sively in a duplex linked-list. The list is sorted by line NO. to locates elements with binary 
search. Line-station search list structure is schematically shown in Fig. 3. 


Figure 3: Schematic view of line-station searching linked-list structure 
Since the time complexity of binary search is O(log(n)) , the total time complexity of the 
line-station duplex searching linked-list is 


O(log(m)) + O(log(max(spl))), (4) 


where m is the count of all lines; max(spl) is the max count of stations that a line runs 
through. 


3.2.4 Station-line searching linked-list 
The station-line searching linked-list is also designed as a duplex linked-list, every of whose 
nodes stores all lines that run through the station. The list is sorted by station NO. to locates 


elements with binary search. Line-station search list structure is schematically shown in Fig. 


4. Goo = 


Figure 4: Schematic view of station- line searching linked-list structure 
Similarly, the total time complexity of the station-line duplex searching linked-list is 
O(log(n)) + O(log(max(Ips))) (5) 
where n is the count of all stations; max(Ips) is the max count of lines that runs through a 
station. 
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3.3. Implementation of algorithm 

The following route searching steps are acquired with integrated consider of search-level 
control strategy and optimization methods: 

With regard to the start station ss, and the destination station S¢dest 

1). Set Dstare and set Luest are obtained from station-line searching linked-list by binary 
searching and select the one with less elements, that is min(Lstart, Lest). 

2). Pick out all elements from min(Lstart, Laest) one by one and search for Sstart O Sdest 
in line-station searching linked-list (If min(Lstart, Laest) = Lstart, search for Start, else search 
for Sdest), record matched lines and store them in set Daten; 

a). If Lmaten # BW, output all elements in Lynatch which are direct lines; 

b). Else if Lmateh = @, deepen the search levels and continue with step 3); 

3). Figure out D = Agtart MN Adest 

a). If D4 @ and D = {di, dg, ...d,}, then d1,do,...d, are transfer stations. Go to step 1), 
perform direct-line searching with 551, and elements in D, as well as elements in D and S¢est 
respectively. 

b). Else if D = @ deepen the search levels and continue with step 4); 


/ 


be / 
4). Pick out Sipare ANd Share 


from Agtare and Agest respectively, then set the searching 
level to 1 and continue with step 1), namely transform the problem to search for direct-line 


/ / 
between Sirar¢ ANd $4.44: 


§4. Summary 


Because of the characteristics of intersecting and overlapping, the commonly used router 
algorithms based on Graph Theory or Genetic Algorithm are inapplicable in the public transport 
network; therefore, a set-based connectivity judging method was proposed in this article. Search 
level control strategy is adapted to resolve the graph problem with huge number of vertexes 
and edges and a hash table and two duplex linked-list are designed to improve the searching 
speed and reduce the searching complexity. 
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Abstract In this paper, some elementary methods are used to study the property of the 


Smarandache-Riemann zeta sequence and obtain a general result. 
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§1. Introduction and result 


For any complex number s, let 
Co 


C(s)= Sok 
k=1 


be the Riemann zeta function. For any positive integer n, let T;, be a positive real number such 
that 


((2n) = 7 (1) 


where 7 is ratio of the circumference of a circle to its diameter. Then the sequence T = {T;, }°2, 
is called the Smarandache-Riemann zeta sequence. About the elementary properties of the 
Smarandache-Riemann zeta sequence, some scholars have studied it, and got some useful results. 
For example, in [2], Murthy believed that T,, is a sequence of integers. Meanwhile, he proposed 
the following conjecture: 

Conjecture. No two terms of 7), are relatively prime. 

In [3], Le Maohua proved some interesting results. That is, if 


ord(2, (2n)!) < 2n — 2, 


where ord(2, (2n)!) denotes the order of prime 2 in (2n)!, then J), is not an integer, and finally 
he defies Murthy’s conjecture. 
In reference [4], Li Jie proved that for any positive integer n > 1, we have the identity 


+00 
ord(2, (2n)!) = ag(2n) = » EB = 2n — a(2n, 2), 


where [2] denotes the greatest integer not exceeding x. 
So if 2n—a(2n,2) < 2n — 2, or a(2n,2) > 3, then T,, is not an integer. 
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In fact, there exist infinite positive integers n such that a(2n,2) > 3, and T,, is not an 
integer. From this, we know that Murthy’s conjecture is not correct, because there exist infinite 
positive integers n such that 7, is not an integer. 

In this paper, we use the elementary methods to study another property of the Smarandache- 
Riemann zeta sequence, and give a general result for it. That is, we shall prove the following 
conclusion: 

Theorem. If 7), are positive integers, then 3 divides T;,,, more generally, if n = 2k, then 
5 divides T,,; If n = 3k, then 7 divides T,,, where k ¥ 0 is an integer. 

So from this Theorem we may immediately get the following 

Corollary. For any positive integers m and n(m # n), if T,, and T;, are integers, then 


(Im; Tn) > 3, (Tam,Tan) > 15, (L3m,T3n) > 21. 


§2. Proof of the theorem 


In this section, we shall complete the proof of our theorem. First we need two simple 
Lemmas which we state as follows: 


Lemma 1. If n is a positive integer, then we have 


nt (277)*" Bon 
2(2n)!? 





¢(2n) = (-1) 


where Bo, is the Bernoulli number. 
Proof. See reference [1]. 


Lemma 2. For any positive integer n, we have 


os (3) 


p—1|2n 


where J, is an integer and the sum is over all primes p such that p — 1 divides 2n. 
Proof. See reference [3]. 


Lemma 3. For any positive integer n, we have 


(2n)!by, 


Tr = 
os 9 
g2n la, 


where a, and b,, are coprime positive integers satisfying 2||b,,3|b,,n > 1. 
Proof. It is a fact that 





Pe i a 
C(2n) = (1) Bay, 1, (5) 
where 
Ban = (1), n>1. (6) 


Using (1), (5) and (6), we get (4). 
Now we use above Lemmas to complete the proof of our theorem. 
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For any positive integer n, from (4) we can directly obtain that if T;, is an integer, then 3 
divides T,,, since (adn, b,) = 1. 
From (1), (2) and (3) we have the following equality 




















talk (27)? Bo (2x)? 1 
i= = ie 2 =(-1)"" md dey 
(en) =a Jt) amy OY Sal Ds, pr]? 
p—1|2n 
Let 
II P= P1P2°""Ps; 
p—1|2n 
where p;(1 <7 < s) is a prime number, and p; < p2--: < ps. 
Then from the above, we have 
—1)rtl . 72n —1)"t1. (2n)! 
— (<1 sn a (an) és 
(27)°" 1 2n—-1 1 
(m-St) ee (n- 5} 
! 
2(2n)! p—ijan? p—lj2n? 
(0 Gal- Tl p 
_ p—1|2n 
es 1 
(tT e- Ie © 3 
p—-1|2n p—1|2n plan” 
- (—1)"t? - (2n)!- pipe +++ ps 
a 1 1 1 
ants (10 pipas-ps = papa Pe: (—+ ~+.2)) 
Pi P2 Ds 
_ (—1)"** - (Qn)! - pipa-++Ds 
2°" «(Tn Papas Pa — Pabla*** Ds — Pipa *** Dy — + —Pipa** P21) 


Then we find that if p;|pypo---ps, 1<%i< s, but 


pit Un + pipe>+*Ps — P2P3°**Ps — P1P3*** Ps —*** — PiP2-**Ds—1)- 


So we can easily deduce that if T;, are integers, when n = 2k, 5 can divide T,,; While 
n = 3k, then 7 can divide Tj. 
This completes the proof of Theorem. 
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Abstract In this paper, a system of generalized set-valued variational inclusions in Banach 
spaces is introduced and studied, which includes many variational inclusions studied by others 
in recent years. By using some new and innovative techniques, several existence theorems 
for the system of generalized set-valued variational inclusions in q-uniformly smooth Banach 
spaces are established, and some perturbed iterative algorithms for solving this kind of system 
of set-valued variational inclusions are suggested and analyzed. Our results improve and 
generalize many known algorithms and results. 
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81. Introduction 


In recent years, variational inequalities have been extended and generalized in different 
directions, using novel and innovative techniques, both for their own sake and for the applica- 
tions. Useful and important generalizations of variational inequalities are set-valued variational 
inclusions, which have been studied by [1-9]. 

Recently, in [1], Chaofeng Shi introduced and studied the following class of set-valued 
variational inclusion problems in a Banach space E. 

For a given m-accretive mapping A, B,C,G : E — CB(E), nonlinear mappings N(-,-), M(-,- 
Ex E — E, single-valued mapping g: H — H, any given f € FE and \ > 0, findu € E,% € 
Au,y € Bu,z € Cu,v € Gu such that 





f € N(%,9) — M(z,v) + Aw(g(u)), (1.1) 


where CB(E) denotes the family of all nonempty closed and bounded subsets of £.Under the 
setting of g—uniformly smooth Banach space, Chaofeng Shi [1] gave the existence and conver- 


gence theorem for the solution of variational inclusion (1.1). 





1 This research is supported by the natural foundation of Shaanxi province of China(Grant. No.: 2006A14)and 
the natural foundation of Shaanxi educational department of China(Grant. No.:07JK421) 
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For a suitable choice of the mappings A, B,C,G,W,N,M,g and f € E, a number of known 
and new variational inequalities, variational inclusions, and related optimization problems can 
be obtained from (1.1). 

On the other hand, X. P. Ding [2] introduced and studied the system of generalized vector 
quasi equilibrium problems in locally G-convex uniform spaces, and the existence results for 
the system of equilibrium problems is obtained. 

Inspired and motivated by the results in Chaofeng Shi [1], S. S. Chang [2], and X. P. 
Ding [3], the purpose of this paper is to introduce and study a class of system of set-valued 
variational inclusions. By using some new techniques, some existence theorems and approximate 
theorems for solving the set-valued variational inclusions in g-uniformly smooth Banach spaces 
are established and suggested. The results presented in this paper generalize, improve and unify 
the corresponding results of Chaofeng Shi [1], S. S. Chang [2], Ding [3,4], Huang [5,6], Zeng [7], 
Kazmi [8], Jung and Morales [9], Agarwal [10], Liu [11,12] and Osilike [13,14]. 


§2. Preliminaries 


Definition 2.1. Let J be a finite or infinite index set. Let {F;};<¢7 be a family of real 

Banach spaces, A, B,C,G : ][~ = [[ CB(Zvbe set-valued mappings, W; : D(W) Cc E; > 
i€l i€L 

2”: be an set-valued mapping, g; : E; — E; be a single-valued mapping,and Nj;(-,-), Mi(-,-) : 

E, x E; — E; be two nonlinear mappings, for any given f; € E; and X > 0, we consider the 


following problem of finding u € ][4.2 € Au,y € Bu,z € Cu,v © Gu such that, for each 
i€l 





(eT 
fi © Ni( Pid, PiG) — Mi( Piz, Piv) + AWi(9:(Piu)), (2.1) 


where P; : [| — E; is a projection operator. This problem is called the generalized set- 
i€l 

valued variational inclusion problem in a Banach space. 

Next we consider some special cases of problem (2.1). 


Qrr={h,ff=f,NMi = N,M; = M,W; = W, 9; = g, then problem (2.1) is equivalent 











to 
f € N(%,9) — M(z,v) + Aw(g(u)), (2.2) 


This problem is introduced and studied by Chaofeng Shi [1]. 


(2) Furthermore, if M = 0,W = A: D(A) — 2” is an m-accretive mapping, A = T, B = F 
and C = G=0, then problem(2.1) is equivalent to finding q € E,w € Tq,v € Fq such that 


f € N(w,v) + AA(g(q)). (2.3) 


This problem was introduced and studied by S. S. Chang [1]. 


For a suitable choice for the mappings A, B,C,G,W;, Ni, Mi, gi, f; and the family of spaces 


EF, we can obtain a lot of known and new variational inequalities, variational inclusions and the 
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related optimization problems. Furthermore, they can make us be able to study mathematics, 
physics and engineering science problems in a general and unified frame [1-9]. 

Definition 2.2. [10] Let & be a real Banach space. The module of smoothness of E is 

defined by 

pxu(t) = sup{ 5 (ile + yll + lle — ll) — [lel] <1, Ilyll < th. 
The space F is said to be uniformly smooth if lim (pe (t) /t) = 0. Moreover, FE is called q- 
uniformly smooth if there exists a constant c > 0 such that px(t) < ct?. 

Remark 2.1. [10] All Hilbert spaces, L,(or l,) spaces (p > 2) and the Sobolev spaces, 
W?P(p > 2) are 2-uniformly smooth, while, for 1 < p < 2, L,(or l,) and W?, spaces are 
p—uniformly smooth. 

Definition 2.3. Let S : EH — CB(E) be a set-valued mapping. S is said to be quasi- 
contractive, if there exists a constant r € (0,1) such that for any p € Sa, q € Sy, 





lp — all Sr max{||a — yl], |e — pl, lle — all, lly — pil ly — all}- 


To prove the main result, we need the following lemmas. 
Lemma 2.1. [15] Let E be a g—uniformly smooth Banach space with g > 1. Then there 
exists a constant c > 0 such that 


[tw + (1 — t)y — 2l|? < [1 — tq — 1] lly — 2]? + tella — 21|¢ 





—t(1 — t7*e)||z — y|4, 


for all x,y,z € E andt € [0,1]. 
Lemma 2.2. [16] Suppose that {en}, {fn}, {gn} and {yn} are nonnegative real sequences 
such that 
enti S (1— fn)en + fagn + Yn, n = 0, 


with {fn} ¢ [0,1], 5° fr = 00, lim gy =0 and 5°, < oo. Then lim e, =0. 
n=0 n=0 
For the remainder of this paper, r and c denote the constants appearing in definition 2.3 
and lemma 2.1. We assume that H(-,-) is the Hausdorff metric on CB(£) defined by 


H(A, B) = max{sup d(z, B), sup d(y, A)}, A, B € CB(E). 
aed yeB 
Lemma 2.3. [17] Let E be a complete metric space, T : E — CB(E) be a set-valued 
mapping. Then for any given € > 0 and any given z,y € E,u € Tax , there exists v € Ty such 
that 
d(u,v) < (1+e)H(Tx, Ty). 


Lemma 2.4. [18] Let X and Y be two Banach spaces, T : X — 2” a lower semi-continuous 
mapping with nonempty closed and convex values. Then T admits a continuous selection, i. e., 
there exist a continuous mapping h: X — Y such that h(x) € Tx, for each x € X. 

Using Lemma 2.3 and Lemma 2.4, we suggest the following algorithms for the generalized 


set-valued variational inclusion (2.1). 
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Algorithm 2.1. For any given x9 € ][ 4.26 € Azo, yp € Bro, 2% € Cx, v6 € Gao , 
ie! 
compute the sequences {x,} and {y,} by the iterative schemes such that, for each 7 € I, 








(2) Pin © = 6a) Pita + Bnl et Peta = Ne Pie Pay, ) + Pe e,) = AWG Pt) 


1 
(3) In € AYn, Pi:Xn = PiEn41 | < (1 + Aa tilPi Aun, PiAyn+1), 





+1 
1 

(4) Gn € Byn, ||Pi9n — PiYn+ill < (1+ Tp HiFi Bun, PiByn+1), (2.4) 
1 

(5) Zn E Cyn, Pen, — Py2n+1| < (1 — )Hi(PiCyn, PiCyn+1), 











(6) Un E Gyn, Pivn — PiUn+1| < (1 T 1 
n+ 


)Hi(PiGyn, PiGyn+41), 





1 
(7) x, € Arn, ||Pix’, — Pyx’,4,|| < (1 + ——)Hi(PAtn, P/Atn+1), 





n+l 
1 
(8) y, € Brn, ||Piy, — Pivnsaill < + ee pHi PiBan, PBen+1), 
1 
(9) z, E Cin, Pe. = Piz,.4| < (1 ba )H;(P;Can, P;C&n+1), 











~ n+1 








1 
(10) wu, € Gz, ||Pio, — Pionill < a+ Ap HiPiGtn, PiGan+1); 


n=0,1,2,:°°. 


The sequence {x,,} defined by (2.4), in the sequel, is called Ishikawa iterative sequence. 


/ 
n 








In algorithm 2.1, if 6, =0, for alln >0, then y, =a, . Take Z, = 2! , tn = Y).,2n = 2 
and v, = v},, for all n > 0 and we obtain the following. 
Algorithm 2.2. For any given ro € ][ 4. 20 € Azo, io € Buy, 2 € Cxo,v0 € Gro, 


i€l 
compute the sequences {2}, {En}, {Gn}, {zn} and {v,} by the iterative schemes such that, for 
each i € J, 


Pitn4i € (1 — An) Pitn + On(fi + Pitn — Ni(PiEn, Pitin) + Mi(Pizn, Pitn) — AWi(gi(Pitn))), 





























En € Atn, || Pin — PiEn4i|| < A+ HP Ato: P;A&n+41), 

Gn € B&n, ||PiGn — Pi¥nsi\| < A+ 1) H(P:Bon, P;B&n41), (2.5) 

Zn € C&n, ||Pizn — Pi2n4i|| < (1+ = (Hi(PiCan, PiCtn41), 

Un © Gan, || Pin — Pivn41|| < (1+ = PHi(PiGtn, PiGan+1), 
n=0,1,2,-:°. 


The sequence {x,,} defined by (2.5), in the sequel, is called Mann iterative sequence. 
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§3. An existence theorem for the system of the generalized 


set-valued variational inclusions 


In this section, we shall establish an existence theorem for solutions of the system of the 
set-valued variational inclusions (2.1). We have the following results. 

Theorem 3.1. Let E be a q-uniformly smooth Banach space, A, B,C,G : ][~ => 

i€l 
[[c8&) be four set-valued mappings,for each i € I, Ni(-,-), Mi(-,-): Ex E> E,g,: E; > 
i€l 
E;,W; : E; — E; be four single-valued mappings. If I—N;(P;A(-), P,B(-))+Mi(PiC(-), PiG(-))- 
AW; (gi (Pi(-))) is quasi contractive, then there exists u € ][4. zE Au,y € Bu,v € Cu, z 
i€l 

Gu, which is a solution of system of the generalized set-valued variational inclusion (2.1). 

Proof. For each i € J, let 





Since S is set-valued quasi contractive, S; is also set-valued quasi contractive. 
It follows from Naddler [17] that S$; has a unique fixed point P;u € E;, that is 


Puwe S (Pu) = f,+Pu- N;(P;A(u), P;B(u)) + M;,(P;C(u), P,G(u)) a AW; (g(P;u)). 





Thus, there exist  € Au,y € Bu, z € Cu,v € Gu such that, for each i € J, 
fi € Ni( Pit, Pig) — Mi(Piz, Piv) + Wi(gi(Piu)), 


that is, u € E is a solution of the generalized set-valued variational inclusion (2.1). 


84. Approximate problem of the system for generalized 


set-valued variational inclusions 


In theorem 3.1, under some conditions, we have proved that there exist u © ][£.2 € 
ier 
Au,y € Bu,z € Cu,v € Gu which is a solution of the system of generalized set-valued varia- 





tional inclusion (2.1). In this section, we shall study the approximate problem for the solution 
of the system of generalized set-valued variational inclusions (2.1). We have the following result. 


Theorem 4.1. Let FE; be a family of g-uniformly smooth Banach spaces, A,B,C,G : 


II Ey, Il CB(£;) be four set-valued continuous mappings, for each i € J, 

i€I i€I 

Ni(-,°), MiC,-) : By x Ej — F;,9; : BE; > E;,w; : BE; — E; be four single-valued continuous 
mappings, and {a,},{8n} be two sequences in [0,1] satisfying the following conditions: 

(i) T— N,(P,A(-), P;B(-)) + Mi(P:C(-), PiGC)) — AWi (gi P;(-)) is quasi contractive, 

(ii) for each ¢ € I, A, B,C,G are M-Lipschitz continuous with the constants j11, f2, W3 and [a, 
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(iii) 0 <A <an,n>0.0< yj <1/2,i=1,2,3,4, 
(iv) Qn(q —1- cr?) <1,er? <q—-1,89"! < 1/e(1 — er), 


col—+ + er?Bn(er? —q+1) <1-er?. 
Then for any given x € || 4.26 € Azo,yo € Bro, 2% € Cxo,v9 € Gao, the sequences 


i€l 
{xn}, {En}, {Gn}, {en} and {v,} defined by algorithm 2.1 strongly converge to the solution 





Uwe II E;,z € Au,y € Bu,z € Cu,v € Gu of the system of generalized set-valued variational 
ie] 
inclusions (2.1) which is given in Theorem 3.1, respectively. 


Proof. For each i € J, in (1) and (2) of (2.4), choose h,, € W;(gi(Pitn)), kn © Wi(gi(Piyn)), 
such that 


Pitnti = (1 — an) Pitn + an(fi + Pin — Ni(Pikn, Pitn) + Mi(Pizn, Pitn) — Akn), 








Let 
C= hha j= Nits Ae) Pa Peta) Ai 





Then 
Piln41 = (1 = Gg) Piling + AnPn, 


Piyn = (1 — Bn)Pitn + Brit n- (4.1) 
Since I — N;(P;A(-), P:B(-)) + Mi(P;C(-), PiG(-)) — AW; (giPi(-)) is quasi contractive, 
lpn — Piul] <r max{||Piyn — Pill, |lpn — Piynll} 
which implies that 
IlPn — Prull? < r9(||Pitn — Piull? + lpn — Pin|l*)- (4.2) 


Note that 
lPn _ rn|| < rd(P;xn, PiYn); n = 0. 


We consider the following cases. 
Case 1: Suppose that d(P;rn, Piyn) = ||Pian — Piyn|| for some n > 0. It follows from (4.1) that 


[Pn — Tall S rll Pian — PiYn|| = 7l]Bn(Pitn — Pn) 


<rBr\|Pitn —Tnll- (4.3) 


Case 2: Suppose that d(P;an, Piyn) = || Pian — 1rn|| for some n > 0. Then we have 
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Case 3: Suppose that d(Pjan, Piyn) = ||Piyn — 1n|| for some n > 0. Using (4.1), we have 
IPn a Tall ps r||Piyn _ Tnll 
=r||(1 — Br)(Pitn — rn)Il 
Case 4: Suppose that d(P;2,, Piyn) = || Pian — py|| for some n > 0. Then we have 
Pn — Pall < r||Pitn — Pall. (4.6) 
Case 5: Suppose that d(P;an, P:yn) = ||Piyn — pn|| for some n > 0. It follows from (4.1) that 
|Pn 7 Tn| < r||Pivn — Pn|| 


rll = EN Ose =P) ae Brltn = Dn) || 
<r(1 = 9pq)||Pte = Pall Pr eallPa= Pall: 


which implies that 





lm — Fall <A) | Py — pall StPitn — Pall (4.7) 
It follows from (4.3)-(4.7) that 
IIPn — Pall? < 74 Pitn — Pall + 74||Pitn — Pall’, 
n>0. (4.8) 
It follows from Lemma 2.1 that 
|Pign — Prull? = | — Bn) Pitn + Ban — Prull? 


= ||(1 — Bn)(Pitn — Pru) + Baltn — Pyu)|l4 
< [1 —Bn(q—-I)]||Pitn — Pyull? + Brellrn — Piull? 





—6,(1 — BI") ||Pitn =f,y||2,n > 0. (4.9) 
Similarly, we have 
—Bn(1 — 82-'e)||Pitn — rn||?,n = 0. (4.10) 


and 
|Pitngi — Pyull? < [1 —-an(q—1)]||Pitn — Piull?+ ene|lpn — Pull? 


—Oy(1 — a2~*e)||Pitn — py||4,n > 0. (4.11) 


By virtue of the condition (iii) and (iv), (4.2) and (4.8)-(4.11), we have 


|Pitn41 — Pyull? < [L—an(q—1)]l| Pian — Prul|? 
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tancr"|||Pign — Prull? + ||Pitn — Pall] 





—an(1 —a%-*e)||Pitn — Pall? 
< [1 — an(q— 1) + erfan(1 — Bala — 1)))||Pitn — Prull? 
+27 r40n Bn|ltn — Pull? + c?rtanBnllpn — Tall? 
+[er4an(1— Bn(q—1)) — a@n(1 — af-*e)]||Pitn — Pall? 


—2cr?anBn(l — eo) ||Pite — Trl? 





< [1 — Ong = 1) Or OL Belg = 1) ) Oro Bn ||| Pian — Pull? 


+2cr%0inGn(er? + (84 te — 1))||Pian — nll? 





+an[e*r?4B, + er4(1 — Br(q—1)) — (1— af *e)]||Pitn — pall? 





< [lL -—an(q—1-er2)(1 + cr26,)]||Pitn — Pyull? 
< [Ll -— an(q -— 1 -cr*)]||Pian — Prull’, 


Set en = ||Pitn — Pyull?, fr = an(q —1— cr"), 
It follows from the conditions (iii), (iv) and Lemma 2.2 that 


lim e, = 0, 
n—co 


that is, 


lim ||Pjv, — Pyul| = 0. 


Note that 
Il?n — Pyul| < rmax{||Pitn — Pyull, |lrn — Pitall}, 


which implies 


T — Pyu(n — ov). 
Thus, from (4.1), we have 


From (2.4) and conditions (ii) and (iii), we have 


_ 7 1 


1 
< (1+ ——)illPitn — Pinal 


which implies that {P,Z,,} is a Cauchy sequence in Fj. 

So, there exists © € Il E; such that P;z,, — P;z. Now we show that Z € Au. 
i€l 

In fact, 
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< ||Piz — PiE_|| + w(14 |Pitn — Pull > 0. 


n+1 ) 
Also, we have 


Yn 2 Y € Bu, zn 3 2 € Cu, vn - v € Gu. 


From (2.4), we have, for each i € J, 
Paee fot Pai NAP a P iy) + Me, Fv) — AW, (a Fu), 


that is, 
fi © Ni(P%, Pig) — Mj (Piz, Piv) + AW; (g;(Pu)). 


The required results. 


Remark 4.1. Theorem 4.1 generalizes Theorem 4.1 in Chaofeng Shi [1], the corresponding 
results in S. S. Chang [2], Agarwal [10], Liu [11, 12], Osilike [13, 14] and others. 


Remark 4.2. Since algorithm 2.2 is a special case of algorithm 2.1, from Theorem 4.1, we 
can obtain the convergence Theorem for algorithm 2.2, the details are omitted. 
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Abstract For generalized reflexive matrices R,S, ie., R’ = R, R? =I, ST = S,S? =T, 
matrix A is said to be (R,S)-symmetric matrix when RAS = A. In this paper, an iterative 
method is established for solving the matrix equation AX B = C with (R,S)-symmetric so- 


lution X. By this method, the solvability of the above matrix equation can be determined 








autonomically, and a solution or the least-norm solution of which can be obtained within 
finite iteration steps. Meanwhile, the optimal approximation solution to a given matrix Xo 


can also be derived by the least-norm solution of a new matrix equation AXB=C , where 





X =X—-—X),C =C—AXoB. F inally, given numerical examples illustrate the efficiency of 
the iterative method. 


Keywords Matrix equation, iterative method, (R,S)-symmetric solution, least-norm solution, 


optimal approximation solution. 


81. Introduction 


Let R;’*” be the set of all m x n real matrices with rank k, SR"*", OR"*” be the 
set of all symmetric matrices, orthogonal matrices in R"*", respectively. Denoted by the 
superscripts T and I, the transpose and identity matrix with order n, respectively. For matrices 
A= (a1, 02,...,@n),B E R™*", a; © R™, &(A) and tr(A) represent its column space and trace, 
respectively; The symbol vec(-) represents the vec operator, i.e., vec(A) = (a}’,a3,...,a2)7; 
A® B stands for the Kronecker product of matrices A and B; Moreover, (A,B) = tr(B7 A) 
is defined as the inner product of the two matrices, which generates the Frobenius norm, i.e., 
Al = VZAAS = Vr. 

Definition 1.1. !!2/For given generalized reflexive matrices R € R?*™, S © RIX", ie., 
RT = R, R? = 1, ST = S, S* =T, we say matrix A € R™*” is (R,S)-symmetric ((R,S)-skew 
symmetric), if RAS = A (RAS = —A). 

The set of all m x n (R,S)-symmetric ((R,S)-skew symmetric) matrices with respect to 
(R, S') is denoted by GSR™*" (GSSR™*"), 

Definition 1.2. Assume M,N € R?*X™, where p,m are arbitrary positive integers, if 
(M,N) = 0, matrices M, N are called to be orthogonal each other. 
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Remark 1.1. 
(a) In this paper, we let R, S be generalized reflexive matrices as in Definition 1.1. 
(b) For VF € GSR™*", Ge GSSR™*", we can verify that (F,G) = 0. 


For above generalized reflexive matrix R, R’ = R implies that the eigenvalues of R belong 
to real field, and the absolute of the eigenvalues equal to 1 since R? = I. Certainly, S has 
similar properties to R. Hence, we have the following assertion. 

Lemma 1.1. For given generalized reflexive matrices R, S as in Definition 1.1, there exist 
unitary matrices U; € OR™*™, and Uz € OR"*"”, such that 

a ns um, g=0/°" * UL. 
0 —-In—r 0 —In-1 

From Definition 1.1 and Lemma 1.1, we generate the following conclusion. 

Lemma 1.2. Suppose X € GSR™*", then 
xX, O 

X=U; Us, 
0 Xe 
where VX; € R™', Xo e ROO), 

From Lemma 1.2, we can easily get a (R,S)-symmetric matrix by choosing different A;, 
but U; (¢ = 1,2) are fixed in Lemma 1.1. 

The problems to be discussed in this paper can be expressed as follows: 

Problem I. For given generalized reflexive matrices R € R™*™, S € R"*", and A € 
RPex™, Be R41, Ce RP*X4, find X € GSR™*", such that 


AXB=C. (1) 


Problem II. When Problem I is consistent, i.e., its solution set Siz is not empty, for given 
matrix Xo € R™*", find X € Sp such that 


X — Xoll = min |X — Xl. 2 
| oll fain | oll (2) 


The linear matrix equation (1) has been discussed widely with unknown Xl, where 
X has particular structures, such as, symmetric, skew-symmetric, reflexive, anti-reflexive, Her- 
mitian etc.. In these literatures, the authors obtained the expressions of associated solutions 
and solvability conditions by general inverse or matrix decompositions: the singular value 
decomposition (SVD), the generalized SVD (GSVD)§! or the canonical correlation decompo- 
sition (CCD)!. In complex field, the inverse eigenvalue problem and procrustes problem on 
(R,S)-symmetric matrices have been studied by William [1,2]. Of course, the (R,S)-symmetric 
solution of matrix equation (1) can also be given by using GSVD or CCD. In addition, by 
establishing iterative algorithm, the symmetric solution, skew-symmetric solution, Hermitian 
least-norm solution of matrix equation (1) have been derived in [10-12], respectively. Moreover, 
the general solution and symmetric solution of matrix equations A,X B, = C,, A,X Bz = C2 
have also been investigated by Sheng [13] and Peng [14], respectively. 

Problem II is the optimal approximation problem, which occurs frequently in experimental 
13] 


design! In equality (2), the matrix Xo may be obtained from experiments, but it is not 
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necessary to be needed form and the minimum residual requirement. The nearness matrix xX 
is the matrix that satisfies the needed form and the minimum residual restriction. 

Motivated by the iterative methods mentioned in [10], in this paper, we will construct an 
iterative method to obtain the (R,S)-symmetric solutions of Problem I and II. This work will 
be accomplished in section 2 and section 3, respectively. In section 4, we will provide some 
numerical examples to illustrate the efficiency of the iterative method. 


§2. The iterative method for Problem I 


In this section, we will propose an iterative algorithm to solve Problem I, some lemmas 
will be given to analyze the properties of this algorithm. For any (R,S)-symmetric matrix, We 
will show that the solution of matrix equation (1) can be obtained within finite iterative steps. 
The iterative method of Problem I is stated as follows: 

Algorithm 2.1. 

Step 1: Input matrices RE R™™, S € R"*", and A € RPX*™, Be RIC € RP*4, 
Choosing arbitrary matrix X; € GSR™*". 

Step 2: Calculate 


Ry =C-—AX,B, 
1 
P= 5 (ATR BT + RAT R,B’S), 


k:=1. 
Step 3: If Ry = 0, stop; otherwise goto next step. 
Step 4: Calculate 








I| Rx |I? 

Xpiy = X Es 

k41 k+ Pe IP k 
Step 5: Calculate 

I| Re ||? 
Regi = C — AXy41 B= Ry — Pele” 
1 ar T T T | Riga ll? 

Prat = 3'A RriiB + RA RrpyiB S) + [Rel ke 


Algorithm 2.1 reveals that X;, P; € GSR™*” (4 = 1,2,...). In the sequel, we prove the 
feasibility of this algorithm. 


Lemma 2.1. The sequences R;, P; (t= 1,2,...) generated by Algorithm 2.1 satisfy that 
|| Rill? 








(Riti, Ri) = (Ri, Ri) — ype re PY 
when j > 1, 
Ral? al? ILR 1 
Riz Rk; — R,, BR; = V2 44 j= : 3 
ee ee aia ee lala o 


Proof. When j > 1, from the Algorithm 2.1, noting that RP;S = P;, we can obtain 
(AP, B, R;) = (P;, ATR; B*) 


_ (p, ABT + RATRBTS\ | (p ART RATRBTS\ 
w) 9 i a) 9 
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2 Tp.pRT TR.BT 
= (PAP, Rj Ps) | (2,4 RB ) (nes, RAT R, 5) 























|| Rj—1 ||? 2 2 
Rll? 
= (PL Pi) — Te a Pe Pia), (4) 
moreover, 
Rill? 
(Rizi, Rj) = ( i TPP (AP,B), BR; 
I| Rill? 
= (Ri, Rj) ~ Tp py (AP.B, Ry) 
= 19 J ||P; ||? ie] ] * 


The proof of (5) implies that (3) holds when j = 1. Furthermore, submitting (4) into (5) 
generates (3) holds. The proof is completed. 
Lemma 2.2. The sequences {R;}, {P;} generated in the iterative process are self-orthogonal, 
respectively. i.e., 
(R;,R;)=0, (P,,P;)=0, i,j =1,2,...,k(R>2), 143. (6) 
Proof. We prove the conclusion by induction. Since (F,G) = (G,F), we only proof (6) 
holds when i > j. 


When & = 2, it is clear from Lemma 2.1 that 
(Re, Ri) = (Ri, Ri) — 


Applying Algorithm 2.1 yields 











ATR,BT + RATR,BTS _ ||Roll? 
P>,P;)= P,P. 
uae ( 2 RAP :) 
| Rall? 
= (ATR2BT, Pi) + Ae (P1, P:) 
I Rol? all 
he APB) Ea 
ae 
IPP I Roll? 
Ried 
= [Re ee a — Be) + OT 
=i (7) 


Assume that (6) holds for k = s, that is, (P;,P;) = 0, (Ps, P;) =0, 7 =1,2,...,s—1. Being 
similar to the proof of (7), we can verify from Lemma 2.1 and the algorithm that (R,.1, Rs) = 0, 


and (P, st+1> Ps) = = 0. 
Now, we can finish the proof if (R.41,R;) = 0, (Ps41, P;) = 0 hold. In fact, when j = 1, 

noting that the assumptions and (4), then 

Rell? 


(Rs 1; Ri) = (Rs, Ri) - 
. IPs? 


(AP,B, R:) 


RIP . 
Spe ees 
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_ Rell? 


P,, P,) = 
[Paz =O 





and connecting with (3), we have 














A? Royi Bt + RAT Ri BTS | ||Rosill? 
P, ae = P;, P 
(Paar P)= ( 2 ia) 
AT Rye B* + RA? Roi BTS ) \|Rs41l? 
= pt par as qe ee 
( 2 I|Fs ||? 
= (AT Roa BTA) 
= (Re41; AP, B) 
I| Pall? 
= Ral eet (Ri — Ra)) 
|| Pill? 
~ ER Ae Ba) =o. = 
Furthermore, when 2 < 7 < s—1, Lemma 2.1 and the assumptions imply that 
|Rs ||? |Rsll?7||R5 ll? 
(Rs41,R;) = (Rs, Rj) — (Ps, P3) + Fae ppp (Ps, Pj-1) = 0. 
ace I|Psll? || Ps ||? || Ryall? : 


Similar to the proof of (8), we can obtain (P;+1,.P;) = 0. Hence, the conclusions hold. 
Remark 2.1. From Lemma 2.2, we know that, R; (i = 1,2,...,pq) can be regarded as 
an orthogonal basis of matrix space R?*? for their orthogonality each other, so a solution of 
Problem I can be obtained by the Algorithm 2.1 at most pq +1 iteration steps. 
Lemma 2.3. When matrix equation (1) is consistent, provided that X is an arbitrary 
solution, then R,, P, generated by Algorithm 2.1 satisfy that 


(X= Xia) Rel? RS 1, 22x: 
Proof. When k=1, from the Algorithm 2.1 and Lemma 2.2, we get 
A? R,BT + ans 





2 


- (x, AES BA BTS (x ABTS BA PTS 
2 ” 2 


= (X, ATR, BT) — (X;, ATR, BT) 


(X —X4,P,) = (x Xi, 





= (A(X — X1)B, Ri) 
= ||Ri||?. 
Assume that the conclusion holds for k = s, ie., (X — X,,P,) =|| Rg ||?, then 











ri de IRs II? 
(X =X pai Ps) (X Xs |Pal20 2 
_ ux — WRsl/? 
—_ (Xx Age) \| Ps ||? (Pais) 
= 1,2 - el Rae 


|| Ps]? 
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= 0, 
and AT L AT 2 2 
v7 v7 Ss B Ss B Ss 
(% — Xeut, Pays) = (x Xeat, R,41B* + RA* Roi BS eo |B. + P.) 
2 I|Rs| 
lis ir T r 7 |Rs+ill? , ¢ 
= 5X — Xs41,4 Rs4iiB + RA RsiiB 5) + “Trae ~ Kota Po) 


= (X — X541, AT Rs41B") 
= (A(X — Xo41)B, Re+1) 
= ||Resill?. 


Therefore, we complete the proof by the principle of induction. 

Remark 2.2. Lemma 2.3 implies that, if the linear matrix equation (1) is consistent, then, 
R; = 0 if and only if P; = 0. In other words, if there exist a positive number ¢ such that R; 4 0 
but P; = 0, then matrix equation (1) is inconsistent. Therefore, the solvability of Problem I 
can be determined automatically by Algorithm 2.1 in the absence of roundoff errors. 


Based on the previous analysis, we have the following conclusion whose proof is omitted. 

Theorem 2.1. When matrix equation (1) is consistent, for arbitrary initial matrix X, € 
GSR™*", a solution of which can be obtained within finite iterative steps in the absence of 
roundoff errors. 

To facilitate the statement of our main results, we cite the following conclusion. 

Lemma 2.4. (see [10]) Suppose that the consistent linear equations My = b has a solution 
yo € R(M7*), then yo is the least-norm solution of which. 

Theorem 2.2. When matrix equation (1) is consistent. If initial iteration matrix X, = 
AT HB?’ + RATHB’S, where arbitrary H € R?*%, or especially, let X; = 0 € R?*4, then the 
solution generated by the iteration method is the unique least-norm solution of Problem I. 

Proof. Algorithm 2.1 and Theorem 2.1 imply that, if let X, = A7HB? + RATHB’S, 
where H is an arbitrary matrix in R?*4, a solution X* of Problem I is of the form X* = 
ATY BT + RATY B'S. So it is enough to show that X* is the least-norm solution of matrix 
equation (1). 


Considering the following linear matrix equations with X « GSR™*” 


AXB=C 


(9) 
ARXSB=C 


Obviously, the solvability of (9) is equivalent to that of matrix equation (1). Denote 
vec(X*) = a*, vec(X) = x, vec(Y) = y, vec(C) = c, then the above matrix equations can be 
transformed into 

BT@A c 
(BTS) @ (AR) c 


In addition, the solution X* can be rewritten as 
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TP T 
BT@A y BT@A 
C= ER ; 
(BT S) @ (AR) y (BTS) @ (AR) 


which implies, from Lemma 2.4, that X* is the least-norm solution of the linear matrix equations 


(9). We complete the proof. 


§3. The solution of Problem II 


Suppose that matrix equation (1) is consistent, i.e., Sz is not empty. It is easy to verify 
that Sg is a closed convex set in matrix space GS'R™*” under Fronbenius norm, so the optimal 
approximation solution is unique. Without loss of generality, we can assume the given matrix 
Xo € GSR™*” in Problem II because of the orthogonality between (R,S)-symmetric and (R,S)- 


skew symmetric matrix. In fact, for X € Sz, we have 
2 

















xX x xX — BY 
|X — Xo? = x - 240 _ Ane 
2 2 
Xt RS? . || Xo— RSI? 
=|—° 9 | ee 














Let X = X — Xo; C=C- AXoB, then Problem II is equivalent to find the least-norm 
solution X* € GSR™*” of the following matrix equation 


AXB=C. (10) 


From Theorem 2.2, if let initial iteration matrix X, = A?HBT + PATHB™P, where 
arbitrary H € R?*4, or especially, let X; = 0 € R™*”, we can obtain the unique least-norm 
solution X* of matrix equation (10) by applying Algorithm 2.1. Furthermore, the unique 
optimal approximation solution X to Xo can be obtained by X = X* 4+ Xp. 


84. Numerical examples 


Example 1. Input matrices A, B, C, and generalized reflexive matrices R, S' as follows: 


5 -3 O 3S O 2 8B if &@ 2 7 @.o 

GO =4 -6 4 =6 0 -1l4 oo <1. o. © 0 

—§ 0 7 @ FT 8B 1 0-1 0 00 0 0 
A= 0 5 -8 -5 -3 0 8 ,R=]0 0 0 00 0 1 f, 

4 -7 @ TF O -8 -3 Oe h Mt we a 

—1 0 +6 0 —5 Go @ 0 0 0 00 -1 0 

GO =-8 0 8 =F O Ff 0 0 0 10 0 0 
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-3 5 -5 -2 5 -12 0 1 0 0 0 0 

0 4 9 9 4 —6 10 0 0 0 0 

6 -l 7 0 -l1 3 0 0 -1 0 O 0 
B= ; S= ’ 

—2 4 0 5 4 i) 0 0 0 1 0 0 

-l -6 -2 O -6 2 0 0 0 0 0 -!1 

0 -9 1 1 -9 2 0 0 0 0 -1 O 


254.7246  —420.4859 —4.5129 —444.0827 —420.4859 95.2199 
—32.1326 -—896.3380 -—642.6184 —414.3636 —896.3380 78.8874 
285.7966 —281.1569 169.3670 —306.3903 —281.1569 472.0991 
C= —219.4309 907.1669 255.3442 601.2957 907.1669 —332.3348 
—63.6234 —937.6031 —462.9456 —449.7040  —937.6031 149.2335 
352.1195 —511.6855 —185.5221 —550.9943 —511.6855 224.8234 
300.0260 —987.7629 —247.5968 —678.6078 —987.7629 862.4447 
Next, we will find the associated solutions of Problems I and II by Algorithm 2.1 and 
Matlab software. For the influence of roundoff errors, R; will usually unequal to zero in the 
Fronbenius norm. Therefore, for any chosen positive number ¢, however small enough, e.g., 
€ = 1.0e — 010, whenever ||R;|| < ¢, stop the iteration, and X;, is regarded as a required 
solution. 
(I) The solutions of Problem I: 


Choose an arbitrary initial matrix 


12.0207 12.0207 0 14.0000 6.3639 —6.3639 
—6.9999 21.9996 —14.1420 10.6065 21.9996 1.0000 
—21.9996 6.9999  —14.1420 —10.6065 1.0000 21.9996 


X= —34.9993 —62.9988 44.5473 —19.0917 —15.9997 3.9999 ) 
9.8994 9.8994 0 16.0000 12.7278  —12.7278 
18.3846 —18.3846 —75.0000 0 —24.0414 —24.0414 


—62.9988 —34.9993 —44.5473 —19.0917 —3.9999 15.9997 
by Algorithm 2.1 and 36 iteration steps, we obtain a solution of Problem I, 


6.3639 6.3639 0 4.0000 2.1213 —2.1213 
0.0221 7.9080  —4.0470 2.7378 0.8945 —8.1973 
—7.9080 —0.0221 —4.0470 —2.7378 —8.1973 0.8945 


X36 = —0.9779 —6.0917 2.3169 —5.7474 1.8945 3.8025 ) 
4.2426 4.2426 0 —0.0000 2.8284 —2.8284 
3.5355  —3.5355 8.0000 0 1.4142 1.4142 


—6.0917 —0.9779 —2.3169 —5.7474 —3.8025 —1.8945 
with ||Rgol| = 3.8765 — 011 <e, and ||_X36|| = 27.3649. 


If let initial matrix X,; = A7HB? + PA’HB’P, where arbitrary H € R’*®, then the 
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solution X; generated by Algorithm 2.1 is the unique least-norm solution of matrix equation 


(1) by Theorem 2.2. Choosing arbitrary matrix 


-4 1 -2 O 1 0 
6 0 -1 1 0 -1 
—2 1 0 -1 1 0 
A= 1 0 0 2 2 1 ; 
0 -1 0 -1 O 0 
2 0 1 0 -1 1 
0 -1l O 1 1 0 
then the least-norm solution X* of Problem I is 
6.3639 6.3639 0 4.0000 2.1213 
—0.3925 8.3779  —4.2023 3.3713  —0.6099 
—8.3779 0.38925 —4.2023 —3.3713 —6.4524 
X* = X37 = —1.38925 —5.6218 2.1616  —5.1139 0.3901 
4.2426 4.2426 0 0.0000 2.8284 
3.53855  —3.5355 8.0000 0 1.4142 
—5.6218 —1.3925 —2.1616 —5.1139 —5.5474 
with || R37|| = 8.4542e — 011 < e, and || X*|| = 26.9134. 
(II) The solution of Problem II 
Suppose that the given matrix Xp € GSR™® is 
—1.4142 —1.4142 0 —1.0000 —2.8284 
2.5000 2.5000 —3.5355 —4.2426 —2.5000 
—2.5000 —2.5000 —3.5355 4.2426  —6.4999 
Xo = —0.5000 —7.4999 6.36389 —3.5355  —1.0000 
0 0 0 9.0000 —1.4142 
4.2426  —4.2426 7.0000 0 2.1213 
—7.4999 —0.5000 —6.3639 —3.5355 —6.9999 


—2.1213 
—6.4524 
—0.6099 
5.5474 ; 
—2.8284 
1.4142 
—0.3901 


2.8284 
—6.4999 
—2.5000 

6.9999 

1.4142 

2.1213 

1.0000 


Compute Cp = AXoB, let X = X — Xo, C = C — Cp, then we can obtain the optimal 
approximation solution X of Problem II by finding the least-norm solution X* of equation (10), 


that is, 
7.7781 7.7781 0 5.0000 4.9497 
—2.8110 5.7857 —0.6363 7.4894 2.1856 
—5.7857 2.8110 —0.6363 —7.4894 —0.2953 
X*= —0.8110 1.7858  —4.1718 —1.7029 1.6856 
4.2426 4.2426 0 —9.0000 4.2426 
—0.7071 0.7071 1.0000 0 —0.7071 
1.7858 —O0.8110 4.1718  —1.7029 = 1.7953 


Hence the solution of Problem II is 


—4.9497 
—0.2953 
2.1856 
—1.7953 |, 
—4.2426 
—0.7071 
—1.6856 
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6.3639 6.3639 0 4.0000 2.1213  —2.1213 
—0.3110 8.2857 —4.1718 3.2468 —0.3144 —6.7952 
—8.2857 0.3110 —4.1718 —3.2468 —6.7952 —0.3144 


X= —1.3110 —5.7141 2.1921 —5.2384 0.6856 5.2046 
4.2426 4.2426 0 0 2.8284 —2.8284 
3.5355  —3.5355 8.0000 0 1.4142 1.4142 


—5.7141 —1.3110 —2.1921 —5.2384 —5.2046 —0.6856 
Example 2. Given matrices A, B,C being the same as example 2 in [8], 





2 -10 1 
1 -1 0 38 20 3 -22 2 
-3 0 1 -1 
A= 1 -3 -4 4 1, B= ,C=]| 24 24 -72 6 . 
0 —-2 4 #1 
3 1 4 2 16 -18 28 -2 
1 -2 1 1 
and 
1 0 0 0 01 0 0 
0 0 0 -1 10 0 0 
R= > = 
0 0 1 0 00 -1 0 
0 -1 0 O 00 0 1 


Choosing positive number ¢ = 1.0e — 005, we will stop the iteration when ||R,|| < ¢, or 


|| Rx|| > €, but ||Px|| < ¢. Let initial iterative matrix X, = 0 € R**4, by the iterative method, 
we obtain that || Rg|| = 2.1149e + 003 > «, and ||P.|| = 7.4838e — 009 < «. Then, Theorem 2.1 


implies that Problem I has no solution in GSR**+. 
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Abstract For any positive integer n, let S(n) denotes the Smarandache function. The main 
purpose of this paper is using the elementary method to study the solvability of the equation 
S(a)-? + $(b) * = S(c)~’, and obtain its all positive integer solutions. 
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§1. Introduction and result 


For any positive integer n, the famous Smarandache function $(n) is defined as the smallest 
positive integer m such that m! is divisible by n. That is, 


S(n) =min{fm: n|m!, me N}, 


where N denotes the set of all positive integers. For example, the first few values of S(n) are: 
S(1) = 1, S(2) = 2, $(3) = 38, S(4) = 4, S(5) = 5, S(6) = 3, S(7) = 7, 5(8) = 4, S(9) = 6, 
S$(10) = 15, S$(11) = 11, $(12) = 4, $(13) = 18, $(14) = 7, $(15) = 5, S(16) = 6, S(17) = 17, 
S(18) = 6, S(19) = 19, S(20) = 5, ---. From the definition of S(n) we can easily get the 
following conclusions: if n = p%, S(n) = ap, a < p; S(n) < ap, a> p. Let n = ptt ps? --- per 
denotes the factorization of n into prime powers, then S(n) = max {S (ps) }. 

About other properties of S(n), many scholars had studied it, and obtained some interesting 
results, see [1], [2], [3], [4] and [5]. For example, in reference [5], Xu Zhefeng proved that for 


any real number « > 1, 





3\ 2 3 
(s(n) — Pin)? = "SE" 40 7 
pe 3lnz In* x 
where P(n) is the greatest prime divisor of n, ¢(s) is the Riemann zeta-function. 

In reference [6], Rongi Chen and Maohua Le proved that the equation $?(n) + S(n) = kn 
has infinitely many positive integer solutions for every positive integer k. 

On the other hand, Bencze [7] proposed the following problem: Find all positive integer 
solutions of the equation 
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About this problem, it seems that none had studied it yet, at least we have not seen any 
related papers before. In this paper, we use the elementary method to study this problem, and 
solved it completely. That is, we shall prove the following: 

Theorem. The equation (1) has infinite positive integer solutions. And each (a, 6, c) 
satisfy the equation (1) if and only if S?(a) + $?(b) = 2z*, which means S(a), S(b), z are 
Pythagorean triples, and GCD(S(a), $(b)) = d>1,z| a2, S(o) = ~”, where GCD(S(a), S(b)) 


|2| 
denotes the greatest common divisor of S(a) and S(b). 
From this Theorem we can obtain many positive integer solutions of the equation (1). For 
example, if we take 5(@) = 3; SO) =4, 7 =5, d= 5, S(¢) = 12, then we-can get a= 5*, b= 5"; 


e=3°, If d= 10, then a= 5", b= 5", c= 3". 








§2. Proof of the theorem 


In this section, we shall complete the proof of our theorem. First we need the following: 

Lemma. For every positive integer m, the equation S(n) = m has positive integer solution. 

Proof. From the definition of S(n) and reference [8] we can easily get this conclusion. 

Now we use this Lemma to prove our theorem. For convenience, we denote S(a) by x, S(b) 
by y, S(c) by w. Then the equation (1) become 


1 1 1 

et po @) 
It is clear that if one of S(a), s(b) and S(c) is 1, then (a, b,) does not satisfy the equation 

(1). So without loss of generality we can assume that x > 2, y > 2, w > 2. Obviously the 





solvability of the equation (1) is equivalent to the solvability of the equation (2). Now we 
consider the positive integer solutions of the equation (2). From equation (2) we can obtain 
wy? = (0? + yyw (3) 
and 
vy? er 
= = 
ao ye 


(4) 


Suppose that the equation (1) has positive integer solutions, then there exist z, y, w 
satisfied the equation (3). Noting that « > 2, y > 2, w > 2, form the equation (3) we can get 
that 2? + y? is a complete square, in other words, there exists an positive integer z such that 
x? +y? = 27, and soz, y, z are Pythagorean triples. Hence we can get 


(7B) =u. (5) 


z 


Suppose GC'D(xz, y) = d, we can get GCD(a, y) = GCD(a, z) = GCD(y, z) = din 
equation (5). Suppose z = da’, y = dy’, z = dz’, then we have 


dax'y! 
(146) <0 : 
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from equation (6), w is integer, then z’ | dx’y’. Note that GCD(2’, 2’) = GCD(y’, z') =1 and 
GCD(a'y’, z') = 1, then we can get z’ | d and d> z' > 2. Because z = dz’, z | d?, we can get 
w=” trom equation (5). 

At last, we get the positive integer solutions of equation (2), and the solutions satisfying 
x? +y? is a complete square, in other words, there exists a positive integer z such that 2?+y? = 
z7, and GCD(a, y) =d > 2,2z|d?, w= “Y From Lemma we can get the positive integer 
solutions of equation (2). Since there are infinite Pythagorean triples x, y, z satisfied the 
above situation, so the equation (2) has infinite positive integer solutions. And the solvability 
of equation (1) is equivalent to the solvability equation (2), so the equation (1) also has infinite 
positive integer solutions. This completes the proof of Theorem. 
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§1. Introduction and result 


Divisibility theory is the fundamental part in mathematics. The study of the divisibility 
theory is very important and interesting contents in number theory, many authors pay much at- 
tention to this topic, and some deeply results have been obtained. For example, F.Smarandache 
[1] proved that for any positive integer m, we have 


m | (a —a)-(m—I)!. 
Maohua Le [3] improved the above result, and obtained the following conclusion: 


“ m 
m | (a™—a)- FE 
Because (m, m— 1) = 1, it is obvious that m | (a™ — a) -(m-— 2)!. In this paper, we using the 
elementary method to study this problem, and give another form of the Smarandache divisibility 
theorem. That is, we shall prove the following conclusion: 
Theorem. If a and m are two integers with m > 1, then we have 


(a™ — a)-(m-— 2)! 
m—-1 
| 4 


where [2] denotes the largest integer less than or equal to x. 


m | 





2. 





§2. Some preliminary lemmas 


In this section, we shall give several simple Lemmas which are necessary in the proof of 
our theorem. They are stated as follows: 

Lemma 1. If a and m are two integers with m > 0, then (a” — a)(m— 2)! can be divided 
by m. 
Proof. (See reference [2]). 
Lemma 2. Let a and 0 are two positive integers, S = {1, 2, 3, ..... a}. Then the number 


of all positive integers in S which can be divided by b is [F] ; 
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Proof. (See reference [1]). 
Lemma 3. Little Fermat’s Theorem: For any integer a and any prime p, we have the 


congruence 


a? =a mod p. 


Proof. (See reference [3]). 


§3. Proof of the theorems 


In this section, we shall complete the proof of our theorem directly. It is clear that we just 


need to prove that 





| | Stee ae (1) 


Now we separate m into three cases: 
(a). Ifm=1, 2, 3, 4, then from Lemma 1 we know that (1) is true. 
m—-1 


(b). Ifm =5, then a =e laa =" 


(a™—a)-(m—2)! (a®—a)3! (a—l)a(a+1)\(a?+ 13! 


mB 5 





We also know that there exists at least an even number among a, (a+ 1) and (a— 1). So 
we must have 4 | (a? — a) - (5 —2)!, note that (4, 5) = 1, we may immediately deduce that (1) 


is correct. 
m—1 





1 


If m= 8, then ee = | 4 |u4 


(a™—a)-(m—2)! (a®—a)-6!  a-(a”—1)-6! 





m 8 7 8 


It is clear that there exists an even number between a and (a’ — 1), so we have 2 | (a® — a), 
4? | 6!. That is to say, our Theorem is true if m = 8. 

(c). If m > 5, and m # 8, then we can discuss it in two cases: 

If m is an odd number, we can write m = 4k 4+ t, where t= 1 or 3, k © N. Then 


m—1 [4k +t—1] 
lees ~al 4 J Jak 


If m = 4k +1, then we have the estimate 


ee] = Pete] 


= (2k—-1)+(k—-1)+-+-> (3k—2) > 2k. (2) 
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If m = 4k + 3, then we also have 


> ar eat araal ee la 


= W+k+--->3k> 2k. (3) 








From (2) and (3) we may immediately deduce that 4” is a divisor of (4k — 2)!. Let (m — 2)! = 
4* . D, where D is a positive integer. Note that (4k +t, 4”) = 1, so (4k +t) | (a™—a)-D. 
Therefore, if m > 5 is an odd number, then (1) is correct. 

If m is an even number, we can write m = 2"-q, where qg is an odd number. It is clear that 
if r = 1, then from the above we know that the (1) is true. So we can assume that r > 2. This 


time, in order to prove our theorem, we just need to prove that 





92") .q—24r | (a?"4 Zz a)(2"q _ 2)! ; 


Similarly, we also have 


[m2 








| 7 al | eS pen | bare 











= Bool 1 Og 1 td pa ewe 

S Bly pe ha egal 

= g(2” +2"? 4---4+1)-r 

= q2"-—q-r. (4) 


Note that r > 2 and (2, g) = 1, so we have (g-2"—q—r)—(2—)).q—24+1r) = 2"-!q—q4+2—2r > 0. 

Thus, 22” '?-2+" is an divisor of (m—2)!. Since q < m—2, we also have q | (m—2)!. Therefore, 

from (2), (3) and (4) we know that (1) is correct, ifm > 5 (m #8 ) is an even number. 
Combining (a), (b) and (c) we may immediately complete the proof of our Theorem. 
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Abstract The main purpose of this paper is using the elementary method to study the 


properties of the permutation sequence, and prove some interesting conclusions. 
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§1. Introduction and result 


For any positive integer n, the F.Smarandache permutation sequence {P(n)} is defined 
as P(n) = 135---(2n — 1)(2n)---42. For example, the first few value of the sequence {P(n)} 
are: P(1) = 12, P(2) = 1342, P(3) = 135642, P(4) = 13578642, ------ . This sequence 
was introduced by professor F.Smarandache in reference [1], where he asked us to study its 
elementary properties. About this problem, many people had studied it, and obtained a series 
valuable results, see references [2], [3] and [4]. In reference [5], F.Smarandache proposed the 
following problem: Is there any perfect power among the permutation sequences? That is, 
whether there exist positive integers n, m and k with k > 2 such that P(n) = m*. 

The main purpose of the paper is using the elementary method to study this problem, 
and solved it completely. At the same time, we also obtained some other properties of the 
permutation sequence {P(n)}. That is, we shall prove the following: 

Theorem 1. There is no any perfect power among the permutation sequence. 

Theorem 2. Among the permutation sequence {P(n)}, there do not exist the number 
which have the factor 2*, where k > 2; There exist infinite positive integers a € {P(n)} such 
that 3? divide a. 


§2. Proof of the theorems 


In this section, we shall use the elementary method to prove our Theorems directly. First 
we prove Theorem 1. For any positive integer n > 2, note that P(n) be an even number and 


P(n) = 135---(2n — 1)(2n) ---42 
= 10% +43 x 10% +--+ (2n—1) x 10% + (2n) x 10% 4--- 441042. (1) 


So from (1) we know that 


P(n) = 10 +3 x 10% +---+(2n—1) x 10% 4 (2n) x 10+ +---4+4x 10+2=0 mod 2 
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and 


P(n) = 10% 43x 10° 4+---+(2n—1) x 10% + (2n) x 10°41 +---+4x 1042 
= 2mod 4. (2) 


That is to say, 2 | P(n) and 4} P(n). So P(n) is not a perfect power. Otherwise, we can write 
P(n) = m*, where k and m > 2. Since P(n) be an even number, so m must be an even number. 
Therefore, P(n) = m* = 0 mod 4. Contradiction with 4 + P(n). This proves Theorem 1. 
Now we prove Theorem 2. For any positive integer n > 2, from (2) we may immediately 
deduce that 4" | P(n), where r > 2 is an integer. 
Now we write P(n) as: 


P(n) = 10% +3 x 10% +--+» + (2n—1) x 10% + (2n) x L0%*+1 4---+4x% 1042. (3) 


Note that 10 = 1 mod 9, from the properties of the congruence we may get 10° = 1 mod 9 for 
all integer t > 1. Therefore, 
P(n) = 10% 43x 10° 4+---+(2n—1) x 10% + (2n) x 10°"+1 +---+4x 1042 
= 142434---4+---4+2n=n-(2n+1) mod 9. (4) 





From (4) we know that P(n) = 0 mod 9, if n = 0 or 4 mod 9. That is to say, for all integers 
n> 2, P(9n) and P(9n + 4) can be divided by 9 = 3?. This completes the proof of Theorems. 
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Abstract For any positive integer n > 1, the Smarandache quotients Q(n) is defined as the 
smallest positive integer k such that n-k is a factorial number. That is, Q(n) = min{k : 
n-k = ml}. The main purpose of this paper is using the elementary method to study 
the properties of the Smarandache quotients sequence, and give an identity involving the 


Smarandache quotients sequence. 
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81. Introduction and results 


In his book “Only problems, not solutions”, professor F.Smarandache introduced many 
functions, sequences and unsolved problems, many authors had studied it, see references [1], 
[2] and [3]. One of the unsolved problems is the Smarandache quotients sequence {Q(n)}, it 
is defined as the smallest positive integer & such that n-k is a factorial number. That is, 
Q(n) = min{k: n-k = m!}, where m is a positive integer. For example, from the definition of 
Q(n) we can find that the first few values of Q(n) are Q(1) = 1, Q(2) = 1, Q(3) = 2, Q(4) = 6, 
Q(5) = 24, Q(6) = 1, Q(7) = 720, Q(8) = 3, Q(9) = 80, Q(10) = 12, Q(11) = 3628800, 
Q(12) = 2, Q(13) = 479001600, Q(14) = 360, Q(15) = 8, Q(16) = 45, ------ ‘ 

In reference [4], professor F.Smarandache asked us to study the properties of the sequence 
{Q(n)}. About this problem, some authors had studied it, and obtained several simple results. 
For example, Kenichiro Kashihara [5] proved that for any prime p, we have Q(p) = (p— 1)!. 

In this paper, we use the elementary method to study the properties of an infinite series 
involving the Smarandache quotients sequence, and give an interesting identity. That is, we 
shall prove the following conclusion: 

Theorem 1. Let d(n) denotes the Dirichlet divisor function, then we have the identity 


- 1 _ 3 m- d(m!) 

24 Onn 24 Gat DE 
Let Qi(2n — 1) denotes the smallest positive odd number such that Qi(2n — 1) - (2n — 1) 
is a two factorial number. That is, Qi(2n — 1)-(2n—1) = (2m — 1)!!, where (2m — 1)!! = 
1x3x5x--++x (2m—1). Then for the sequence {Q1(2n — 1)}, we can also get the following: 
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Theorem 2. Let d(n) denote the Dirichlet divisor function, then we also have the identity 


cs i ((2m — 1)!!) 


2 OiQn= 1) Gna) ee 2 — 








§2. Proof of the theorems 


In this section, we shall use the elementary method to complete the proof of our theorems 
directly. For any positive integer n, let Q(n) = k, from the definition of Q(n) we know that 
k is the smallest positive integer such that n-k = m!. So k must be a divisor of m!. This 
implies that k | m!, and the number of all & (such that & | m!) is d(m!). On the other hand, if 
k | (m—1)!, then Q(n) # k, and the number of all k& (such that & | (m—1)!) is d((m—1)!). So 
the number of all k (such that k | m! and k + (m—1)!) is d(m!) — d((m—1)!). That is means, 
the number of all & (such that Q(n) = & and Q(n)-n = m! ) is d(m!) — d((m— 1)!). From this 


we may immediately get 





= 1 = 1 T8 H{[Q(n)-n=m! 
ae 3 4 _ $F Hin) n= mt 








eats m-—1)! te d(m! te d(m! 
si25- Bay ey ee 


m! m-+1)! 


m=2 m=1 m=1 





TS (m+1)-d(m!) —d(m!) — SS m-d(m!) 
- >! =e 


a m+ 1)!’ 


m=1 m=1 


where d(m) is the Dirichlet divisor function, and #{Q(n) -n = m!} denotes the number of all 
solutions of the equation Q(n)-n = m!. This proves Theorem 1. 
Similarly, we can also deduce Theorem 2. This completes the proof of Theorems. 
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Abstract Murthy [1] introduced the Smarandache Bisymmetric Determinant Natural Se- 


quence. In this paper, we derive the sum of the first n terms of the sequence. 


Keywords The Smarandache bisymmetric determinant natural sequence, the n-th term, 


the sum of the first n terms. 


81. Introduction 


The Smarandache bisymmetric determinant natural sequence (SBDNS), introduced by 
Murthy [1], is defined as follows. 
Definition 1.1. The Smarandache bisymmetric determinant natural sequence, {SBDNS(n)}, 


12 3 4 
1 2 3 
1 2 2 3 4 3 
1], | aoe Se Fells ; 
2 3. 4 3 2 
3.2 1 
4 3 2 1 


A first few terms of the sequence are 1, -3, -8, 20, 48, -112, -256, 576, ---. 
The following result is due to Majumdar [2]. 
Theorem 1.1. Let a, be the n-th term of the Smarandache bisymmetric determinant 


natural sequence. Then, 


1 2 3 ss m—-2 n-il n 
2 3 4 se onal n n—-1 
n 
n-2 n-1l n 5 4 3 
n-1 n n—-1 4 3 2 
n n-1 n-2 3 2 1 
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Let {S;,} be the sequence of n — th partial sums of the sequence {a,,}, so that 


a ee 


k=1 


This paper gives explicit expressions for the sequence {S,,}. This is given in Theorem 3.1 in 
Section 3. In Section 2, we give some preliminary results that would be necessary for the proof 
of the theorem. We conclude this paper with some remarks in the final section, Section 4. 


§2. Some preliminary results 


In this section, we derive some preliminary results that would be necessary in deriving the 
expressions of S;, in the next section. These are given in the following two lemmas. 


Lemma 2.1. For any integer m > 1, 


S- 92(k-1) _ 1 (gam _ Tg 


15 
k=1,3,...,(2m—1) 
Proof. Since 


S- QUT) — 1 4 244 28 +. 4 oAlm—D) 
k=1,3,...,(2m—1) 


is a geometric series with common ratio 24, the result follows. 


Lemma 2.2. For any integer m > 1, 


1 
225 





1 
x ke 220k-1) — Bm jo? (perl 2 1%), 


Proof. Denoting by s the series on the left above, we see that 








PS 1B oD ears Om 1 ee, (*) 
so that, multiplying throughout by 2+, we get 
Oa = ee eo oe i = 3) Os Oe De (+x) 


Now, subtracting (**) from (*), we have 
(1—24%)s=142-24 [a+ ott ----4.240"—2) — (2m — 1) -24 


94(m—1) _ 4 
=1428-{ \ (2m — 1). 24” 





24-1 


= = (2441 — 17) — (2m —1)-24™" 


which now gives the desired result. 
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§3. Main results 


In this section, we derive the explicit expressions of the n-th partial sums, S,, of the 
Smarandache bisymmetric determinant natural sequence. 


From Theorem 1.1, we see that, for any integer & > 1, 


a2k + G2k41 = (—1)*(6k + SV oleae 


Qop+2 + Gops3 = (—1)*** (6k + 11)2*, 


so that 
Gop + don41 + Gor+2 + doe+3 = 3(—1)"*1(6k + 13)224-9), (1) 


Letting 
Sp = a, + a2 +-+++4n, 


we can prove the following result. 


Theorem 3.1. For any integer m > 0, 








1). Sami = 57m gry de ll 7 = 7 {(60m + 31)2*"~* — 3}; 

2). Samia = =; mm 24mts — = POTN) ~ - - (10m + 11)24"*7 + 3}; 
3): Saunas = -: is AD) = _ Yaoi ~ = -4 { (60m + 61)24"*? + 31. 
4). Sam44 = : mm gamers + = Yai ~ = 7 { (5m + giorerh) 3h. 


Proof. To prove the theorem, we make use of Lemma 2.1 and Lemma 2.2, as well as 
Theorem 1.1. 


1). Since S441 can be written as 


Sam41 = 41 +2 +°++°+ 44,4, = a1 + s (dak + Gor41 + G2r42 + G2k43); 
k=1,3,...,(2m—1) 





by virtue of (1), 


Sam41 = a1 +3 S- (1) 6+ 192 


=0,+326 Ss k22*-)) 4.13 S- 92(k—-1) 
k=1,3,...,(2m—1) k=1,3,...,(2m—1) 


Now, appealing to Lemma 2.1 and Lemma 2.2, we get 


2 13 
oe Qe 17) 


: (oe 1) 








Sam+1 =1 t [scm ij2”" 
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which now gives the desired result after some algebraic manipulations. 
2). Since 


Sam+2 = S4m+1 + A4m+25 


from part (1) above, together with Theorem 1.1, we get 





3 31 6 
ee =. Ee. : 92(2m+1) : 94m 4 . 94m 
Sam+2 E oe 25 ag ree 
which gives the desired expression for S442 after algebraic simplifications. 
3). Since 


Ln games _ 


11 
Sam+3 = Sam+2 + A4m4+3 = |-5 2 


6 
bearer 92(2m+1) 3 le 4 4 . 9g4m+1 


we get the desired expression for S4,,43 after simplifications. 
4). Since 


6 ee A eee 6 ” 
Sam+4 = Sam+3 + G4am44 = [-$-m-24 ie al +) F] — an5)-2! a? 


the result follows after some algebraic simplifications. 
The case when m = 0 can easily be verified. 


Hence, the proof is complete. 


§4. Remarks 


Theorem 3.1 in the previous section gives the n-th term of the sequence of partial sums, 
{S,,}, in all the possible four cases. The following lemmas prove that, in each case, S,, is indeed 
an integer. 

Lemma 4.1. For any integer m > 0, 24”~!(60m + 31) — 3 is divisible by 25. 

Proof. The result is true for m = 0,1. So, we assume its validity for some positive integer 


m. Now, since 
[2443 (60(m + 1) +31} — 3] — {24-1 (60m + 31) — 3} = 25(36m 4+ 57)24"-1, 


it follows, by virtue of the induction hypothesis, that 24*?{60(m-+ 1) +31} —3 is also divisible 
by 25. Thus, the result is true for m+ 1 as well, completing induction. 

Lemma 4.2. For any integer m > 0, 24”*1(10m + 11) +3 is divisible by 25. 

Proof. is by induction on m. The result is clearly true for m = 0,1. Now, assuming its 


validity for some positive integer m, since 





[24*5 {10(m + 1) + 11} + 3] — {244 (10m + 11) + 3} = 25(6m + 13)20"41, 


is divisible by 25, it follows that the result is true for m+ 1 as well. This completes the proof. 
Lemma 4.3. For any integer m > 0, 24*1(60m + 61) + 3 is divisible by 25. 
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Proof. is by induction on m. The result is clearly true for m = 0,1. We now assume that 


the result is true for some positive integer m. Then, since 





[24+5{60(m + 1) + 61} + 3] — {24"+1(60m + 61) + 3} = 25(36m + 75)24"+1, 


this, together with the induction hypothesis, shows that the result is true for m + 1 as well. 
This completes the proof by induction. 

Lemma 4.4. For any integer m > 0, 240+) (5m + 32) — 3 is divisible by 25. 

Proof. We first assume that the result is true for some positive integer m. Now, since 





(ater) sa + 1) +13} — 3) — {2 (Gm + 13) — 3} = 258m +8)2™», 


this together with the induction hypothesis, shows that 24(+?){5(m + 1) + 13} — 3 is also 
divisible by 25. This, in turn, shows that the result is true for m-+ 1 as well. To complete, we 
have to prove the validity of the result for m = 0,1, which can easily be checked. 

The Smarandache bisymmetric arithmetic determinant sequence, introduced by Murthy 


[1], is 
a+d a+d2d 
a, 7 ae! a+d a+2d a+d}), --. 
ore a+2d a+d a 





The n-th term of the above sequence has been found by Majumdar [2] to be 


n 


lal («+ 1a) (24)""?. 


Note that the Smarandache bisymmetric determinant natural sequence is a particular case of 
the Smarandache bisymmetric arithmetic determinant sequence when a = 1 and d= 1. 
Open Problem: To find a formula for the sum of the first n terms of the Smarandache 


bisymmetric arithmetic determinant sequence. 
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Abstract For any fixed natural number k > 2 and any positive integer n, we call b,(n) as a 
k-power complements function of n, if b,(n) denotes the smallest positive integer such that 
n- be(n) is a perfect k-power. In this paper, we use the elementary method to study the 


asymptotic properties of b,(n!), and give an interesting asymptotic formula for In bx (n!). 


Keywords k-power complements function, Standard factorization, Asymptotic formula. 


§1. Introduction and results 


For any fixed natural number k > 2 and any positive integer n, we call b(n) as a k-power 
complements function of n, if b,(m) denotes the smallest positive integer such that n- bz (n) is 
a perfect k-power. Especially, we call b2(n), b3(n), b4(m) as the square complements function, 
cubic complements function and quartic complements function respectively. In reference [1], 
Professor F.Smarandache asked us to study the properties of the k-power complements function. 
About this problem, there are many authors had studied it, and obtained a series results. For 
example, in reference [2], Wenpeng Zhang calculated the value of the series 


+00 1 
a (n= be(n))>” 


where s is a complex number with Re(a) > 1, k = 2,3,4, and prove that 


— 1, 2s) 
d, (n= b2(n))> — C(4s) ” 





n=1 


where ¢(s) is the Riemann-zeta function. 
Ruiqin Fu [3] discussed the asymptotic properties of In bg(n!), and proved that 


In bo(n!) = nIn2+O[n-exp oo , 
(Inlnn)s 


where n > 2 is an positive integer, A > 0 is a constant. 

But for the properties of the k-power complements function of n, we still know very little 
at present. In this paper, we use the elementary method to study the asymptotic properties 
of b,(n!), and give an interesting asymptotic formula for In b;,(n!). That is, we shall prove the 


following general conclusion: 
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Theorem. For any fixed natural number & > 2 and any positive integer n, we have the 
asymptotic formula 


Co 


1 —Aln? n 
In by. (n!) = k—- —— O : —— 
nriod=a(t-F apazy)+° (e-em Set) 


where A > 0 is a constant, and exp(y) = e¥. 


§2. Two simple lemmas 


To complete the proof of the theorem, we need the following two simple lemmas: 
Lemma 1. Let n! = p{' pS? --- p%s denotes the factorization of n! into prime powers. Then 
we have the calculate formula 


by (n!) bk (pt* )bk (po?) +++ be (PS) 


d d 
po” (p1) - per (p2) pr dPs) 


l| 


where the ord function is defined as: 


k-1, ifa;=km+1, 


k-2, ifa;=km+2, 
ord(p;) = 





0, ifa; =k(m+1). 


where m = 0,1,2,-:-. 
Proof. See reference [3]. 


Lemma 2. For any real number x > 2, we have the asymptotic formula 


—Aln? x 
O(x) =) mp=2+0 (0 (ae*)) ; 


psa 


where exp(y) = e”, and A > 0 is a constant. 
Proof. See reference [4] and [5]. 


§3. Proof of the theorem 


In this section, we shall complete the proof of Theorem. First from Lemma 1 we have 


In b(n!) = Im (be (py) be (D2?) - + - be (DS*)) 
= S- ord(p) In p 
pen 
ktord(p) 
= (k-1) S- Inp+ (k — 2) S- Inp+-++++ ye Inp 
5 <pgn 3 SPL EPL ER 





+ (k-1) >> Inp+---+ SO Inpt+---+O(1). 


Eps SPS ET ah <P Et 
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Let n be a positive large enough, if a prime factor p of n! in the interval G n|, then the power 
of p is 1 in the standard factorization of n!. Similarly, if a prime factor p of n! in the interval 
(5. os then the power of p is 2 in the standard factorization of n!; If a prime factor p of n! in 
the interval es =, then the power of p is 3 in the standard factorization of n!,---. 


Then from the above formula and Lemma 2 we have 

(1) (0(m) -0(5)) + 2 (0(5) -9(5)) +--+ (0 (Gz) -0(G)) 

+ «-1)(0(,2)) 6(ce5)) ++) (x) ++ 00) 
( oe i))teni( (scitepat ta) 


In by (n!) 














II 
3 
aN 
> 
| 
a 
the 
SS 
> 
a —_ 
+ 
a 
Sie eae” 
+ 
) 
aN 
3 
fa) 
ta 
Ko} 
aN 
Ba 
= — 
S.\"5 
ar! 3 
YY 
YS 


where A > 0 is a constant, and exp(y) = e¥. 
This completes the proof of Theorem. Absolutely, this theorem can be used to study 


k-power residues numbers. 
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Abstract This paper gives an alternative approach to the solution of the LP problem with 
equality constraints. In this scheme, the constraint system of linear equations is first solved 
by elementary row operations (if such a solution exists). In the second stage, the solution 


of the system of linear equations is directly plugged in the objective function and then it is 


minimized. This scheme thus avoids the use of the simplex method. 


Keywords Linear programming, equality constraint, elementary row operations. 


81. Introduction 


In a recent paper, Ru, Shen and Xue [1] considered the problem of finding an initial basic 
feasible solution (bfs) of the LP problem of the form 


minz=C'X st: AX=B, X > 0, (1) 
where 
ai Q12°** Qin 
C1 a21 a22°** Gan by Ly 
C2 : be Z2 
C = ,A —4 ,B = xX = 
ail Qj2°°* Ain 
om ; bm Tn, 


Gm1 Am2°*** Amn 


In (1), without loss of generality, we may assume that B > 0; for otherwise, if some b; < 0 
(1 <i <™m), then multiplying the i-th constraint throughout by (-1), we get one with (—b;) > 0. 

To find an initial bfs of the LP problem (1), the procedure suggested by Ru et al. [1] is as 
follows : 

Step 1: Starting with the augmented matrix (A | B), reduce it to a form with an identity 
matrix of desired order, by elementary row operations only. 
Note that, the elementary row operations applicable in this case are 

(a) multiplying any row of (A | B) by a positive constant, 
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(b) multiplying any row of (A| B) by a constant and adding it to a second row. 

Throughout this paper, the first row operation would be denoted by R; — kR; (which 
means that the i-th row R; is multiplied by the constant k, where k > 0), and the second row 
operation would be denoted by R; — R;+kR,; (which means that the j-th row R; is multiplied 
by the constant k and the resulting row is added, term-by-term, to the 7-th row R; to get the 
new i-th row). Note that the latter row operation is allowed only if b; + kb; > 0. 

If rank(A) = rank(A | B) = I(< m), then using the above elementary row operations, it 
is possible to find a unit matrix of order 1, which may be exploited to find an initial bfs of the 
LP problem (1). It may be mentioned here that, if 1 << _m, then only / of the m system of linear 
equations AX = B are independent, and the remaining / — m equations can be ignored. It may 
also be mentioned here that, if rank(A) 4 rank(A | B), then the system of linear equations 
AX = B is inconsistent, and as such, the LP problem (1) is infeasible. 

Step 2: With the initial bfs found in Step 1 above, form the initial simplex tableau in the 
usual manner. 

Step 3: Proceed in the usual way of the simplex method to find an optimal solution. 

The method outlined above is nevertheless not a new one. In fact, the same procedure has 
been followed by Papadimitriou and Steiglitz [2]. 


In the next section, we propose an alternative approach. 


§2. An alternative approach 


In this section, we suggest an alternative approach to the solution of the LP problem (1). 
The procedure is quite simple, and can be done in the following two steps. 

Step 1: Reduce the system of linear equations AX = B to a form containing the unit 
matrix of requisite order, using elementary row operations only. 

Step 2: With the solution found in Step 1, consider the problem of minimizing z = CtX. 

Thus, the method works with the solution of the system of linear equations AX = B and 
the objective function only, avoiding the simplex method completely. 

We illustrate our method with the help of the following examples, due to Ru et al. [1]. 

Example 2.1. Consider the following LP problem : 





min z = 4%, + 322 8.t. 
1 1 2 

tT +22 +~2£3 ta =2 
3 ¢ 3 
971 “gta = 3 
32, —b62r2 +4¢7,  =0 


U1, 02,3, 04 2 0. 


To solve the above LP problem, we first solve the constraint system of linear equations, 
using elementary row operations only. We thus start with the augmented matrix (A | B) and 
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proceed as follows, using the indicated row operations: 





bor og -B 2) pip -R 2 
(A|B)=] 2 oO -2 0 3 | Rs > R3+6R,| 3 -1 9 3 
3-6 0 4 0 6 0 3 0 12 
es 210 -2 5\ R3;—>iR; P11 4 

Ry — Ry + Re 4 
> 0 0 0 5 Ri — Ri — eRe 2 0 0 0 5 

Rz — Ro + §Ro 4 
6 0 3 0 12 fig — + fig Ap 0 0 1 0 0 

The final tableau gives the following solution : 








2 2 
v =2,%3=0,%2=1+ 34 =1+ 3h (t > 0 is any real number), 


and the objective function is z = 4%, + 3%. = 84+ 3(14+ af =11+ a 
Since the objective is to minimize z, it is clear that z is minimized when t = 0, and the optimal 
solution of the given LP problem is thus xj = 2,275 = 1,23 = 0,24 = 0; min z* = 11. 

Example 2.2. To solve the LP problem 











min z = —34, + %9+ 43 st. 

Ly 222 @3 +24 =11 
—Ax, +2X2 +2x3 5 = 3 
—221 TL => 1 





T1,%2,%3,04,05 2 0, 


we first solve the constraint system of linear equations by elementary row operations. This is 
done below. 


2.47 ft Dt ae ae ee ee 
(A|B)=|] -4 120 -1 3 |R—R+2R.| 4120 -1 3 
<5 =f 20 © 4 20 io © 4 


Ry, — RR, +5R3 0100 -1 1 
—2 0 1 0 0 1 


The final tableau is equivalent to the following system of linear equations 


321 t+a4 —2¢45 =12 
r2 —£5 =2 
291 +23 =1 


with the solution 








vy = 1,85 = te, vo = 14te, x3 = 142t), 24 = 12—3t14-2tg; (t1 > 0,t2 > 0 are any real numbers), 
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and the objective function is z = 2 — t) + to. 

Since the objective is to minimize z, it is evident that t2 should be as small as possible and t, 
should be as large as possible. Now, since x4 = 12 — 3t, + 2t2 > 0, we see that, we must have 
to = 0, t; = 4. Hence, the desired optimal solution is 





x} =4,05 = 1,073 = 9,05 =0, 2% = 0; min z* = —3. 





Example 2.3. To solve the LP problem 





min z = 42, + 323 S.t. 
2 

2 +2%2 "3° ~3%4 = 2 

371 14% =3 

32, —b62r2 +24 =0 


1, 02,3, 04 2 0, 


we start with the augmented matrix (A|B) and proceed as follows, using the indicated row 
operations : 


a Dg A eg 
(A| B) = 3 0 3 ; 3 | R3 — R3+ 6R, 3 0 3 : 3 
D) 4 D) 4 
3 6 0 40 603 01 
La 2 oe i e442 22 4 
aS , 4 3 R, — RR, - =Re 4 3 
R3—>R3-4Ri| 2 09 2 9 3 oe ios 22 
2 4 Ro —> =Ro 2 
00 0 00 3 00 1 0 0 


The final tableau is equivalent to the following system of linear equations 





with the solution 


1 1 2 
%3 = 1,0, =2— gti @4 =to,% =1- qa + 3f23 (ti; > 0,t2 > 0 are any real numbers), 


and the objective function is z= 8+ ty. 


Since the objective is to minimize z, we must have t; = 0. Hence, the desired solution is 
gy = 2,25 =1+ gts = 0,23, = t(t > 0 is any real number); min z* = 8. 


The above solution shows that the given LP problem has infinite number of solutions. 
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§3. Some observations and remarks 


In the previous section, we have illustrated our method with the help of three simple 
examples. However, in large-scale problems, the success as well as the computational efficiency 
of the method remains to be checked out. 

Clearly, the method is useful when the correct basis has been found. This point is illustrated 
with the help of the following example, due to Papadimitriou and Steiglitz [2]. 

Example 3.1. Consider the LP problem below : 


ming =27,+%2+%3+%4+ 25 8.0. 
34, +2% +23 =] 
524 +%Qq +23 +24 =3 
221 +522 TL3 +25 = 4 








U1, 02,%3,%4, 05 2 0. 


We first solve the constraint system of linear equations by elementary row operations, using the 
indicated row operations. 





321001) R,—4R, 24. 2 oo 2 
(A|B)=) 5 1110 3 +) R—R-R fo ¢10 3 
25 101 4/ R3—R3-5Ri \ -# 0 -2 01 3 
The final tableau is equivalent to the following system of linear equations 
3 1 1 
i +22 +273 = 2 
571 +573 +24 = 
2 3 
igs ie = Sg 
with the solution 
Ly = t1, 23 = to(ti > 0,t2 > 0 are any real numbers), 
1 - 1, 5 C i eis i, ies 
=>-— — _-_ =>--—- — _ Ce => = — — 2 
r2 2. 91 204 2 gt 9°25 = 5 gq 1 T 9° 


9 3 
and the objective function is z = = + —(ti + t2). 
Since the objective is to minimize z, it is evident that we must have t, = 0, t2 = 0, giving the 
desired optimal solution 





xi =0,25 





* «9 ,_ 3 
93 0,24 97% 2? 


On the other hand, the solution below 


3 2 100 1 es 3 21 0 0 41 
Ry — Ro - Ry 

(A|B)= 5 1 1 10 3 2 -1 0 1 0 2 
R3 — R3—-— Ry 

25 10414 —1 3.0 0 1 3 


94 A.A.K. Majumdar No. 3 





leads to the system of equations 





32, +2%2 +23 =| 
224 —2X2 +2%4 =2 
TX +329 +25 =3 


whose solution is 
Ly = t1, Lo = to(ti > 0,t2 > 0 are any real numbers), 


v3 = 1—- 3ty — 2te,%4 = 2 — 241 + t2,%5 =3 +t, — 3t9; 


and the objective function is z = 3[2 — (t; + t)]. 
Thus, the problem of finding the solution of the given LP problem reduces to the problem of 
solving the LP problem 


max t, + to S.t. 
3t, +2t2 <1 
2t; —tzg <2 
—t; +38t2 <3 





t1,t2 > 0. 


1 
The optimal solution of the above LP problem is t, = 0, te = 5 
Hence, the optimal solution of the original LP problem is 
» 9 


ca . j _— 
Us, se = 








xi =0,25 ae 0,24 5 

The above example shows that, in the alternative approach proposed in this paper, as well 
as the method suggested by Ru et al. [1], the choice of the unit matrix (of desired order) during 
the solution of the system of linear equations AX = B (by elementary row operations) may 
save the number of iterations considerably. However, no particular criterion is specified in this 
respect. It thus still remains to check how far these methods can compete with the existing 
method. 
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Abstract Two triangles T(a, b,c) and T(a, b F c) are said to be pseudo Smarandache related 
if Z(a) = Z(a), Z(b) = Z(b), Z(c) = Zc), where Z(.) is the pseudo Smarandache function, 
and T’(a,b,c) denotes the triangle with sides a,b and c. This paper proves the existence of 
an infinite family of pairs of dissimilar Pythagorean triangles that are pseudo Smarandache 
related. 
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81. Introduction and result 


The Smarandache function, denoted by S(n), is defined as follows. 
Definition 1.1. For any integer n > 1, ( Z* being the set of all positive integers), 


S(n) =min{m:me Zt, nimi}. 


The following definition is due to Sastry [6]. 
Definition 1.2. Two triangles T(a, b,c) (with sides of length a,b and c) and T(a’,b,c ) 
(with sides of length a,b’ and c’ ), are said to be Smarandache related if 


Definition 1.3. A triangle T(a,b,c) is said to be Pythagorean if and only if one of its 
angles is 90°. 

Thus, a triangle T(a,b,c) (with sides of length a,b and c) is Pythagorean if and only if 
a? + b? = c?. 

Sastry|6] raised the following question : Are there two distinct dissimilar Pythagorean 
triangles that are Smarandache related? Recall that two triangles T(a, b,c) and T(a’ ,b ,c’) are 
similar if and only if the corresponding three sides are proportional, that is, if and only if 


Otherwise, the two triangles are dissimilar. 
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The following result, due to Ashbacher [1], answers the question in the affirmative. 

Theorem 1.1. There are an infinite family of pairs of dissimilar Pythagorean triangles 
that are Smarandache related. 

The proof of Theorem 1.1 is rather simple : For any prime p > 17, the two families 
of dissimilar Pythagorean triangles T (3p, 4p,5p) and T(5p, 12p, 13p) are Smarandache related, 
since S(3p) = S(4p) = S(5p) = S(12p) = S(13p) = p. 

Ashbacher [1] introduced the concept of pseudo Smarandache related triangles, defined as 








follows. 

Definition 1.4. Two triangles T(a, b,c) and T(a’,b',c ) are said to be pseudo Smarandache 
related if Z(a) = Z(a_), Z(b) = Z(b), Z(c) = Zc). 

In Definition 1.4 above, Z(.) denotes the pseudo Smarandache function. Recall that the 
pseudo Smarandache function is defined as follows. 


Definition 1.5. For any integer n > 1, Z(n) is the smallest positive integer m such that 
1 
142t--+m= mn + | ag. divisible byr: Thus, 








(m+ 1) 
2. 


Z(n) = minf{fm:m Zt nl in > 1. 


Ashbacher [1] used a computer program to search for dissimilar pairs of Pythagorean 
triangles that are pseudo Smarandache related. He reports some of them, and conjectures 
that there are an infinite number of pairs of dissimilar Pythagorean triangles that are pseudo 
Smarandache related. 

In this paper, we prove the conjecture of Ashbacher [1] in the affirmative. This is done 
in Theorem 3.1 in Section 3. We proceed on the same line as that followed by Ashbacher. 
However, in the case of the pseudo Smarandache related triangles, the proof is a little bit more 
complicated. The intermediate results, needed for the proof of Theorem 3.1, are given in the 
next Section 2. 


§2. Some preliminary results 


The following result, giving the explicit expressions of Z(3p), Z(4p) and Z(5p), are given 
in Majumdar [5]. 
Lemma 2.1. If p> 5 is a prime, then 


Z(3p) = 
oP) D, if 3\(p + 1). 


Lemma 2.2. If p> 5 is a prime, then 


5 


p= 1, if 8|(p — 1) 
D, if 8|(p+ 1) 
3p — 1, if 8|(3p—1), 
3p, if 8|(3p +1) 


5 
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Lemma 2.3. If p> 7 is a prime, then 





p=, if 10|(p—1), 

Z(5p) = P; if 10|(p +1), 
do — 1, if 5|(2Qp—1), 

2p, if 5|(2p+ 1). 


The explicit expressions of 7(12p) and Z(13p) are given in the following two lemmas. 


Lemma 2.4. If p > 13 is a prime, then 





p—-1, if 24|(p — 1) 

DP, if 24|(p + 1), 

35 = 1, if 8|(3p — 1), 

Z(12p) = p |(3p — 1) 
3D, if 8|(3p + 1) 

7p—1, if 24|(7p—1), 

7p, if 24|(7p+ 1). 


Proof. By definition, 


(m+ 1) 


1 
Z(12p) = min{m : 12p|— 5 ss 


} = min{m : p| “ST }. (1) 


If plm(m + 1), then p must divide either m or m+ 1, but not both, and then 24 must divide 
either m+ 1 or m respectively. In the particular case when 24 divides p — 1 or p+ 1, the 
minimum m in (1) may be taken as p—1 or p respectively. We now consider the following eight 
cases that may arise : 
Case 1: pis of the form p= 24a + 1 for some integer a > 1. 
In this case, 24|(p — 1). Therefore, 7(12p) = p—1. 
Case 2: p is of the form p = 24a + 23 for some integer a > 1. 
Here, 24|(p+ 1), and so, Z(12p) = p. 
Case 3: p is of the form p = 24a + 5 for some integer a > 1, so that 8|(3p + 1). 
In this case, the minimum m in (1) may be taken as 3p. That is, Z7(12p) = 3p. 
Case 4: p is of the form p = 24a + 19 for some integer a > 1. 
Here, 8|(3p — 1), and hence, Z7(12p) = 3p — 1. 
Case 5: p is of the form p = 24a +7 for some integer a > 1. 
In this case, 24|(7p — 1), and hence, Z(12p) = 7p. 
Case 6: p is of the form p = 24a +17 for some integer a > 1. 
Here, 24|(7p +1), and hence, Z7(12p) = 7p+ 1. 
Case 7: p is of the form p = 24a +11 for some integer a > 1. 
In this case, 8|(3p — 1), and hence, Z(12p) = 3p - 1. 
Case 8: p is of the form p = 24a + 13 for some integer a > 1. 
Here, 8|(3p +1), and hence, Z7(12p) = 3p. 
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Lemma 2.5. For any prime p > 17, 











p-, if 13|(p— 1), 

P, if 13\(p + 1), 

2p —1, if 13|(2p — 1), 

2p, if 13|(2p +1), 

3p — 1, if 13|(3p — 1), 

Hike 3p, if 13|(3p +1), 

4p — 1, if 13|(4p — 1), 

Ap, if 13|(4p +1), 

5p —1, if 13|(5p — 1), 

5p, if 13|(5p +1), 

6p — 1, if 13|(6p — 1), 

6p, if 13|(6p + 1). 

Proof. By definition, 
Z(13p) = min{m : 13pm) = min{m pL DY, (2) 


We have to consider the twelve possible cases that may arise : 

Case 1: p is of the form p = 13a +1 for some integer a > 1. 
In this case, 13](p — 1), and so, Z(13p) = p— 1. 

Case 2: p is of the form p = 13a +12 for some integer a > 1. 
Here, 13|(p+ 1), and hence, Z7(13p) = p. 

Case 3: p is of the form p = 13a + 2 for some integer a > 1. 
In this case, 13|(6p + 1), and hence, Z(13p) = 6p. 

Case 4: p is of the form p = 13a +11 for some integer a > 1. 
Here, 13|(6p — 1), and hence, Z(13p) = 6p — 1. 

Case 5: p is of the form p = 13a +3 for some integer a > 1. 
In this case, 13|(4p + 1), and hence, Z(13p) = 4p. 

Case 6: p is of the form p = 13a + 10 for some integer a > 1. 
Here, 13|(4p — 1), and hence, Z(13p) = 4p — 1. 

Case 7: p is of the form p = 13a +4 for some integer a > 1. 
In this case, 13|(3p +1), and hence, Z(13p) = 3p. 

Case 8: p is of the form p = 13a +9 for some integer a > 1. 
Here, 13|(3p — 1), and hence, Z(13p) = 3p — 1. 

Case 9: p is of the form p = 13a +5 for some integer a > 1. 
In this case, 13|(5p +1), and hence, Z(13p) = 5p. 

Case 10: p is of the form p = 13a + 8 for some integer a > 1. 
Here, 13|(5p — 1), and hence, Z(13p) = 5p — 1. 

Case 11: p is of the form p = 13a + 6 for some integer a > 1. 
In this case, 13|(2p +1), and hence, Z(13p) = 2p. 
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Case 12: p is of the form p = 13a+7 for some integer a > 1. 
Here, 13|(2p — 1), and hence, Z7(13p) = 2p — 1. 


§3. Main result 


We are now state and prove the main result of this paper in the following theorem. 

Theorem 3.1. There are an infinite number of pairs of dissimilar Pythagorean triangles 
that are pseudo Smarandache related. 

Proof. We consider the pair of dissimilar Pythagorean triangles 


T(3p,4p,5p) and T(5p, 12p, 13p), (3) 
where p is a prime of the form 
p= (2°-3-5-13)n+1=1560n+1,neZ*. (4) 


By Lemma 2.1 - Lemma 2.5, 





Z(3p) = Z(4p) = Z(5p) = Z(12p) = Z(13p) = p — 1, 





so that the triangles T(3p,4p,5p) and T'(5p,12p,13p) are pseudo Smarandache related. Now, 
since there are an infinite number of primes of the form (4) (by Dirichlet’s Theorem, see, for 
example, Hardy and Wright [3], Theorem 15, pp. 13), we get the desired infinite number of 
pairs of dissimilar Pythagorean triangles that are pseudo Smarandache related. 

It may be mentioned here that, the family of pairs of triangles (3), where p is a prime of 
the form 


p = 1560n —1,ne Zt, (5) 


also forms (dissimilar) pseudo Smarandache related Pythagorean triangles. 
§4. Some remarks 


The pseudo Smarandache function Z(n) is clearly not bijective. However, we can define 


the inverse Z~!(m) as follows : 
Z'(m) ={n€Z*+:Z(n) =m} for any integer m> 3, (6) 


with 
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As has been pointed out by Majumdar [5], for any m € Z*, the set Z~!(m) is non-empty 
1 
aid tounier gi 


following two conditions are satisfied : 


(1) n divides mine) ; 


(2) n does not divide las Yee any with 1<€<m-1. 


We can look at (6) from a different point of view : On the set Z* , we define the relation 


as its largest element. Clearly, n € Z~!(m) if and only if the 





R as follows : 
For any n1,n2€Z+t,ni Rng ifandonlyif Z(n1) = Z(nz). (8) 


It is then straightforward to verify that # is an equivalence relation on Zt . It is well- 
known that an equivalence relation induces a partition (on the set Zt ) (see, for example, Gioia 
[2], Theorem 11.2, pp. 32). The sets Z~!(m),m € Z* , are, in fact, the equivalence classes 
induced by the equivalence relation # on Z* , and possesses the following two properties : 

(1) Ss 2M) = 2+ 

(2) Z-2(nma) 9 Z-1(ana) = Of ms f mo. 

Thus, for any n € Z*, there is one and only one m € Z* such that n € Z~!(m). 

A different way of relating two triangles has been considered by Ashbacher [1], which is 
given in the following definition. 

Definition 4.1. Given two triangles, T(a, b,c) and T(a ,b',c ), where 


at+b+c=180=a +b +c, (9) 


they are said to be pseudo Smarandache related if 


The difference in Definition 1.4 and Definition 4.1 is that, in the former, the sides of the 
pair of triangles are pseudo Smarandache related, while their angles, measured in degrees, are 
pseudo Smarandache related in the latter case. Note that the condition (9) merely states the 
fact the sum of the three angles of a triangle is 180 degrees. 

Using a computer program, Ashbacher searched for pseudo Smarandache related pairs of 
triangles (in the sense of Definition 4.1). He reports three such pairs. 

However, in this case, the equivalence classes Z~!(m),m € Z*, might be of some help. 
The condition (9) can be dealt with by considering the restricted sets Z~(m|z) : 


Z- (mix) = {ne Z+ : Z(n) =m,1 <n < 178}. (10) 


Table 4.1 gives such restricted sets related to our problem. 
Clearly, two pseudo Smarandache related triangles T(a, b,c) and T(a ,b,c) must satisfy 
the following condition : 


aa € Z*(my|7); 0,0 € Z~+(ma|r), ¢, ce Z\(m3|m) for some m1,m2,m3 € ZT. 
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Thus, for example, choosing 


a,a € {2,6} = Z-1(3|r);b,b © {44, 48, 88, 132,176} = Z~!(32|r), 


we can construct the pseudo Smarandache related triangles T(2,48,130) and T(6, 44, 130). 
Again, choosing a,a’,b, b! € {25,50, 75, 100,150} = Z~1(24|7), we can form the pair of pseudo 
Smarandache related triangles T(25, 150,5) and 775, 100,5) with the characteristic that 








Choosing a,a ,b,b € {8, 20, 24, 30, 40,60, 120} = Z—1(15|7), we see that the equilateral triangle 
T (60, 60, 60) is pseudo Smarandache related to the triangle T(20, 40, 120)! 


Ashbacher [1] cites the triangles T(4, 16, 160) and T(14, 62, 104) as an example of a pseudo 
Smarandache related pair where all the six angles are different. We have found three more, 


given below : 
(1) T(8, 12,160) and T(40, 36, 104), 
(with Z7(8) = Z(40) = 15, Z7(12) = Z7(36) = 8, Z(160) = 7(104) = 64), 
(2) T(16, 20,144) and T(124, 24, 32), 
(with 7(16) = 7(124) = 31, 7(20) = 7(24) = 15, 7(144) = Z(32) = 63), 
(3) T(37, 50,93) and T(74, 75, 31), 


(with Z(37) = Z(74) = 36, Z(50) = Z(75) = 24, Z(93) = Z(31) = 30). 


Table 4.1. Values of Z~!(m|r) = {n € Zt : Z(n) =m,1 <n < 180} 
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m | Z~1!(m|r) m | Z—-1(m|r) m | Z~-+(m|r) 
1 | {1} 35 | {90} 80 | {81, 108, 162, 180} 
2 | {3} 36 | (37,74, 111} 82 | {83} 

3 | {2,6} 39 | {52, 130, 156} 83 | {166} 

4 | {5,10} 40 | 41,82, 164} 84 | {170} 

5 | {15} 41 | {123} 87 | {116,174} 
6 | {7,21} 42 | {43,129} 88 | {89,178} 
7 | {4,14, 28} 43 | {86} 95 | {152} 

8 | {9, 12,18, 36} 44 | {99,110,165} 96 | {97} 

9 | {45} 45 | {115} 100 | £101} 

10 | {11,55} 46 | {47} 102 | {103} 

11 | {22, 33,66} 47 | {94,141} 106 | {107} 

12 | {13, 26,39, 78} 48 | {49,56,84, 98,147,168} | 108 | {109} 

13 | {91} 49 | {175} 111 | {148} 

14 | {35,105} 51 | {102} 112 | {113} 

15 | {8, 20, 24, 30, 40, 60,120} | 52 | {53,106} 120 | {121} 

16 | {17, 34,68, 136} 53 | {159} 124 | {125} 

17 | {51,153} 54 | {135} 126 | {127} 

18 | {19,57,171} 55 | {140,154} 127 | {64} 

19 | {38,95} 56 | {76, 114, 133} 128 | {172} 

20 | {42,70} 58 | {59} 130 | {131} 

21 | {77} 59 | {118,177} 136 | {137} 

22 | {23} 60 | {61,122} 138 | {139} 

23 | {46, 69,92, 138} 63 | {32,72,96, 112, 144} 148 | {149} 

24 | {25, 50,75, 100, 150} 64 | {80, 104, 160} 150 | {151} 

25 | {65} 65 | {143} 156 | {157} 

26 | {27,117} 66 | {67} 162 | {163} 

27 | {54, 63, 126} 67 | {134} 166 | {167} 

28 | {29,58} 69 | {161} 168 | {169} 

29 | {87,145} 70 | {71} 172 | {173} 

30 | {31,938,155} 71 | {142} 178 | {179} 

31 | {16,62, 124} 72 | {73,146} 255 | {128} 

32 | {44, 48, 88, 132, 176} 78 | {79} 

34 | {85,119} 79 | {158} 


























Ashbacher [1] reports that, an exhaustive computer search for pairs of dissimilar search for 
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all pseudo Smarandache related triangles T(a, b,c) and T(a ,b',c’) (a+b+c = 180 =a +b +c) 
with values of a in the range 1 < a < 178, revealed that a cannot take the following values : 


1,15, 23,35, 41, 45, 51, 59, 65, 67, 71, 73, 77, 79, 82, 83, 86, 87, 89, (11) 
90, 91,97, 101, 102, 105, 107, 109, 113, 115, 116, 118, 121, 123, 125, 126, 127, 
131, 134, 135, 137, 139, 141, 142, 143, 148, 149, 151, 152, 153, 157, 158, 159, 

161, 163, 164, 166, 167, 169, 170, 171, 172, 173, 174, 175, 176, 177, 178. 


Can the table of the sets Z~'(m|m), Table 4.1, be utilized in explaining this observation? Of 
course, if a is too large, then there is a possibility that no two dissimilar triangles exist: Very 
large value of a very often forces the two triangles to be similar. For example, if a = 174 € 
Z~'(87|7), then the three possible pairs of values of (b,c) are (1,5), (2,4) and (3,3). Note that 
the two sets Z—!(1|r) = {1} and Z~!(2|r) = {3} are singleton. Thus, if a = 174, then the two 
triangles T(a, b,c) and T(a ,b ,c ) must be similar. On the other hand, if a = 116 € Z~1(87|z), 
then we cannot find b,c with b +c = 64. Now, we consider the case when a = 164. This case 
is summarized in the tabular form below : 






































a=164 a =164 a =82 a’ =41 Remark 
(1,15) | Similar Triangle | Not Possible | Not Possible | Both belong to singleton sets | (b,c) 
(2,14) | Similar Triangle | Not Possible | Not Possible b+e <34 
(3,13) | Similar Triangle | Not Possible | Not Possible | 3 belongs to a singleton set 
(4,12) | Similar Triangle | Not Possible | Not Possible b+c <64 
(5,11) | Similar Triangle | Not Possible | Not Possible b+c <65 
(6,10) | Similar Triangle | Not Possible | Not Possible b+ce <16 
(7,9) | Similar Triangle | Not Possible | Not Possible b+e¢ <57 
(8,8) | Similar Triangle | Not Possible | Not Possible 












































However, if a = 165, we can get two dissimilar pseudo Smarandache related triangles, namely, 
T(165,7,8) and T(81, 21,60). 

A closer look at the values listed in (11) and those in Table 4.1 reveals the following facts: 
(1) The numbers not appearing in any 180 degrees triplets are all belong to singleton sets, with 
the exception of 3 (though we conjecture that 3 cannot appear in any triplet), 47, 64, and 103. 
For the last three cases, we can form the following examples : 

(a) T(47, 76,57) and T(47, 114,19) (with 7(76) = 56 = Z7(114), Z7(57) = 18 = Z(19)), 

(b) T(64, 12,104) and T(64, 36,80) (with 7(12) = 8 = Z(36), 7(104) = 64 = Z(80)), 

(c) T(103, 11,66) and T(103, 55, 22) (with Z(11) = 10 = Z(55), Z7(66) = 11 = Z(22)). 
(2) In most of the cases, if a value does not appear in the triplet, all other values of the corre- 
sponding Z~+(m|z) also do not appear in other triplet, with the exception of 145 € Z~1(29|7) 
(87 € Z~1(29|m) does not appear in any triplet), and 146 € Z~'(72|r)(73 € Z—1(72|7) does 


not appear in any triplet). In this connection, we may mention the following pairs of triplets : 
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(a) T(145, 14,21) and T(145, 28,7) (with 7(14) = 7 = Z(28), Z(21) = 6 = Z(7)), 

(b) T(146, 4,30) and 7(146, 14,20) (with 7(4) = 7 = Z(14), 7(30) = 15 = Z(20)). 

(c) Both 54,63 € Z~1(72|) can appear in one or the other triplet, but 126 € Z~'(72|r) 
cannot appear in any triplet. 

It is also an interesting problem to look for 60 degrees and 120 degrees pseudo Smarandache 
related pairs of triangles, in the sense of Definition 4.1, where a 60 (120) degrees triangle is one 
whose one angle is 60 (120) degrees. We got the following two pairs of pseudo Smarandache 
related 60 degrees triangles : 

(1) T(8, 60, 112) and T(24, 60,96), (with 7(8) = 7(24) = 15, Z(112) = Z(96) = 63), 

(2) T(32,60,88) and T(72,60,48), (with 7(32) = Z(72) = 63,27(88) = Z(48) = 32), 

while our search for a pair of pseudo Smarandache related 120 degrees triangles went in vain. 


We conclude the paper with the following open problems and conjectures. The first two 
have already been posed by Ashbacher [1]. 

Problem 1. Are there an infinite number of pairs of dissimilar 60 degrees triangles that 
are pseudo Smarandache related, in the sense of Definition 1.4? 

Problem 2. Is there an infinite family of pairs of dissimilar 120 degrees triangles that are 
pseudo Smarandache related, in the sense of Definition 1.4? 

Ashbacher [1] reports that a limited search on a computer showed only four pairs of dis- 
similar 60 degrees psedo Smarandache related triangles, while the number is only one in the 
case of 120 degrees triangles. 

In the formulation of Definition 4.1, the number of pairs of dissimilar Smarandache related 
or pseudo Smarandache related triangles is obviously finite. In fact, since the number of integer 
solutions to 


a+b+c= 180 
a>lbel1,ce>1 


is C(179,2) (see, for example, Johnsonbaugh [4], Theorem 4.5, pp. 238), the number of such 
triangles cannot exceed C'(179,2). But this is a very crude estimate. The next problems of 
interest are 

Problem 3. Is it possible to find a tight or better upper limit to the number of pairs of 
dissimilar triangles that are pseudo Smarandache related, in the sense of Definition 4.1? 

Problem 4. Is it possible to find a good upper limit to the number of pairs of dissimilar 
60 degrees triangles that are pseudo Smarandache related, in the sense of Definition 4.1? 

Problem 5. Is it possible to find another pair of 60 degrees pseudo Smarandache related 
triangles, in the sense of Definition 4.1, where all the six angles are acute? 

By a limited search, we found only two pairs of dissimilar 60 degrees pseudo Smarandache 
related triangles, the ones already mentioned above. 

Conjecture 1. There is no pair of dissimilar 120 degrees triangles that are Pseudo Smaran- 
dache related, in the sense of Definition 4.1. 

Conjecture 2. There is no pair of dissimilar 90 degrees triangles that are Pseudo Smaran- 
dache related, in the sense of Definition 4.1. 
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Abstract For any positive integer n, the famous F.Smarandache function S(n) is defined as 
the smallest positive integer m such that n | m!. That is, S(n) = minf{fm: n|m!, n € N}. 
The main purpose of this paper is using the elementary method to study the number of the 
solutions of the congruent equation 19(°—) + 25@-D +...4 (n—1)8(-) +1 =0 ( mod n), 


and give its all prime number solutions. 
Keywords F. Smarandache function, divisibility, primitive root. 


§1. Introduction and result 


For any positive integer n, the famous F. Smarandache function S(n) is defined as the 
smallest positive integer m such that n | m!. That is, S(n) = min{fm: n|m!, n € N}. For 
example, the first few values of S(n) are S(1) = 1, S(2) = 2, S(3) = 3, S(4) = 4, S(5) = 5, 
S(6) = 3, S(7) = 7, S(8) = 4, S(9) = 6, S(10) = 5, S(11) = 11, S(12) = 4, ---. About the 
elementary properties of S(n), many authors had studied it, and obtained a series interesting 
results, see references [1], [2], [3] and [4]. For example, Xu Zhefeng [2] studied the value 
distribution problem of S(n), and proved the following conclusion: 

Let P(n) denotes the largest prime factor of n, then for any real number x > 1, we have 
the asymptotic formula 


 (S(n) - Pw? = 22" 5 (=| , 


2 
nee ln* x 


where ¢(s) is the Riemann zeta-function. 
Lu Yaming [3] studied the solutions of an equation involving the F. Smarandache function 
S(n), and proved that for any positive integer k > 2, the equation 


S(mz + m2 +--+ + mp) = S(m 1) + S(mg) +--+ + S(me) 


has infinite group positive integer solutions (m1, m2, --: , mz). 
Rongji Chen [4] proved that for any fixed r € N with r > 3, the positive integer n is a 


solution of the equation 


S(n)" + S(n)" 1 4---4+ S(n) =n 
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if and only if 
n= pip" +p"? +---+1), 


where p is an odd prime such that 





On the other hand, in reference [5], C. Dumitrescu and V. Seleacu asked us to study the 
solvability of the congruent equation 


(2) 4 Fe dh (ge yr) A =) mod Hi, (1) 


About this problem, it seems that none had studied it yet, at least we have not seen any 
related papers before. In this paper, we use the elementary method and the properties of the 
primitive roots to study the solvability of the congruent equation (1), and obtain its all prime 
solutions. That is, we shall prove the following conclusion: 


Theorem. Let n be a prime, then n satisfy the congruent equation 
jo —2) ptt) pack (1) 41 = 0 mod n 


if and only ifn = 2, 3 and 5. 
It is clear that our Theorem obtained all prime solutions of the congruent equation (1). 
About the general solutions of the congruent equation (1) is still an unsolved problem. 


2. Proof of the theorem 


In this section, we shall complete the proof of our Theorem directly. We only discuss the 
prime solutions of (1). 

(1) For n = 2, since 2 | 19@) + 1 = 2, so n = 2 is a prime solution of (1). 

(2) For n = 3, since 3 | 192) + 25) + 1 = 6, son =3 is a prime solution of (1). 

(3) For n = 5, since 5 | 15) + 25) 4 35 4 45@ +1 = 355, so n = 5 is also a prime 
solution of (1). 








(4) For prime n = p > 7, it is clear that p has at least a primitive root. Let g be a primitive 
root of p, that is to say, (g' —1, p) =1 for all 1 < i < p—2, the congruences 


g?-' =1mod p and g™®-) =1 mod p (2) 


hold for any positive integer m. 
Then from the properties of the primitive root mod p we know that g®°, g!, ---, g?~? isa 
reduced residue class. So we have the congruent equation 


(OCD) PD ose ge eee (3) 
= came ig) fageeee ® gP-2)'S@-1) 
= g®-8@-)) af 


dp. 
gen zy OCP 
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It is clear that for any prime p > 7, we have S(p—1) < p—2, so (gS@-D -1, p) = 1. Therefore, 
from (2) and (3) we have 


p-1)-S(p-1) _ 
gS@-) —I] 





( 
p81) 4 9S(m-I) 4 (n—- 1)8@-)) =/ = 0 mod p. 


So from this congruence we may immediately get 
fe Be aa eS aed i 


Thus, if prime p > 7, then it is not a solution of (1). This completes the proof of Theorem. 
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Abstract We introduce the notion of o-finite Ponomarev-system (f,M,X,{Pn}) to give a 
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X from some relatively compact metric space M. As applications, we systematically obtain 
characterizations on mssc-images of relatively compact metric spaces under certain covering- 
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§1. Introduction 


An investigation of relations between spaces with countable networks and images of sep- 
arable metric spaces is one of interesting questions on generalized metric spaces. Related to 
characterizing a space X having a certain countable network P by an image of a separable 
metric space M under some covering-mapping f, many results have been obtained. In the past, 
E. Michael [18] proved the following results. 

Theorem 1.1. ([{18], Proposition 10.2) The following are equivalent for a space X. 

(1) X is a cosmic space. 
(2) X is an image of a separable metric space. 

Theorem 1.2. ({18], Theorem 11.4) The following are equivalent for a space X. 

(1) X is an No-space. 
(2) X is a compact-covering image of a separable metric space. 

Also, he posed the following question, on page 999. 

Question 1.3. Find a better explanation for images of separable metric spaces? 
In [26], Y. Tanaka and Z. Li proved the following result. 

Theorem 1.4. ({26], Corollary 8 (1)) The following are equivalent for a space X, where 
“subsequence-covering” can be replaced by “pseudo-sequence-covering” . 


(1) X has a countable cs*-network (resp., cs-network, sn-network). 
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(2) X is a subsequence-covering (resp., sequence-covering, 1-sequence-covering) image of a sep- 
arable metric space. 

Recently, Y. Ge [8] sharpened Theorem 1.2 as follows. 

Theorem 1.5. ([8], Theorem 12) The following are equivalent for a space X. 

(1) X is an No-space. 
(2) X is a sequence-covering, compact-covering image of a separable metric space. 
(3) X is a sequentially-quotient image of a separable metric space. 

Taking these results into account, the following question naturally arises. 

Question 1.6. When P is a countable network (cs-network, sn-network, so-network) for 
X; how nice can the mapping f and the metric domain M be taken to? 

Around this question, the author of [2] has shown that, when P is a countable cs-network, 
f and M in Theorem 1.5 can be an mssc-mapping and a relatively compact metric space, 
respectively. More precisely, the following has been proved. 

Theorem 1.7. ([2], Theorem 2.5] The following are equivalent for a space X. 

(1) X is an No-space. 

(2) X is a sequence-covering, compact-covering mssc-image of a relatively compact metric 
space. 

(3) X is a sequentially-quotient image of a separable metric space. 

In the above result, f and M can not be any compact mapping and any compact metric 
space, respectively; see [2, Example 2.11 & Example 2.12]; and the key of the proof is to 
construct a sequence-covering, compact-covering mssc-mapping f from a relatively compact 
metric space M onto X. 

By the above, we are interested in the following question. 

Question 1.8. Find a consistent method to construct an mssc-mapping with covering- 
properties onto a space from some relatively compact metric space? 

In this paper, we answer Question 1.8 by introducing the notion of o-finite Ponomarev- 
system (f,M,X,{P,}) to give a consistent method to construct an mssc-mapping f with 
covering-properties onto a space X from some relatively compact metric space M. As applica- 
tions, we systematically obtain characterizations on mssc-images of relatively compact metric 
spaces under certain covering-mappings, which gives an answer for Question 1.6, particularly, 
for Question 1.3. 

Throughout this paper, all spaces are regular and 7,, N denotes the set of all natural 
numbers, w = NU {0}, and a convergent sequence includes its limit point. Let P be a family 
of subsets of X. Then UP, and ()P denote the union U{P : P € P}, and the intersection 
(\{P:P €P}, respectively. A sequence {z, :n € w} converging to xp is eventually in AC X, 
if {vp :n > no}U{ao} C A for some no € N, and it is frequently in Aif{r,, :k E N}Uf{ao} cA 
for some subsequence {2y, : k € N} of {a, :n € N}. 

Definition 1.9. Let P be a family of subsets of a space X. 

(1) For each x € X, P is a network at x in X, if x € ()P, and if x € U with U open in X, then 
there exists P € P such that re P CU. 

(2) P is a cs-network for X [10] if, for every convergent sequence S converging to « € U with 
U open in X, there exists P € P such that S is eventually in P C U. 


Vol.4 Notes on mssc-images of relatively compact metric spaces 111 





(3) P is a cs*-network for X [5] if, for every convergent sequence S' converging to « € U with 
U open in X, there exists P € P such that S is frequently in P C U. 

(4) P is a k-network of X [19] if, for every compact subset kK C U with U open in X, there 
exists a finite F C P such that K CUF CU. 

Definition 1.10. Let P = U{Pz : x € X} be a family of subsets of a space X satisfying 
that, for each x € X, P, is a network at x in X, and if U,V € Pz, then W C UNV for some 
W €P,. 

(1) P is a weak base for X [21], if G C X such that for each x € G, there exists P € Py 
satisfying P C G, then G is open in X. 

(2) P is an sn-network for X [14], if each member of P, is a sequential neighborhood of x in 
Xx. 

(3) P is an so-network for X [12], if each member of P, is sequentially open in X. 

(4) The above P,, is respectively a weak base, an sn-network, and an so-network at x in X [12]. 

Definition 1.11. ([4]) Let P be a subset of a space X. 
(1) P is a sequential neighborhood of « if, for every convergent sequence S converging to x in 
X, S is eventually in P. 
(2) P is a sequentially open subset of X if, for every « € P, P is a sequential neighborhood of 
x. 

Definition 1.12. Let X be a space. 
(1) X is relatively compact, if X is compact. 
(2) X is an No-space [18] (resp., sn-second countable space [7], so-second countable space, g- 
second countable space [21], second countable space [3]), if X has a countable cs-network (resp., 
countable sn-network, countable so-network, countable weak base, countable base). 
(3) X is a sequential space [4], if every sequentially open subset of X is open. 

Remark 1.13. ({15]) (1) For a space, weak base > sn-network => cs-network. 
(2) An sn-network for a sequential space is a weak base. 

Remark 1.14. (1) It is easy to see that “compact metric = relatively compact metric > 
separable metric”, and these implications can not be reversed from Example 2.14 and Example 
2.15. 

(2) It is well-known that a space X is an No-space if and only if X has a countable k-network 
(cs*-network), see [23, Proposition C], for example. 
Remark 1.15. Let f : X — Y be a mapping. 


(1) f is an mssc-mapping[13], if X is a subspace of the product space II X,, of a family 
neN 
{X, : n € N} of metric spaces, and for every y € Y, there exists a sequence {Vy :n € N} of 


open neighborhoods of y in Y such that each pn(f—!(Vyn)) is a compact subset of X,,, where 
Dn: [[% — X,, is the projection. 

(2) fis an 1-sequence-covering mapping [14] if, for every y € Y, there exists x, € f~'(y) 
such that whenever {y, : n € N} is a sequence converging to y in Y there exists a sequence 
{xy :n € N} converging to xy in X with each x, € f~'(yn). 

(3) f is a 2-sequence-covering mapping [14] if, for every y € Y, x € f(y), and sequence 
{yn : n € N} converging to y in Y, there exists a sequence {x,, : n € N} converging to x in X 
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with each x, € f~' (Yn). 

(4) f is an 1-sequentially quotient mapping [17] if, for every y € Y, there exists ry € f~'(y) 
such that whenever {y, : n € N} is a sequence converging to y in Y there exists a sequence 
{x, :k © N} converging to x, in X with each xz € f~'(Yn,)- 

(5) f is a sequence-covering mapping [20] if, for every convergent sequence S' in Y, there exists 
a convergent sequence L in X such that f(L) = S. 

(6) f is a compact-covering mapping [18] if, for every compact subset K of Y, there exists a 
compact subset L of X such that f(L) = K. 

(7) f is a pseudo-sequence-covering mapping [11], if for every convergent sequence S of Y, there 
exists a compact subset AK of X such that f(K) = S. 

(8) f is a sequentially-quotient mapping [1] if, for every convergent sequence S' in Y, there 
exists a convergent sequence L in X such that f(Z) is a subsequence of S. 


For terms are not defined here, please refer to [3, 24]. 


§2. Results 


Definition 2.1. Let P = {P, : n € N} be a countable network for a space X. Because 
X is T, and regular, we may assume that every member of P is closed. For each n € N, let 
Pn = {X}U{P, : i <n} = {Py : a € Ay}, where A, is a finite set, and let every A, be 
endowed with the discrete topology. Setting 


M = {a = (a;) € [| 4: : {Pa, : t © N} forms a network at some point EX}. 
ieN 


Then M, which is a subspace of the product space II A;, is a metric space. Since X is T; and 
regular, %q is unique for each a € M. We define f:M — X by f(a) = x, for eacha € M. 
Then f is an mssc-mapping from a relatively compact metric space onto X by the following 
Theorem 2.3. The system (f,M,X,{P,,}) is a o-finite Ponomarev-system. 

Remark 2.2. For more details on Ponomarev-systems and images of metric spaces, see 
(9,16,25], for example. 

Theorem 2.3. Let (f,M,X,{Pn}) be the system in Definition 2.1. Then the following 
hold. 
(1) M is a relatively compact metric space. 
(2) f is an mssc-mapping. 

Proof. (1) For every i € N, since A; is finite, A; is compact metric. Then II A; is compact 

ieN 

metric, and so Mc II A; is. It implies that M is relatively compact metric. 
(2) It suffices to es following facts (a), (b), and (c). 

(a) f is onto. 
Let « € X. Since P is a countable network for X, (P)z = {P € P : x € P} is a countable 
network at 2 in X. We may assume that (P), = {Pr : j © N}, where P,,; € Pj;) for some 
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i(j) € N satisfying i(j) < i(j +1). For each i € N, taking Py, = Pj;) if i = i(7) for some j € N, 
and otherwise, Py, = X. Then (a;) € [| 4 and {Py, : i € N} = (P), forms a network at x 
ieN 
in X. Setting a = (a;), we get a € M and f(a) =z. 
(b) f is continuous. 
Let x = f(a) € U with U open in X anda€é M. Setting a = (a;) € [| 4 where {P,, : i € N} 


forms a network at x in X. Then there exists n € N such ial a € Py, C U. Setting 
M, = {b = (8;) € M: By, = an}, then M, is an open neighborhood of a in M. For each 
b € Mg, we get f(b) € Ps, = Pa, CU. It implies that f(M.) CU. 

(c) f is an mssc-mapping. 

Let « € X. For each i € N, taking V,,, = X. Then {V,; : 7 © N} is a sequence of open 
neighborhoods of x in X. Since A, is finite, A; is compact. Then p;(f—!(Vz)) = pi(f7'(X)) C 
A; is compact. It implies that f is an mssc-mapping. 

Remark 2.4. In the view of (a) in the proof (1) = (2) of [2, Theorem 2.5], f is onto by 
choosing a = (a,) satisfying f(a) = x, where x € Py, € P;. Unfortunately, this argument is 
not true. In fact, we may pick « € Py, = X € P;. Then {Py, : i € N} = {X} does not form 
a network at « in X. It is a contradiction. The wrong argument is corrected as in the proof 
(2).(a) of Theorem 2.3. 

Theorem 2.5. Let (f,M,X,{Pn}) be a o-finite Ponomareyv-system. Then the following 





) 

) If P is a cs*-network for X, then f is pseudo-sequence-covering. 
4) If P is an sn-network for X, then f is 1-sequence-covering. 

) If P is an so-network for X, then f is 2-sequence-covering. 

Proof. (1) Let S = {a : m € w} be a convergent sequence converging to zo in X. 
Suppose that U is an open neighborhood of S in X. A family A of subsets of X has property 
cs(S,U) if: 

i) A is finite. 
ii) foreach PE A,OAPNSCPCU. 
iii) for each x, € S, there exists unique P,,, € A such that 2, € Pz,,. 
iv) if aj € PEA, then S — P is finite. 
For each i € N, since A = {X} C P; has property cs(S,X) and P; is finite, we can assume 
that 


= Se Ss 


{A cP; : A has property cs(S, X)} = {Aig 29 = mi-1 +1,..., 4}, 


where no = 0. By this notation, for each j € N, there is unique i(7) € N such that Aj(;) has 
property cs(S, X). Then for each j € N, we can put Ajj) = {Pa : a € Ej}, where E; is a finite 
subset of Aj. 
For each 7 € N, m € w and z,, € S, it follows from (iii) that there is unique ajm € E; 
such that t, € Po,,, € Aj. Let am = (Qjm) € [4 Cc [[ 4. We shall prove that 
jen jeN 


{P. : 7 © N} is a network at x,, in X. In fact, let x, € U with U open in X, we consider 


Aj,m 
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two following cases. 

(a) If m = 0, then S is eventually in P,, C U for some P,, € P. For each x € S— P,,, 
let c € P, C X —(S— {x}) for some P, € P. Then G = {P,,}U{P,: «2 ¢ S—P,,} CP 
has property cs(S,X). Since G is finite, G C P; for some i € N. It implies that G = Ajv;) 
for some i(j) € N with j € {mj_14+1,..., mj}. Since tm = Zo € Pa;o, Pajo = Pro. Hence 
Im = Xo € Pa, CU. 

(b) If m £0, then S — {x,,} is eventually in P,, C X — {%m} for some P,, € P. For each 
x € (S — {am}) — Pi, let « € P, C X — (S — {x}) for some P, € P, and let tp, € Py, C 
UN(X —(S'—{am})) for some P;,, € P. Then H = {P,,}U{Pz,, }U{ Pr : 2 € (S—{@m}) — Par, } 
has property cs(S,X). Since H is finite, H C P; for some i € N. It implies that H = Aj,;) 
for some i(j7) € N with 7 € {mat 1,...,ms}. Since am € Po, 45 Posen 
Im © Pajm CU. 

By the above, for each m € w we get dm = (Qj,m) € M satisfying f(a@m) = 2m. For each 


= P,,,- Hence 


j €N, since families H and G are finite, there exists mj; € N such that aj,m = aj,9 ifm > m,;. 
Hence the sequence {a;,m :m € N} converges to aj,9 in Aj. Thus, the sequence {a,, :m € N} 
converges to ag in M. Setting L = {am : m € w}, then L is a convergent sequence in M and 
f(L) = S. This shows that f is sequence-covering. 
(2)Let K be a compact subset of X. Suppose that V is an open neighborhood of K in X. A 
family B of subsets of X has property k(K,V) if: 

(i) B is finite. 

(ii) PAK #0 for each PE B. 

(iii) K CUBCYV. 

For each i € N, since B = {X} C P; has property k(K, X) and P; is finite, we can assume 
that 


{B Cc P; : B has property k(K,X)} = {By :j =m-14+1,..., ni}, 


where no = 0. By this notation, for each j € N, there is unique i(j) € N such that B,;) has 


property k(,X). Then for each 7 €¢ N, we can put By) = {Py : a € Fj}, where F; is ae 
subset of Aj. 

Setting L = {a = (aj) € [4 : (\(« A Py,) # 0}, we shall prove that L is a compact 
subset of VM satisfying f(L) _'K. fine 7 is compact-covering, by the following facts (a), (b), 
and (c). 

(a) L is compact. 

Since L C II F; and Il F;; is compact, we only need to prove that L is closed in Il F;. Let 


ieN ieN ieN 
a= (a;) € Il F,;—L. Then () (KO Pa,) = 0. Since KN Py, is compact, there exists ig € N such 
icN ieN 
that () (KN Py,) = 9. Setting W = {b = (8;) € [[“ : GB; = a; if i < ig}. Then W is an open 
i<io ieN 
neighborhood of a in II F, and WL = 9. If not, there exists b = (3;) € WNL. Since be L, 
ieN 


( (Kn Po,) #0, hence (| (KN Pg,) # 0. Since be W, (| (KN Pa.) = (| (KN Py,) #0. 


ieN i<io i<io i<io 
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This is a contradiction of the fact that () (KM Pa,;) = 0. 
i<io 
(b) LC Mand f(L) CK. 
Let a = (a;) € L, then a € [[“ and altsa Po,) #0. Taking x € ( (E20 Pa,). If 


ieN iceN icN 
{Po, : i € N} is a network at x in X, then a € M and f(a) = 2, hence LC M and f(L) C K. 


So we only need to prove that {P,, : i € N} is a network at 7 in X. Let V be an open 
neighborhood of x in X. There exist an open subset W of K such that « © W, and compact 
subsets W and K —W such that W" CV and K—W c X — {x}, where W” is the closure 
of W in K. Since P is a k-network for X, there exist finite families F; C P and Fo C P such 
that W CUFF, CV and K—W CUFo C X— {x}. We may assume that PAK # 0) for each 
Pe F,UFy. Setting F = F, UFo, then F has property k(K, X). It implies that F = Bj;) for 
some 2(7) € N with 7 © {mi +1,...,mi}. Since sm € Py, € Bi), Po, € Fy. Then 2 € Py CV. 
This prove that {P., : 7 € N} is a network at # in X. 

(c) Kc f(L). 
Let « € K. For each i € N, there exists a; € F; such that x € P,,. Setting a = (a;), then 
a € L. Furthermore, f(a) = « as in the proof of (b). So K Cc f(L). 
(3) Let S be a convergent sequence of X. By using [26, Lemma 3] and as in the proof of (2), 
here S plays the role of the compact subset K, there exists a compact subset L of M such that 
f(L) = S. Then, f is pseudo-sequence-covering. 
(4) Let P = U{Pz : « © X}, where each P, is an sn-network at x in X. Let x € X. We may 
assume that Pz = {Pz,j : 7 € N}, where P,,; € P,(;) with some i(7) € N satisfying (7) < i(j+1). 
For each i € N, take Py, = Pj) if i = i(j) for some j € N, and otherwise, Py, = X. Then 
(ai) € ][4. and {P,, : i © N} = P, forms a network at x in X. Setting a, = (a;), then 
ae Ne and a, € f—1(x). For each n EN, setting B, = {b = (6;) € M: B = a ifi < n}. 
Then {B,, : n € N} is a decreasing neighborhood base at a, in M. For each n € N, we get 
f(Bn) = (| Pa: In fact, for each b = ((;) € Bn, f(b) = (| Pa, € Po, = Pa, ifi <n. Then 


i<n 1EN 
f(b) € (| Po,. It implies that f(Bn) C (| Pa,- Conversely, if 2 € () Pa,, then x = f(b) for 


some b = (6;) € M. Setting c = (7%), where y; = a; ifi <n, and y% = Bn if i > n. Then 
{P,, : 7 © N} is a network at x in X. For each i > n, since Bin € Aj-n C Ai, y; € Aj. Then 
c € By and f(c) = «. It implies that (| Pa, C f(Bn). 
i<n 

Suppose that the sequence {x,; : 7 € N} converges to # in X. For each i EN, since Py, is 
a sequential neighborhood of x in X, there exists 7; € N such that x; € P,, for each j > jj. 
For each n €N, setting j(n) = max{j; : i < n}, then x; € () Po, = f(Bn) for each j > j(n). 

i<n 

It implies that f~'(2;) NB, 4 0 for each j > j(n). We may assume that 1 < j(n) < j(n +1). 
For each j € N, taking 


f-* (x3) if j < j(1), 


aj e 
fo\(aj)N By if j(n) <j <j(n+1). 


Then the sequence {a; : 7 € N} converges to a, in M. In fact, let U be a neighborhood of 
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a,, then there exists n € N such that a, € B, C U. Obviously, a; € B, C U for each j > j(n), 
so {a; : 7 © N} converges to a,. Note that f(a;) = x; for every 7 € N. It implies that f is 
1-sequence-covering. 

(5) Let « € X and a, € f~'(x). Setting az = (aj), then a; € A; for each i E N, and {Py, :i€ 
N} forms a network at x in X. For each n €N, setting B, = {b = (8;) € M: 8; =a; ifi <n}. 
Then {B, : n € N} is a decreasing neighborhood base at a; in M, and for each n € N, 
f(Bn) =Mi<n Pa; as in the proof of (4). 

Suppose that sequence {x; : 7 € N} converges to in X. For each i € N, since Py, is 
sequentially open, there exists j; € N such that x; € Py, for each j > j;. As in the proof (4), 
we get a sequence {a,; : 7 € N} converges to a, in M such that f(a;) = 2; for each j € N. It 
implies that f is 2-sequence-covering. 

By using Theorem 2.3 and Theorem 2.5 above, we obtain characterizations of mssc-images 
of relatively compact metric spaces under certain covering-mappings, which sharpen results in 
[18,26] and more. Firstly, we sharpen Theorem 1.1 as following. 

Corollary 2.6. The following are equivalent for a space X. 

(1) X is a cosmic space. 
(2) X is an mssc-image of a relatively compact metric space. 
(3) X is an image of a separable metric space. 

Proof. We only need to prove (1) = (2), other implications are routine. Since X is a cosmic 
space, X has a countable network P. Then the o-finite Ponomarev-system (f,M,X,{P,,}) 
exists. Therefore, f is an mssc-mapping and M is relatively compact metric by Theorem 2.3. 
It implies that X is an mssc-image of a relatively compact metric space. 

Remark 2.7. Corollary 2.6 is an answer of Question 1.3. 

Next, we extend partly Theorem 1.4 as follows. 

Corollary 2.8. The following are equivalent for a space X. 

(1) X is an sn-second countable space. 

(2) X is an 1-sequence-covering, compact-covering mssc-image of a relatively compact metric 
space. 

(3) X is an 1-sequentially-quotient image of a separable metric space. 

Proof. (1) = (2). Since X is an sn-second countable space, X has a countable closed 
sn-network P, and a countable closed k-network P2. Set P = P, UP2, then P is a countable sn- 
network and k-network for X. It follows from Definition 2.1 that the o-finite Ponomarev-system 
(f, M, X,{P,}) exists. Therefore, f is an 1-sequence-covering, compact-covering mssc-mapping 
from a relatively compact metric space M onto X by Theorem 2.3 and Theorem 2.5. 

(2) = (8). It is obvious. 

(3) = (1). Let f : M — X be an 1-sequentially-quotient mapping from a separable 
metric space M onto X. Since M is separable metric, M has a countable base 6B. We may 
assume that B is closed under finite intersections. For each 2 € X, there exists a, € f~'(x) 
such that whenever {x,, : n € N} is a sequence converging to x in X there exists a sequence 
{ax : k € N} converging to a, in M with each ax, € f~'(an,). Let B, = {B € B: az € B} and 
P. = f (Bz). We shall prove that P = U{Pz : x € X} is a countable sn-network for X. Since 
B is countable, P is countable. It suffices to prove the following facts (a), (b), and (c) for every 
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cEX. 

(a) Pz is a network at x in X. 

It is straightforward because 6, is a neighborhood base at a, in M. 

(b) If Pi, Po € Px, then P C P, M Pp for some P € Py. 

Let P, = f(Bi), Po = f(B2) with Bi, Bo € B,. Then B, NM By € B,. Setting P = 
f(Bi.N Bz), we get PE P, and Pc Pi NA Pp. 

(c) Each P € P,, is a sequential neighborhood of «. 

Let P = f(B) with B € By, and let {x, : n € N} be a sequence converging to x in 
X. For each subsequence {zp, : 7 € N} of {x, : n € N}, since f is 1-sequentially-quotient, 
there exists a sequence {a, : k € N} converging to az in M with each a, € f~'(an,,)- 
Since {a, : k € N} U {a,} is eventually in B, {x,, : i © N}U {2} is frequently in P. Then 
{ty :n © N}U {2} is frequently in P. It follows from [6, Remark 1.4] that P is a sequential 
neighborhood of a. 

Corollary 2.9. The following are equivalent for a space X. 

(1) X is a g-second countable space. 

(2) X is an 1-sequence-covering, compact-covering, quotient mssc-image of a relatively compact 
metric space. 

(3) X is an 1-sequentially-quotient, quotient image of a separable metric space. 

Proof. (1) = (2). It follows from Corollary 2.8 that X is an 1-sequence-covering mssc- 
image of a relatively compact metric space under some mapping f. Since f is an 1-sequence- 
covering mapping onto a g-second countable space, f is a quotient mapping by [15, Lemma 
3.5]. It implies that X is an 1-sequence-covering, quotient mssc-image of a relatively compact 
metric space. 

(2) = (8). It is obvious. 

(3) = (1). Since X is a quotient image of a separable metric space, X is a sequential space 
by [4, Proposition 1.2]. Then X is a sequential space having a countable sn-network P by 
Corollary 2.8. It follows from Lemma 1.13 that P is a weak base for X. Then X is a g-second 
countable space. 

In the following, we prove a result on spaces having countable so-networks. 

Corollary 2.10. The following are equivalent for a space X. 

(1) X is an so-second countable space. 

(2) X is a 2-sequence-covering, compact-covering mssc-image of a relatively compact metric 
space. 

(3) X is a 2-sequence-covering image of a separable metric space. 

Proof. (1) = (2). Since X is an so-second countable space, X has a countable closed 
so-network P; and a countable closed k-network P2. Set P = P| UP2, then P is a countable so- 
network and k-network for X. It follows from Definition 2.1 that the o-finite Ponomarev-system 
(f, M, X,{P,}) exists. Therefore, f is a 2-sequence-covering, compact-covering mssc-mapping 
from a relatively compact metric space M onto X by Theorem 2.3 and Theorem 2.5. 

(2) = (8). It is obvious. 

(3) => (1). Let f : M — X be a 2-sequence-covering mapping from a separable metric 
space M onto X. Since M is separable metric, M has a countable base 6. For each x € X, 
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let B, = {BE B: f-'(xz) NB F O}, and P, be the family of finite intersections of members 
of f(B,). We shall prove that P = U{P. : x € X} is a countable so-network for X. Since B 
is countable, P is countable. It suffices to prove the following facts (a), (b), and (c) for every 
cEX. 

(a) Pz is a network at x in X. 

It is straightforward because B, is a base of f~+(x) in M. 

(b) If P,, Po € Pz, then P C P, M Pp for some P € Py. 

It is straightforward by choosing P = P, M Pp». 

(c) Each P € P,, is sequentially open. 

For each y € f(B) with B € Bz, we get f-'(y) NBO. Let {yn : n € N} be a sequence 
converging to y in X. Then there exists a sequence {a,, :n € N} converging toa € f~!(y)N B 
such that f(a@n) = yn for each n € N. Since B is open, {an :n € N}U {a} is eventually in B. It 
implies that {y, :n € N}U {y} is eventually in f(B). Hence f(B) is a sequential neighborhood 
of y in X. Therefore, f(B) is sequentially open. Since each P € P, is some intersection of 
finitely many members of f(6,), P is sequentially open. 

Corollary 2.11. The following are equivalent for a space X. 

(1) X is a second countable space. 

(2) X is a 2-sequence-covering, compact-covering, quotient mssc-image of a relatively compact 
metric space. 

(3) X is a 2-sequence-covering, quotient image of a separable metric space. 

Proof. (1) = (2). It follows from Corollary 2.10 that X is a 2-sequence-covering mssc- 
image of a relatively compact metric space under some mapping f. Since f is a 2-sequence- 
covering mapping onto a second countable space, f is a quotient mapping by [15, Lemma 3.5]. 
It implies that X is a 2-sequence-covering, quotient mssc-image of a relatively compact metric 
space. 

(2) = (8). It is obvious. 

(3) > (1). Since X is a quotient image of a separable metric space, X is a sequential 
space by [4, Proposition 1.2]. Then X is a sequential space having a countable so-network P 
by Corollary 2.10. For each P € P, since P is sequentially open and X is sequential, P is open. 
Then P is a countable base for X, i.e., X is a second countable space. 

Remark 2.12. By using Theorem 2.3 and Theorem 2.5, we also get the proof of Theorem 
1.7 again, which was presented in [2]. 

From the above results, we obtain preservations of certain spaces under covering-mappings 
as follows. 

Corollary 2.13. Let f : X —> Y bea mapping. Then the following hold. 

(1) If X is an sn-second countable space and f is 1-sequentially-quotient, then Y is an sn-second 
countable space. 

(2) If X is a g-second countable space and f is 1-sequentially-quotient, quotient, then Y is a 
g-second countable space. 

(3) If X is an so-second countable space and f is 2-sequence-covering, then Y is an so-second 
countable space. 

(4) If X is a second countable space and f is 2-sequence-covering, quotient, then Y is a second 
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countable space. 

Proof. We only need to prove (1), by similar arguments, we get the other. Since X is an 
sn-second countable space, X is an 1-sequentially-quotient image of a separable metric space 
under some mapping g by Corollary 2.8. It is easy to see that fog is also an 1-sequentially- 
quotient mapping. Then Y is an 1-sequentially-quotient image of a separable metric space 
under the mapping f og. Therefore, Y is an sn-second countable space by Corollary 2.8. 

Finally, we give examples to illustrate the above results. 

Example 2.14. ({2], Example 2.8) A relatively compact metric space is not compact. 

Example 2.15. ([2], Example 2.9) A separable metric space is not relatively compact. 

Example 2.16. A 2-sequence-covering, compact-covering mapping from a separable met- 
ric space is not an mssc-mapping. 

Proof. Let f be the mapping in [2, Example 2.10]. Then f is also 2-sequence-covering. 
This complete the proof. 

Example 2.17. A second countable space is not any image of a compact metric space. It 
implies that “relatively compact metric” in the above results can not be replaced by “compact 
metric”. 

Proof. Let R be the set of all real numbers endowed with the usual topology. Then R is a 
second countable space. Since R is not compact, R is not any image of a compact metric space. 

Example 2.18. A g-second countable space is not any 1-sequence-covering, compact- 
covering compact image of a separable metric space. It implies that “mssc-image” in Corollary 
2.8 can not be replaced by “compact image”. 

Proof. See [22, Example 2.14(3)]. 
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